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This thesis deals with the dynamics of irreversible processes within the context of the 
general theory of relativity. In particular, we address the problem of the 'infinite' speed 
of propagation of thermal disturbances in a dissipative fluid. Although this problem is 
not new, its best known solution - the Israel and Stewart second order expansion - has an 
effective, rather than fundamental, character. 

The present work builds on the multi-fluid variational approach to relativistic dissipation, 
pioneered by Carter, and provides a dynamical theory of heat conduction. The novel 
property of such approach is the thermodynamic interpretation associated with a two-fluid 
system whose constituents are matter and entropy. The dynamics of this model leads to 
a relativistic generalisation of the Cattaneo equation; the constitutive relation for causal 
heat transport. A comparison with the Israel and Stewart model is presented and its 
equivalence is shown. This discussion provides new insights into the not-well understood 
definition of a non-equilibrium temperature. 

A crucial feature of the multi-fluid approach is the interaction between its constituents. It 
is a well known fact that when two, or more, fluids interact, instabilities may occur. Within 
this work, the two-stream instability analysis is extended to the relativistic domain. As 
far as the author is aware, such extension has not been discussed in the literature. The 
analysis allows to assess the stability and causality of relativistic models of matter and 
their linear deviations from thermodynamic equilibrium directly from their equations of 
state or, equivalently, their Lagrangian densities. For completeness, a brief digression on 
a consistent (stable and causal) 'first-order' model is also included. 

Finally, the road to follow is laid by posing some physical applications together with some 
future perspectives and closing remarks. 

To sum up, the variational approach to heat conduction presented in this thesis consti- 
tutes a mathematically promising formalism to explore the relativistic evolution towards 
equilibrium of dissipative fiuids in a dynamical manner and to get a deeper conceptual 
understanding of non-equilibrium thermodynamic quantities. Moreover, it might also be 
useful to explore the more fundamental issues of the irreversible dynamics of relativity and 
its connections with the time asymmetry of nature. 
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Introduction 



It has become fashionable to justify one's scientific endeavours by means of 
their immediate impact in shaping modern society in a manner that seems 
to accommodate research into two categories: the self-evident applicable to 
technological development kind and the interesting but 'why do we care?' type, 
to which most 'curiosity-driven' research belongs. The work presented in this 
thesis is in the latter group. Therefore, this chapter is written for those readers 
whose interest in this work was driven by the mere curiosity of knowing what 
is this thesis about? but whose main concern really is why they should care?. 

Does a moving body appear cold? This remarkably simple question, raised by the late 
Professor P. T. Landsberg some forty years ago, highlights a profound missing link in our 
current understanding of classical thermodynamics and relativity. Some might interpret 
this as the lack of a relativistic transformation law for temperature measurements (see 
|45l Hll I43|). Some others, as we do here, believe that this is a misunderstood concern 
about the nature of thermal motion in a relativistic settings. Such a view shifts the 
attention from a mere exercise in velocity transformations to an argument about thermal 
equilibrium between moving bodies and the inertial properties of heat |37| , l26 t[TFlll4| . Thus, 
it promotes the discussion about the missing Lorentz transformation for temperature to the 
more elaborate - and physically meaningful - problem of heat exchange and irreversibility 
in a relativistic context. 

Irreversible processes in general relativity represent an outstanding challenge in our 
understanding of the macroscopic dynamics of a wide range of physical systems. Although 
the problem of relativistic dissipation found a pragmatic solution through the Israel & 
Stewart formalism (c./. Chapter 4), in recent years there has been a resurgence of interest 
in the topic. This thesis is motivated by phenomena whose leading mechanism of dissipation 
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is due to heat flow. Such is the case of the dynamics of super-fluid neutron stars, where 
thermal and dissipative effects are expected to impact on observations; or high energy 
gases, with photons providing the dominant pressure contribution. 

1.1 The empirical nature of Thermodynamics 

Progress in each individual aspect of science is normally preceded by a revolution. The 
twentieth century is the closest and most vivid example of this, not only because the early 
years saw the rise of Relativity and Quantum Theory - our two most prominent scientific 
achievements - but also due to the urgent need to reconcile the problems they pose in our 
quest of a unified vision of the universe. This manuscript is concerned with the motion 
of systems whose macroscopic dynamics follow the rules prescribed by the general theory 
of relativity, and where the interactions between their microscopic components become 
apparent only through a process of averaging, i.e. through a thermodynamic description. 

The central concept in thermodynamics is that of thermal equilibrium. We will define 
precisely what it means for a system to be in such a state. For the time being it will suffice 
to say that it corresponds to the highly idealised (macroscopic) state in which an observer 
would not be able to distinguish past from future solely from the dynamics of the system, 
i.e. the system's motion is time symmetric. The vast majority of physical situations 
we experience occur away from thermal equilibrium e.g. the Earth's climate, the Sun's 
internal reactions, black hole formation, the expansion of the universe, life itself... That 
a system's evolution occurs in such a manner that its final state is thermal equilibrium 
follows from a set of empirical laws; it is not a consequence of the dynamics. In this sense, 
thermodynamics cannot be considered as a 'fundamental' theory, nevertheless its scope is 
universal. 

What is the role thermodynamics plays in a theory of macroscopic dynamics? The 
Laws of Thermodynamics are statements about nature which stem from the observation 
that certain phenomena - although allowed by the available theory of motion - simply do 
not occur. We need to impose these laws 'on top' of our more fundamental dynamical 
basis. Thus, the correct way to understand the role of thermodynamics is by regarding its 
laws as auxiliary hypotheses which rule out entire classes of dynamical processes. 

If we have found a fundamental theory whose elegance and scope satisfy our needs of 
beaut}Q, supplying the auxiliary thermodynamic hypothesis cannot do anything but to 
make our theory stronger. To illustrate this point, let us suppose that we have a theory 
which allows us to describe precisely the motions of the microscopic constituents of a given 



^This is indeed a very subjective notion. For the present purpose of discourse it will suffice to consider 
any theory with a non-empty physical content and a minimal set of hypotheses. 



1.1. The empirical nature of Thermodynamics Introduction 

material such that the total energy and linear and angular momentum are conserved. A 
collective measure of such motions corresponds to the temperature of the material. What 
are the possible dynamical outcomes of bringing into contact two pieces of the same material 
at different temperature ? We cannot decide the answer to this kind of question solely from 
the consequences of our fundamental theory, i.e. its conservation laws. We need more 
hypotheses. For example, if we call heat the energy and momentum transferred from one 
body to the other, we can introduce the following observation: 



(Second Law of Thermodynamics) In every natural process involving two 
bodies at different temperatures heat never flows from the colder to the warmer. 

The addition of this new statement to our fundamental theory rules out every dynamical 
possibility which fails this empirical criterion. Indeed, this is a part of what we observe as a 
universal fact. However, we cannot yet decide a unique outcome from the set of possibilities 
allowed by this 'restricted' dynamics. A less obvious hypothesis is the one which asserts 
that: 

(First La'w of Thermodynamics) There is no perpetual motion machine. 

This new hypothesis prevents our system from transferring heat from one body to the 
other in a periodic manner, analogous to a harmonic oscillator. This, however, still does 
not give us conclusive information about any 'final' state. A closer look at the question 
we originally asked shows that we have made a tacit assumption, that we know when two 
bodies are at different temperatures. Thus, in order to clarify such a subtlety, we introduce 
one more hypothesis 

(Zero-th Law of Thermodynamics) If two bodies are brought into thermal 
contact and there is no heat transfer, they are at the same temperature. 

From these three additional fact-like laws, we can conclude that whenever two bodies at 
different temperatures are brought into thermal contact the warmer will transfer energy 
to the colder. If the size of the bodies is 'small enough', there may be 'thermal bounces' 
from their boundaries, however those will be damped for otherwise we would be able to 
construct a perpetual motion machine. Finally, their temperatures will equilibrate when 
no further heat transfer occurs. The auxiliary hypotheses do no tell us how the process of 
temperature equilibration occurs. All the processes involved are driven by the dynamics 
of our original theory. These auxiliary hypotheses, the Laws of Thermodynamics, simply 
rule out every solution that does not seem to correspond to an observed fact. 
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We may wonder if this process of adding new hypotheses wiU ever come to an end. 
Perhaps we have been sloppy while interpreting the predictions of or our theory, or maybe 
it was just not good enough to reproduce the observed reality. These are serious concerns 
for which every proposed solution would lie far beyond the scope of this treatise. The best 
we can do is to keep a critical eye in every conclusion we can reach from the principles of our 
proposed theories. We are granted by experience that the three laws of thermodynamics 
we have introduced are universal and do not contradict each other. Should we find an 
internal inconsistency while incorporating these hypotheses into our theoretical construct 
we cannot do anything else but to withdraw our theor5Q If on the contrary, the addition of 
these laws to our dynamical principles poses no contradiction, we have made our dynamics 
less ambiguous; our theory is more complete. 

What is the empirical basis of this work? 

In this thesis, we use general relativity as our fundamental theory of motion. In the next 
chapter we will present the principles upon which the general theory of relativity is based 
together with its implications as a dynamical theory. Here, we should answer what do we 
mean by a fundamental theory? Throughout the entire manuscript, whenever we claim 
that a theory is fundamental, it means that there is a variational principle and a minimal 
set of empirical hypothesis from which the dynamics of the theory can be obtained. In 
this sense, general relativity does qualify as a fundamental theory whose basic hypotheses 
consist of the principle of equivalence and general covariance. A detailed introduction to the 
empirical basis of general relativity is presented in the next chapter and the consequences 
of adding the thermodynamic auxiliary hypotheses to the material sector of the theory is 
the central theme of the present work. 

The marriage of thermodynamics and general relativity is a very active debate. Here 
we address the problem of heat in the context of the relativistic dynamics of a non-perfect 
fluid. To this end, we will invoke the principle of equivalence to justify the study of the 
inertial properties of heat. In this sense, the dynamics of heat will be traced down to the 
variational ingredients used to derive the specific matter model we use. We will explain 
the precise meaning of this words in due time. First, let us verify that there is indeed an 
problem of heat and relativity. 



^In the words of Eddington: 

'If someone points out to you that your pet theory of the universe is in disagreement with 
MaxweU's equations - then so much the worse for Maxwell's equations. If it is found to be 
contradicted by observation - well these experimentalists do bungle things sometimes. But if 
your theory is found to be against the second law of thermodynamics I can give you no hope; 
there is nothing for it but to collapse in deepest humiliation.' 



1.2. Historical problems of relativistic heat conduction Introduction 

1.2 Historical problems of relativistic heat conduction 

Let us begin our discussion of heat transfer from a non-relativistic point of view. Consider 
some metallic object or any other good thermal conductor. How long will it take to feel 
the heat propagated if we were holding one of the ends? According to Fourier's theory 
of heat - the standard transport theory used in most 'practical' applications - the effect 
would be felt instantaneously all across the object. This is due to the parabolic nature of 
the partial differential equation describing the heat transport within the body. 

The heat equation is obtained by assuming the validity of the first and second laws of 
thermodynamics together with Fourier's law relating heat with temperature gradients in 
an un-relaxed manner. This implies that when two bodies at different temperatures are 
put in thermal contact, heat spontaneously flows from the warmer to the colder without 
any delay. Such conclusion cannot be satisfactory. One would expect the speed of thermal 
disturbances to be bound by the internal structure of the media they travel on. Cattaneo 
saw a way around this problem by introducing, in a reasonable (but arbitrary) manner, a 
modification to Fourier's law which takes into account the characteristic time a material 
takes to react to thermal stimuli [T8l 



q=-KVr -)■ Tq+q=-KVr. (1.1) 

Here q denotes the heat flux, T the temperature, k the thermal conductivity, r represents 
the relaxation time of the medium and the dot denotes time differentiation. Such an 
amendment leads to the telegrapher equation for the propagation of heat signals, with a 
finite bound on the speed they can reach. It is worth mentioning that Fourier's law is 
the simplest, but not the most general, possible ansatz to ensure that the second law is 
satisfied. This is done by explicitly making the change in entropy a quadratic function of 
the heat flow. We should keep this in mind in the forthcoming discussion. 

In the relativistic case, the unbounded speed of thermal disturbances implied by the 
parabolic nature of the heat equation is more than mere inconvenience, it is indeed a 
fundamental problem. To give context to our discussion, let us consider the case where 
matter's motion is represented by a fluid whose normalised four- velocity iq^ given by n". 
The first attempt of a relativistic extension of the heat equation was introduced by Eckart 
[25j . The kind of model he proposed has become a stereotype of a class of theories referred 
as first-order relativistic theories of dissipation. In this class, the entropy current of the 
model, denoted by a vector field s"" which is generally not aligned with the matter four- 
velocity, may only depend on terms which are linear in deviations from equilibrium, namely 
the heat fiow or the shear viscosity. For this class of theories, the simplest way to impose 



We will elaborate on the choice for such four-velocity in Chapters 3 and 4. 
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the second law, which locahy takes the form s".„ > 0, leads to a relativistic version of 
Fourier's law 

q- = -nh'^b [Tfc + Tub] , (1.2) 

where q"" represents the heat flux, h is a projector orthogonal to the matter flow and the 
semi-colon and dot denote covariant and proper time differentiation, respectively. Being 



essentially identical to (1.1), equation (1.2) inevitably produces a non-causal theory. Fur- 
thermore, as shown by Hiscock and Lindblom [35], it also suffers from stability problems, 
i.e. some thermal disturbances may grow without bound. However, we also note that Eck- 
art's proposal exhibits an extra piece of information which is missing in the non-relativistic 
treatment; the acceleration term. This purely relativistic effect can be interpreted as being 
due to an effective 'mass' per unit entropy given precisely by the temperature (see discus- 
sions in |26l [TBI HSj [6] ) . It is worth mentioning that this term has recently been suggested 
to be the origin of the afore mentioned instabilities (see j31j). However, here we adopt the 



view that the appearance of the four acceleration in (1.2) is an inevitable, and physically 
meaningful, feature of any relativistic theory of dissipation. 

The failure of first order theories to produce a theory of heat conduction compatible 
with the principle of causality can be tracked down to their definition of the entropy 
current. In an effort driven by simplicity, it is not permissible to prematurely drop higher 
order terms in deviations from equilibrium, since they may give rise to linear terms after 
the differentiation required by the second law. This point of view is at the heart of the 
class of second order theories of heat conduction whose key contribution is the widely 
known Israel & Stewart model. Here, the entropy current used by Eckart is extended to 
include all the possible second order combinations of dissipative effects [TQl [371 ES] . This 
strategy, analogous to the Grad's 14-moment method, is firmly ground on kinetic theory 
and provides a causal and stable account of relativistic dissipation. It is not our intention 
here to further explain this particular framework, but to note that, in spite of its success 
regarding stability and causality tests, the price to pay is the introduction of a set of second 
order couplings that, in principle, can be measured but which cannot be obtained within 
the realm of the theory. A detailed description of the Israel and Stewart formalism is given 
in Chapter 4. 

Why do we care? 

Owing to an increasing interest in hydrodynamic descriptions of high-energy relativistic 
plasmas, attention to relativistic theories of dissipation of the Israel & Stewart type has 
been renewed. In many such applications, whose main source of dissipation is due to 
viscosity, heat can effectively be considered as a secondary effect. Indeed, it was this kind 



6 



1.3. The uaulti-Ruid approach to thermal dynamics Introduction 

scenario, in which the ratio of viscous to thermal resistivity can be pushed down to order 
unity, that was the main source of motivation for the development of Israel &: Stewart 
second order theory. Our interests, however, follow from quite a different motivation. It 
ranges from recent efforts to model the dynamics of super-fluid neutron stars [21 [1] to high 
energy gases with photons providing the dominant pressure contribution [16]. 

1.3 The multi-fluid approach to thermal dynamics 

In this brief account of a long standing problem in relativity, it is our wish to shed some 
light upon the more fundamentally satisfactory variational approach to relativistic heat 
conduction |161 [T31 [H]- Historically, due to a simplistic view of the components of the 
theory, this model, pioneered by Carter, failed the tests of stability and causality. However, 
it was later shown by Prion [65] that the predictions of Carter's complete theory are 
equivalent to, and at second order physically indistinguishable from, those of Israel &: 
Stewart. 

Perhaps the most attractive feature of a variational construction of relativistic heat 
conduction is that, once the equation of state (or, equivalently, the Lagrangian density) of 
the system is known, the theory contains no free parameters. The price for this however 
is the inclusion of transport quantities in the "Lagrangian" density of the system. This is 
nevertheless inevitable in any circumstance where one's aspiration is to describe, at least 
to linear accuracy, situations departing from local thermal equilibrium. 

In Carter's variational approach to heat conduction, one considers a multi-fluid system 
whose species are represented by a particle number density current n"" and an entropy flux 
s°. These two currents, together with the spacetime metric, constitute the fundamental 
fields of the matter sector for the Einstein-Hilbert action. General covariance requires the 
Lagrangian to be a proper scalar, therefore, it should depend only on covariant combina- 
tions of its fundamental fields. If we consider the metric as a passive field, the Lagrangian 
density can only depend on combinations of the two fiuxes, which includes the relative flow 
between them. This is precisely what we mean by the inclusion of transport quantities in 
the Lagrangian density. 

A constrained variation of the matter action, whose Lagrangian density has the charac- 
teristics described in the preceding paragraph and which is constrained by the conservation 
law of the particle number density flux, allows us to write the local conservation of en- 
ergy and momentum as a "force balance" equation. It is worth noting that each of the 
individual 'forces' appearing in such balance, takes a form completely analogous to the 
Lorentz force for electromagnetism in a general relativistic setting (c./. Section 5.2). The 
requirement for the local energy conservation law to follow as a Noether identity of the 
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variational principle, only needs one of the currents to be strictly conserved, n'^.^ = say. 
This allows an "extra" freedom to allocate the second law of thermodynamics by taking 
into account the possibility in which the production term s":^ takes values different from 
zero. 

One of the central problems to be faced before giving a real thermodynamic inter- 
pretation to the variational construction, lies in the correct interpretation of temperature 
measurements. As stressed in the opening sentence, this is a non-trivial problem in re- 
lativity and, therefore, a choice of frame is forced upon us. In this simple case, where we 
have only one conserved current, the choice is somewhat natural. In the spirit of phys- 
ical interpretation, we can chose to equip each infinitesimal part of the matter fluid with 
a thermometer and consider all physical measurements to be taken in, and with respect 
to, the matter frame. In such a case one can show that the projection of the canonical 
conjugate momentum to the entropy flux in the matter frame, which we denote here by 
6", corresponds to the thermodynamic temperature in the Gibbs sensqj Having clarified 
this, it is not difficult to show that the equation for the heat flux relative to the matter 
frame takes the form |48j 

f [r + u^'^c] + g" = -K/i"^ [ol + Ohb] . (1.3) 

Here all the quantities, with the exception of f and k which correspond to the effective 
relaxation time and thermal conductivity, have been introduced earlier. This result, which 
is a relativistic generalization of Cattaneo's equation [the modified Fourier's law ( |1.1[ )], 
follows directly from the variational principle together with the simplest assumption (in 
the sense previously discussed) to make the entropy production satisfy the second law of 
thermodynamics . 



The structure of (1.3) combines all the features present of both (1.1) and (1.2), includ- 
ing the acceleration term which in the variational context arises as a consequence of the 
equations of motion. The "check" marks indicate that these quantities depend on devi- 
ations from thermal equilibrium which are higher than second order. A comparison with 
an analogous expression obtained from the Israel & Stewart theory is beyond the scope of 
present means of experimental verification. 

Is this a covariant theory of the caloric? 

The answer to such a question is of ontological, rather than physical, character. As a 
reminder, the theory of caloric was the standard explanation of the mechanism of heat 

■'The rate of change of energy with respect to entropy with all other independent thermodynamic quant- 
ities fixed. 
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conduction until mid-nineteenth century. It posits the existence of a fluid-Uke substance - 
the caloric - which is self-repulsive and attracted to other forms of matter. In this view, 
temperature is a measure of the caloric density of a body. Thus, the previous observation 
that heat flows from warmer to colder bodies is a mere consequence of the natural tendency 
of a fluid to flow from from high to low density regions. It is now known that there is no 
such a substance, that the temperature of a material corresponds to the kinetic energy 
of its microscopic constituents and that it is overwhelmingly unlikely that the transfer 
of this kinetic energy occurs from 'colder' to 'warmer' bodies. Nevertheless, the caloric 
theory still provides accurate descriptions of thermal processes, provided it is applied in 
the right context. This is nothing new in the practice of science, it is simply the regular 
process of recovering older theories from newer ones in special limits. For example, in many 
practical calculations a Newtonian treatment would give quite a satisfactory description of 
the planets' orbits as if there was a gravitational pull from the Sun or the tides of ocean 
as if the Moon would exert a force on the water. This is better explained by Norton in the 
following quote: 

...In returning the older theories, the [reduction] relations revive a defunct 
ontology. More precisely, they do not show that heat is a fluid, or gravity a 
force; rather they show that in the right domain the world behaves as if they 
were. ..[551. 



While reading through the forthcoming chapters, it may become quite a tempting 
thought that entropy can be considered as a 'real' fluid. Being an ingredient of a variational 
construction, it is indeed no less real than the matter density current n" or the spacetime 
metric gab- In this sense, it is clear that the multi-fluid approach to heat conduction is 
an effective model to deal with thermal disturbances under specific circumstances; those 
in which entropy behaves as if it were a fluid. For some of us, this may seem a rather 
un-natural feature of entropy, even though we may be quite happy to use a fluid description 
for the collective motion of a set of particles. This is because there is a conceptual barrier 
associated with the 'kinetic meaning' of entropy (c./. section 3.3) and the properties of 
a fluid. However, soon we will show that it is, in fact, kinetic theory the one responsible 
for the fluid interpretation of entropjjj used in our variational formalism. Although the 
caloric theory was once considered to be equivalent to the kinetic theory of gases, it was 
eventually overthrown. Kinetic theory became the foundation of modern thermodynamics 
and the caloric has been regarded as a dead theory ever since. 



^The modern idea of entropy as "a measure of a system's disorder" makes it hard to visualise a situation 
in which such a disorder would have fluid properties. Remarkably, the thermodynamic limit of relativistic 
kinetic theory includes an entropy current, whose description can be given in terms of an 'entropic' four- 
velocity. We will address this point towards the end of Chapter 3. 
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In the multi-fluid approach to heat conduction, which shah be developed in chapter 
5, the relative flow between the matter and entropy currents plays a central role in the 
dynamical setting of a system away from thermal equilibrium. Moreover, it is through the 
inertial properties of heat that the entropy current can carry momentum from one part 
of the system to another, driving the dynamics towards equilibrium as dictated by our 
empirical laws. Therefore, in our view, treating the entropy density current as a 'real' 
fluid in the matter Lagrangian allows us to use the principle of equivalence to associate 
some inertial properties to the heat flux. Hence, allowing a hyperbolic description for the 
propagation of thermal perturbations. 

In conclusion, one could say that, physically, this is not a modern, covariant version of 
the caloric theory. However, it is one that certainly shares its main ontological properties. 

Past objections and our solution 

Let us emphasize that the main objection to Carter's approach was based on issues about 
stability and causality. This is a consequence of the simplistic spirit of the original model, 
ignoring a crucial effect present in almost every multi-fluid system; entrainment. This 
effect tilts the momenta due to the coupling between their associated currents. 

Compared to the standard single-fluid analysis, a multi-fluid system has more degrees 
of freedom; we need to account for relative flows between them. When two or more fluids 
are allowed to interact, instabilities may occur. Such instabilities are known to exist for a 
variety of configurations. For example, in shearing motion at an interface, this corresponds 
to the well known Kelvin-Helmholtz instability. One of the central results of this work, 
corresponds to a relativistic generalisation of the two-stream instability analysis. This is 
a generic phenomenon that does not require particular fine-tuning to be triggered, nor is 
limited to any specific physical system. The only requirements are that there must be a 
relative (background) fiow and a coupling between the fluids. 

It is worth noting that thermodynamic instabilities are the source of pattern and struc- 
ture formation and that our analysis provides a stability criteria directly from the pre- 
scribed (interacting) Lagrangian density or equation of state. The dynamical role of this 
instability has not been explored beyond the linear regime. 

1.4 The state of the art 

This thesis originally had a simple goal: to understand heat in a wide range of physical 
situations. To this end, we needed to provide a satisfactory account of dissipative processes 
in general relativity. This will remain a challenge with many interesting directions emerging 
from the work presented here. 
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So far, we have only scratched the surface of what is in fact a very difficult problem. 
In the following Chapters, the reader will find a summary of the most significant efforts to 
obtain a physically sound theory of heat conduction compatible with the tenets of relativity. 
A point to note is that in spite of the recent stir made on first order theories, in our view, 
physical indicators such as the presence of second sound in materials hint that the correct 
physics to model departures from thermal equilibrium cannot be less than second order. 

Although the Israel &: Stewart model is the most prominent and widely used tool to 
describe dissipative systems, the additional couplings - necessarily introduced in their ex- 
pansion - give the theory an effective, rather than fundamental, character. In this sense, 
the variational formalism not only provides us with an alternative to the Israel & Stew- 
art model, it centres the attention in the dynamical actors of any canonical theory; the 
canonical conjugate momenta. This observation makes the multi-fluid construction a very 
powerful tool which may help us tackle deeper problems in non-equilibrium thermodynam- 
ics of relativistic systems. 
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The question we address in this thesis is whether we can find a formulation of non- 
equihbriuni thermodynamics compatible with the tenets of general relativity from a fun- 
damental point of view. In particular, we want to study the case of heat conduction for 
matter models described by fluid dynamics. As highlighted in the introduction, this has 
been a long standing problem, partially solved by Israel and Stewart some thirty years ago 
|37] . It is one of our aims to explain the philosophy and achievements of their theory as a 
standard to compare our developments from a variational principle, continuing the work 
started by Carter [TilfTH] . 

We start this chapter with a primer of general relativity. Although the reader may 
find these derivations in most of the references cited in this work, we try to make explicit 
many assumptions and calculations often omitted in a general reading on the matter. 
This will also serve to introduce the conventions and notation used in the forthcoming 
chapters. The second part of this chapter is devoted to irreversible processes. It is a brief 
exposition of some of the developments in non-equilibrium thermodynamics relevant to 
this work, therefore many of the ideas will not be presented with the detail they deserve. 
The interested reader may wish to have a look at the following references |5H \38[ [50] . 

2.1 Spacetime and matter 

The aim of this section is to obtain the dynamical equations which govern the interplay 
between matter and spacetime. Although the derivation presented here is now standard 
knowledge, the construction and principles upon which relativity is based and the precise 
meaning of them is still an active debate among relativists [53]. Therefore, it would be 
negligent to leave aside a discussion of these matters. 
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2.1.1 The principles of relativity 

The elegance and power of any theory, not necessarily physical, can be judged objectively 
by the siniplicit}!^ of its assumptions and the scope of its validit}Q In this way, Einstein's 
theory of relativity is, to the author's taste, the most remarkable achievement in our efforts 
to understand the most familiar state of nature: motion. 

Up to the present date, general relativity has stood the most rigorous experimental 
tests. An excellent account of this was given by Clifford Will in |76j . and a more recent 
account can be found in [75] . The elegance of general relativity, however, can be questioned 
on the grounds of its principles. It is not an easy task to write a version of them which is 
not contentious. Strong arguments and objections have been put forward for over almost a 
century regarding their true physical content. It is not the intention here to participate in 
such debate, but to make a clear exposition of the tenets of relativity in a manner relevant 
to our attempts to reconcile thermodynamics with them. 

Einstein built his special theory of relativity on the empirical evidence that the speed of 
light in vacuum has the same constant value for all inertial observers - those moving with 
uniform relative velocity between each other - and the mathematical principle of Lorentz 
covariance, which states that any physical law should attain the same mathematical form 
in every frame related to another by a Lorentz transformation of coordinates, that is, 
independent of the inertial state of motion of the observer. This is known as the principle 
of relativity. 

It is tempting to enlarge the group of transformations which keep the covariance of 
physical laws to include a larger class of possible observers, such as those moving in an 
uniformly accelerated manner. This is the programme which gave rise to Einstein's theory 
of gravitation. The connection between accelerated motion and gravity is highly non-trivial 
and represents the fundamental physical assumption of the general theory of relativity. Just 
like the invariance of the speed of light for all observers, the equivalence between inertial 
and gravitational mass is an empirical fact which was promoted to a general principle, the 
principle of equivalence. 

The information encoded in the equivalence principle is purely physical. However, it 
has a far reaching consequence in Einstein's programme: it suggests a possible connection 
between the covariance of the physical laws for accelerated observers and gravity. 

The extension of the relativity principle requires physically meaningful statements to 
be independent of the general state of motion of the observer, not just inertial as in the case 



^Simplicity is a highly non-trivial concept within the philosophy of science. Here we follow the convention 
by Popper [B2] which states that a the degree of simplicity of a theory is in correspondence with its 
susceptibility to be 'easily' falsified. 

^Its universality. 
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of special relativity. That is, the laws of physics should preserve their form under a general 
coordinate transformation. This is a formulation of the principle of general covariance. It is 
a statement of mathematical nature whose connection with gravity through the equivalence 
principle form the core of the general theory of relativity. 

It has been argued that every physical theorjjj can be made to satisfy the principle of 
general covariance, and therefore it has been objected that its physical content is empty 
[53j . One should be aware of how this conclusion was reached - the fact that intrinsically 
there is no physics in it. It is together with the equivalence principle that we can learn 
what general covariance is for relativity: 

... it is the language for describing a world without distinction between the 
spacetime entities and the dynamical entities... [66] 

Dynamics and geometry 

To understand better the interplay between the mathematical language of general covari- 
ance and the physical principle of equivalence, it is worth discussing the simpler case of 
the electromagnetic field and Lorentz transformations. This exercise also makes explicit 
the underlying physical path to the general theory. 

The road to build the special theory of relativity combines the experimental results on 
the non-existence of an absolute reference framaj (which motivated the formulation of the 
relativity principle), together with the mathematical fact that Maxwell's equations do not 
change their form under Lorentz transformations, those which leave the Minkowski metric 
invariant. 

The combination of Lorentz covariance with the relativity principle leads to the physical 
notion of spacetime. In special relativity spacetime is fixed and plays no role in the dynam- 
ics of physical processes. This is a consequence of the fact that Lorentz transformations 
only relate inertial observers, in this sense we consider them as defining the observer kin- 
ematics of relativity. Thus, Einstein's programme is that of finding the observer dynamics 
which, through the equivalence principle, is coupled with the dynamics of spacetime. 



''We consider a theory as physical if it satisfies the empirical criteria prescribed by the underlying 
empirical laws. If in addition it can be derived from a variational principle we regard such a theory as 
fundamental. 

''The rigid aether assumption is a consequence of the finite value of the speed of light and the preferred 
Newtonian frame it singles out. This introduced the notion of conformal structure which is defined, even 
in Newtonian theory, by a field of light cones. See section 1.3 in |26j . 
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2.1.2 The symmetry of general relativity 

Now we proceed to understand how the symmetry imphed by general covariance merges 
with the principle of equivalence to produce the dynamics of matter and geometry. The 
construction presented here follows that of [HI [52l [201 133] • 

The mathematical model for space-time we use consists of a four-dimensional pseudo- 
Riemannian manifold M, whose metric g has the signature {— , -|-, -|-, -|-}. The aim should 
be then to obtain a set of equations for the components of the metric, gab, which in the non- 
relativistic limit reproduce the Newtonian theory of gravity as described by the Poisson 
equation, a linear second order partial differential equation for the Newtonian potential. 

For the gravitational theory, Hilbert proposed the following axioms: 

1. The field equations should follow from a variational principle. The independent 
variables in the action integral should be the components of the metric tensor. 

2. The action functional should be a scalar, 

3. The Euler-Lagrange equations that follow should be of second order in the space-time 
metric gab- 

Before introducing Hilbert 's form of the action for general relativity, a reminder is in 
order. We are looking for a condition in which all dynamical entities are in equilibrium with 
the background geometry. In terms of a variational principle this entails considering a single 
action for matter and geometry, which is invariant under the local group transformations 
imposed by general covariance, and whose equations of motion fully reflect the equivalence 
principle. To this end, we consider the following action 



s{n) = Sg{n) + s^n) = f{^g + ^m)dn. (2.1) 

Jn 



where the Lagrangian densities for geometry and matter, ^g and ^m respectively, are to 
be integrated over an open set il. of the spacetime M with compact closure. Now, we need 



to ensure that all the quantities involved in (2.1) are well defined. 



We consider first the part of the action describing the dynamics of the geometry. First 
we note that since the volume element of space-time dil is a scalar density of weight -|-1, 
it follows from axiom 2 that =Sfg should be a scalar densitjrjof weight —1. Also, from the 



^A tensor density differs from the corresponding tensor in that when transformed from one coordinate 
system to another, it gets multiplied by a certain power w of the Jacobian transformation, i.e. \/—g". The 
exponent of the Jacobian is called the weight of the density. A proper scalar should have density 0. The 
weight is additive under tensor product, that is: given two scalar densities of weight w and w respectively, 
their product is a scalar density of weight w ^ w (See Chapter 3 in [61] for a full discussion on general 
tensor densities). 
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third axiom above, ^g cannot be just ^/—g, it should also contain the derivatives of the 
metric. Thus, we could write the spacetime Lagrangian density as 

^g = V^Ag, (2.2) 

where Ag is a proper scalar constructed from covariant combinations of the metric and its 
derivatives only. 

The order of the Euler-Lagrange equations is twice the order of the highest derivative 
in the action integral, however it is not possible to construct a proper scalar out of com- 
binations of the metric and its first partial derivatives since the Christoffel symbols do not 
transform as tensors. Thus we need to include second order derivatives of the metric in A^. 
These terms will not contribute to the equations of motion if they can be collected into a 
divergence of an expression which vanishes on the boundary dQ, as it is shown below. The 
geometry part of the action is postulated to be 



Sg{n) = / ./^Rdn, (2.3) 

Jn 

where R is the curvature scalar. 

In a similar way, the Lagrangian density for the matter ^^ can be written as a product 
of \/—g with a proper scalar built out of covariant combinations of all the matter fields we 
observe. Notice that the metric itself may play a role in these combinations by producing 
the appropriate Lorentz scalars. In this sense, the metric enters the matter sector of the 
action in an algebraic way only. Thus, the matter part of the action is written as 

Smm= [ V^Arni^,^. a, g)dn (2.4) 

Jn 

where we denote collectively all the matter fields by the letter ^ and the semi-colon denotes 
metric compatible covariant differentiation. 

General covariance and spacetime coincidences 

The equations of motion are obtained in the usual way by requiring the action Smi^) to 
be stationary under variations of the fields in the interior of il. The following argument 
has been taken from Hawking and Ellis |33j . It is known as 'the hole argument' for reasons 
that will shortly become apparent, and it plays a central role while discussing the meaning 
of general covariance ( c.f. discussions in |54L [66] ) . 

A variation of the fields, including the metric, is realised by a one-parameter family of 
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diffeomorphisms from spacetime into itself, ifs G Diff(A4), which induces the maps 

f*s^ = '^{s,x) and ip*g = g{s,x), 



(2.5) 



where the parameter s takes its values in a sufficiently small open interval of the real 
numbers centred at the origin, and a; is a point in Ai, not necessarily inside il. such that 

1. ^(0,x) = ^(x), g(0,x) = g(x), i.e. ips=o = id^ 

2. ^(s, x) = ^(x), g(s, x) = g(x) when x is not in the interior of Q. 
Thus, the fields variations are denoted by 






S^ = — v?*^ 



s=0 



and (5g 



ds 



^s 



s=0 



(2.6) 



The term 'hole' now becomes evident. We are looking at variations which leave the 
value of the fields unchanged at each point of M , except inside an open region - the hole 
- Q. This implies, in particular, that the variations of the fields and the metric vanish on 
the boundary, dQ. 

General covariance is equivalent to the requirement of the invariance of the action under 
a general transformation of the spacetime manifold 



Smin;^,g) = Sm{n;^:^,^lg). 



(2.7) 



A crucial point to be noted is that a transformation ips pulls back both, the fields and 
the metric. This removes any physical significance from spacetime points. This fact is a 
statement in modern language of what Einstein referred as spacetime coincidences. 

2.1.3 Equations of motion 
Matter equations 

The requirement of invariance of the action under diffeomorphisms of spacetime imposed 



by general covariance, equation (2.7), makes evident what the matter dynamics is. From 



the definitions (2.6), the variation of the action is 



5S„ 



S^Sn^dift^ , 6^Smd^tg 



+ 



S^"^ ds (Jttg ds 



J s=0 



ds 



Smm 



s=0 



(2. 



This should hold for arbitrary variations of the fields and of the metric. Therefore, 



each term of the variation (2.8) should vanish independently. 
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The first term inside the brackets in (2.8) corresponds to variations of the action which 



keep the metric fixed, whilst the second term corresponds to variations in the metric which 
leave the values of the fields unchanged. We denote with a 'flat' superscript the first kind 
and with a 'sharp' the second. Thus, the flat variation of the action is 



6^S„ 



^^"'S^^i,+ ^S^^i,. 



a^ 



d^-a 



-gdn. 



(2.9) 



Since the parameter s is completely independent of spacetime, it follows from (2.5) that 

d 



(sna 



ds 



v:^ 



ds 



ifl^.a = S^,a, 



(2.10) 



we can integrate by parts the second term in (2.9) 



r dK 



""S'^^a^ 



-gdn 



dA 



m .b 



d^;a 



6'^ 



dK 



d^., 



S'-^ 



-ffdQ, 



(2.11) 



and write the 'flat' variation of the action for matter as 



5'Sr, 



dKn 



dArn 



6^^^/^dn + 






-gdn. 



(2.12) 



We can use Stokes' theorem in the last term to transform the total divergence into 
a surface integral over the boundary dQ where the field variations are identically zero, 
as previously argued. Therefore, the change in the action Sm inside fi under arbitrary 
variations of the fields ^, 






dAm 

1)¥ 



dAm 



-gdn = 0, 



(2.13) 



if they satisfy the Euler-Lagrange equations, i.e. if (2.13) is identically zero. These are the 
equations for the matter fields. 

The variation of Sm with respect to the metric is slightly more elaborate since in addi- 
tion to the field variation we also need to take into account the changes in the density \/—g 
and, possibly, in the covariant derivative of the matter fields induced by the transforma- 
tions ifs- However, since the connection is a coarser geometric object, its variations can 
be considered to be independent of the metric. The variation of the action with respect to 
the connection will not produce an equation of motion for it since the connection itself is 
not an observable field. Instead, it will fix the Levi-Civita connection of the metric g. We 
will follow Palatini' s approach and consider the connection variations as independent so 
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that 



s^s„ 



dAm 
dg 



r\ I 

'f7 5g 



does not include variations of the covariant derivatives of the matter fieldq£] . 



Let us start with the second term in (2.14). The change of the density 
respect to the metric is given by 






1 



dg 



1 



ab 



-gdg, 



ab 



--gg 



ab 



1 



;5 



ab 



-g, 



(2.14) 



-g with 



(2.15) 



so the 'sharp' variation (2.14) becomes 
where 



na6 



T'^'d'gabV^dn, 



T 



ab 



,9A„ 
^dg' 



ab 



Amg 



ab 



(2.16) 



(2.17) 



is a symmetric tensor which represents the energy and momentum distribution of the 
matter fields with respect to the undetermined background geometry. 

Finally, let us assume that the matter fields satisfy the Euler-Lagrange equations so 



that (2.13) vanishes. By the fundamental theorem of differential equations, the integral 



curves of ips correspond to the flow generated by some vector field ^, therefore, the variation 



of the metric, (2.6) can be written as 



(5g 



ds 



s=0 



^^gab — S.{a 



(2.18) 



where C^ denotes Lie-differentiation with respect to the vector field ^. Substituting this 



into (2.16) yields 



I I T-'^(a;b)V^dn, 



(2.19) 



and using the symmetry of the energy momentum tensor, we can integrate (2.19) by parts 
to obtain 



Jn 



-gdn 



lab 



n^aV^dQ 



IS^"'')* 



-gdn 



(2.20) 



where we can use Stokes' theorem again to transform the last term into a surface integral 
over the boundary of the hole - where the vector field ^ vanishes. Thus, the variation of 



^See chapter 2 in [33: for a discussion where these variations are included. 
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(2.16) becomes 



6^S„ 



'^ab 



T'^^CaV^dn. 



(2.21) 



n 



The equations of motion follow from the vanishing of (2.21) for the arbitrary vector field 



^ inducing the metric variation (2.18). These are equivalent to the local conservation law 

rpab 



0. 



(2.22) 



Before we continue with the derivation of the equations of motion for general relativ- 
ity, let us summarize the consequences of Hilbert's variational approach for matter. The 
requirement of general covariance in the formulation of physical laws is equivalent to the 



invariance of the action Sm under diffeomorphisms of spacetime, equation (2.7). This sym 



metry implies the dynamics of matter, given by (2.13) and (2.22) as Noether identities 



and, therefore, they do not need to be postulated independently. However, this is nothing 
particular about general relativity, as discussed in the previous section. So far we only 
have re- written classical dynamics in a generally covariant manner. The next step is to 
incorporate the principle of equivalence to identify the full dynamics of spacetime and 
matter. 

Geometry equations 

The last stage of the variation formalism is to obtain the dynamics of the geometry implied 
by general covariance. The procedure is straightforward. As previously stated, we follow 
Palatini's formalism by considering the connection as an independent field. We stress again 
that no equations of motion for the connection would arise from its variations, those will 
only fix the Levi-Civita connectiorj^ 

Following the previous convention, we denote with a sharp super-script a variation of 
the metric which leaves the other fields (the connection) fixed, thus, the geometry equations 
of motion are those for which 



d^ 
di 



Sgm 



<5»5„ 



s=0 



-g5^R + R5^^^ 



dO = 0. 



(2.23) 



We have already calculated the variation of the density ^/—g in (2.15). Hence, we only 



need to compute the variation of the curvature scalar. We do this in two steps. First, we 
consider the variation 

5^R = Rab5^g'''' + g''H^Rab. (2.24) 

We notice that the Ricci tensor depends only on the connection coefficients, and therefore 



''The proof of this can be found in Section 5.3.1 of [41] 
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its 'sharp' variations are identically zero. Thus we only need to compute the covariant 
components of the metric variations 



"ff 



ab 



-9'"'g'''s^gcd. 



(2.25) 



Substituting the derivative of y/—g given by (2.15), the variations of the curvature 
scalar (2.24) and the covariant components of the metric (2.25) into the variation of Sg, 



equation (2.23) becomes 



S^S„ 



~iab 



G^'6^gabV^dn, 



where 



f^ab P'*^ 



Rg' 



ab 



(2.26) 



(2.27) 



is called the Einstein tensor. Following arguments analogous to those leading to equation 



(2.22) for the stress-energy tensor, (2.27) satisfies the Bianchi identities^ 



G' 



ab 



0. 



(2.28) 



This completes the variational program. It just remains to read the dynamics encoded 
in the symmetry of general relativity. 

Einstein equations and the equivalence principle 

The complete dynamics of general relativity follows solely from the invariance of the 



Einstein- Hilbert action (2.1) under the group of diffeomorphisms of spacetime. Assum- 



ing that the matter fields satisfy the Euler-Lagrange equations, general covariance implies 



d^ 
ds 



s{n) 



s=0 



G 



ab 



_rpab 

2 



''gab^ 



-gdQ. = 0, 



(2.29) 



therefore, from the arbitrariness of the diffeomorphisms ips-, the equations of motion for 
general relativity are 

(2.30) 



Rab — 7,gabR — SvrTaf,. 



These are ten, coupled, non-linear, second order partial differential equations for the metric 
coefficients as required by Hilbert's axioms and obtained as a direct consequence of the 
general covariance of physical laws. Encoded in them is the principle of equivalence. Matter 
and geometry appear on an equal footing. Their interaction does not distinguish any 
specific set of physical fields encoded in the stress-energy tensor provided they have the 

*This was the fact that originaUy motivated Einstein to postulate the tensor G for the right hand side 
in equations (|2.30||. 
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same energy-niomentum distribution [20j. This is the most dramatic insight of general 
relativity; spacetime curvature is gravity and energy is its source. 

The class of solutions admitted by the system of Einstein's equations is very large. 
In principle, any metric which is differentiable up to second order in a given region of 
spacetime could be used to define an energy-momentum tensor. This would leave the task 
of finding which kind of 'matter dynamics' is defined by the geometry prescribed by the 
chosen metric. The matter models of physical interest, however, are not that arbitrary. 
There is a set of conditions which must be satisfied by the energy-momentum tensor in 
order to consider a theory as physical. In practice, one proposes matter models which 
satisfy a physically relevant condition and proceeds to find the geometry associated with 
them. 

2.2 Irreversible processes 

It may seem that all there is left for classical theories of physics is to write down a Lag- 
rangian density or an energy-momentum tensor which satisfies some energy condition and 
then solve Einstein's equations for the particular matter model of interest. Here is a sharp 
distinction between the interpretation of a fundamental theory of nature and the observable 
phenomena it was abstracted from. 

General relativity is the product of the mathematical language of general covariance 
and the physical principle of equivalence whose dynamics is obtained from the Einstein- 
Hilbert action. As a dynamical theory it is complete. Then, the question is if dynamics is 
a faithful description of all natural motions, whether every process described by dynamics 
can be found in nature. 

We understand a natural process as a consecutive set of motions changing the properties 
describing a physical system from one configuration to another. Implicit in this reasoning 
is the notion of an evolution which makes a clear distinction between the directions of 
the past and the future, i.e. it presupposes an arrow of time. In the vast majority of 
mechanical contexts, such a distinction is merely conventional. All the dynamical laws are 
symmetric with respect to a parameter which we identify with time. Thus, we normally 
follow the convention that a system naturally evolves in the positive direction of this time 
parameter. If the goal of any dynamical theory is to give an account of the natural evolution 
of physical situations we encounter in the universe, we are faced with the empirical fact that 
most processes we experience do not possess such a time symmetry - they are irreversible. 

The aim of this section is to highlight the questions raised by incorporating the Laws 



of Thermodynamics (c./. Section 1.1) to the empirical basis of a dynamical theory. 
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2.2.1 Equilibrium thermodynamics 

The detailed dynamics of a system composed of a large number of interacting particles 
in an arbitrary manner is prohibitively complicated. Nevertheless we can still extract a 
significant amount of physically relevant information if certain restrictions on the collective 
motion of the individual components are imposed. In this manner, we can get some insight 
into the most complex scenarios by considering progressively more elaborate steps. 

We have argued that thermal equilibrium corresponds to a highly idealised state in 
which the dynamics of a system are time symmetrical. In this section we will clarify this 
notion from a macroscopic point of view, leaving the micro-physics details for the following 
chapter. Therefore, for any macroscopic system, we say: 

Thermal equilibrium is that state in which the total energy is conserved and 
its interna/Qenevgy is maximally distributed. 

Notice that this definition does not make any assumption about the material nature of 
the system. As far as thermal equilibrium is concerned, it makes no difference if we are 
dealing with atoms, molecules, photons or any other exotic type of matter. This is what 
gives thermodynamics its universal character. 

The assumption of thermal equilibrium allows us to reduce the very large number of 
degrees of freedom of a system to a handful of the relevant variables that describe its state. 
In general, not all such variables are independent, being connected through an equation 
of state. Such equations are not part, nor can be deduced from, the thermodynamic 
hypotheses. They contain the material information of the system we try to analyse. We 
shall see that not every relation between the state variables of a given thermodynamic 
system corresponds to a real matter model, for they can lead to dynamical inconsistencies 
within the underlying fundamental theory. We will explore these matters in further detail 
in Chapter 6. 

As pointed out in [6], the number of independent variables can also be reduced if 
the system has an overall additivity property. To illustrate this, let us consider the case 
of a fluid consisting of a single species of particles whose total energy of E is related 
to its entropy S, volume V and number of particles N through and equation of state 
E = E{S, N, V). In this case, we can express the change in the total energy of the system 
as 

dE = TdS + fidN - pdV (2.31) 



^Here, we understand total energy as a global quantity, while for internal energy we mean that possessed 
by the system's sub-parts. 
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where, by definition, the temperature T, chemical potential fi and pressure p are 

dE dE , dE 






dV 



(2.32) 



In this example, the variables E, S, N and V satisfy the following scaling property. Let 
us denote with a tilde the scaled variables 



S = XS, N = XN and V = XV, where A G M. 



Then, the total energy of the system becomes 



E{S,N,V) = XE{S,N,V). 



(2.33) 



(2.34) 



whilst the variables (2.32) remain unchanged. The local form (2.31) of the scaled energy 
takes the form 

dE = XdE + EdX = TdS + ^xdN-pdV. (2.35) 



Comparing the last two equalities it easy to check that the total energy is 



E = TS + nN- pV. 



(2.36) 



This result is known as the fundamental Euler relation of equilibrium thermodynamics. If 



instead of working with absolute quantities we use densities, we can rewrite (2.36) in the 
form 

p + p = Ts + fin, (2.37) 

where p = E/V, s = S/V and n = N/V are the energy, entropy and particle number 



densities, respectively. Now, if we let the scale factor to be A = 1/V, the energy (2.34) 
becomes 

E = (5, N, V) = E{s, n, 1) = p{n, s), (2.38) 



and the differential form (2.31) reduces to 



dp = Tds + p,dn. 



(2.39) 



Therefore, whenever a system admits an additivity property of the form (2.34), we can 



effectively reduce the number of independent variables by one. In our case, equation 



(2.39) will henceforth be referred as the equilibrium Gibbs' relation. It will serve as our 



thermodynamic reference point in the relativistic discussions of the forthcoming chapters. 
Finally, let us note that thermal equilibrium implies that the rate of change of the total 
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energy with respect to time vanishes identically. Thus, in the case of a single species fluid, 
the entropy and particle productions should be in balance. In particular, if the particle 
number is conserved, the entropy, as required by the definition, should be maximal. We 
will give a clearer explanation of the meaning of maximally distributed energy in the next 
chapter, when we introduce the microscopic motivation for the concept of entropy. For the 



time being, let us note that the Gibbs' relation (2.39) completely characterises the state 



of thermal equilibrium for a fluic Therefore, for a fluid with constant particle number 



density, the above definition of thermal equilibrium is equivalent to the weaker version: 

Thermal Equilibrium (isentropic) is that state for which the entropy dens- 
ity of the system is a conserved quantity, i.e. 

ds = 0. (2.40) 

Notice that this definition of thermal equilibrium is indeed a time symmetric one since 
the requirement of maximal energy distribution implies that there are no hotter or colder 
regions; there is no net energy transport from one part of the fluid to the other. In this 
sense, equilibrium thermodynamics gives an account of collective phenomena in a time 
reversible fashion. 

The above isentropic definition of thermal equilibrium differs from the general one in a 
crucial point. Suppose we have a thermodynamic system in thermal equilibrium. There are 
essentially two different ways we can perturb the system away from equilibrium: reversibly 
and irreversibly. The difference between the two is that, in the former, the total energy of 
the system is modified without changing its entropy while, in the latter, entropy is not a 
conserved quantity, even when the energy itself may be conserved. 

Here, we have been implicitly assuming that time indeed 'flows' in one direction, i.e. 
towards the future. Thus, when we speak about reversibility we mean that, from the 
present state of a system, we have enough information to retrodict the system's previous 
state. This is obviously the case for thermal equilibrium, although it may come as a 
surprise that it is also the case of some systems away from it. In such a case, the only 
relevant auxiliary hypothesis are the zero-th and first laws of thermodynamics discussed 
in the introduction (c./. Section |1.1[ ); the 'time-symmetric' laws. It is the second law 
that describes systems which, internalljj^ can distinguish past from future; we gradually 
lose our ability to retrodict the past while the final state becomes progressively more 
predictable. Moreover, if the final state of thermal equilibrium is reached, we have lost all 



^°The extension to the case with more than one species is straightforward. 

^^At the level of description, the distinction between past and future is purely a convention on the 
direction of the future, not a true feature of the dynamics. 
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previous information about the systems evolution; from the system's point of view, it has 



never departed from it We can give a more precise statement of the second law for a 



system which has been left to evolve without any external influence: 

Second law of thermodynamics The entropy of a system, which internally 
can distinguish the orientation of time, always increases towards the future, i.e. 

f>0. (2.41) 



Let us denote the rate of change of entropy (2.41) by Fg and call this quantity the 
entropy production. To give account of irreversible processes we need a detailed study of 
situations when the entropy production does not vanish. We cannot be satisfied with the 
simple inequality stated by the second law. It is our task to relate entropy production to 
a well-defined physical processes. 

2.2.2 Linear Irreversible Thermodynamics 

This subsection is devoted to the first step in the study of irreversible processes. The 
physics of systems far away from thermal equilibrium remains a topical problem. Thus, 
the natural approach towards a theory of irreversible processes would seem to require an 
expansion in deviations from the state of equilibrium. Here, we give an account of the 
Linear Irreversible Thermodynamics (LIT) programme pioneered by Onsager and later 
developed by Prigogine |64^ I23j . The central result of this subsection are the Onsager 
reciprocity relations, of which Fourier's law is a consequence. We will argue that, although 
it has been experimentally verified to a good level of accuracy, LIT will necessarily fail 
when the relaxation time scale of thermal phenomena is long compared with the dynamical 
one - e.g. in superfluid helium - even when the deviations from equilibrium are not large. 
Therefore, although LIT will not play a significant role in this thesis, this brief introduction 
will serve as a motivation for the Extended Irreversible Thermodynamics (FIT) programme 
which permeates the multi- fluid approach to heat conduction presented in Chapter 5. 

A simple evaluation of the entropy production becomes possible if we assume that en- 
tropy away from equilibrium depends on the same variables as it does while in equilibrium. 
We will introduce no restriction as to the nature of the system except that we shall assume 
that, although the system as a whole will not be in equilibrium, there exists at every point 



^^Determinism is lost! One could argue that a system takes an infinitely long time to completely relax 
to equilibrium and, therefore, its evolution is deterministic for all times. Surprisingly, we can find purely 
mechanical examples for which determinism is also lost (c./. Norton's slippery dome [55] )■ The solution to 
this apparent loss of predictability can be traced down to the differentiability of the system's constraints 
(in the case of Norton's dome the surface is analytic everywhere except at the summit). 
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a state of local equilibrium for which the entropy is given by the classical Gibbs formula 



(2.39). This restricts us to systems whose evolution may be described in terms of mac- 
roscopic thermodynamics and hydrodynamics without explicit introduction of molecular 
concepts. 



The central quantity is the entropy production per unit time (2.41). We should notice 
that this will become an issue when dealing with relativistic systems, when we will have 
to refer to a particular frame to have a definite interpretation. For the time being we shall 
consider the Newtonian framework where thermodynamic systems are Galilean invariant, 
and hence have a universal notion of time. In this case, the entropy production per unit 
time may be expressed as the volume integral 

— = J adV>0, (2.42) 

where a denotes the entropy production per unit time and volume. 

It is very rare that thermodynamic processes occur in complete isolation. It was noted 
by Onsager that when two or more irreversible processes take place simultaneously in a 
thermodynamic system, they may interfere with one another [5^: leading to the discovery 
of the celebrated reciprocity relations. The central assumption in Onsager's reasoning was 
that of microscopic reversibility, which guarantees an independent balancing of different 
classes of molecular processes maintaining statistical equilibrium. This balance is a crucial 
point if we want to extend the LIT programme to physical regimes where non-inertial 
effects may not preserve the statistical equilibrium of the system. 

The LIT programme posits that the entropy production o" is a bilinear form of the 
generalised forces Xi (affinities, gradients of temperature, gradients of chemical potentials, 
etc.), and fluxes Jj representing the rates of the various irreversible processes taking place 
(chemical reactions, heat flow, diffusion, etc.) 



Y^ fX^ > 0. (2.43) 



It should be noted that the validity of the Gibbs relation (2.31) is assumed in (2.42). 
The Gibbs relation, however, can only be established in some neighbourhood of equilibrium 
which defines the region of local equilibrium. At thermal equilibrium the fluxes and forces 
vanish identically, i.e. 

/ = 0, X, = (2.44) 

for all irreversible processes simultaneously. This is the justifying fact that allows us to 
assume linear homogeneous relations between the forces and the fluxes for conditions near 
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equilibrium. 

In this representation, the fluxes are the unknown quantities, and the forces are given 
functions of the state variables or their gradients. Thus, we can then write 

j' = r{Xi,X2,..,Xn). (2.45) 

Therefore, in the neighbourhood of equilibrium, we can make the expansion 

•^^ = E(|^)^^ + ^(^^-^'^)' (2.46) 



and neglect second order terms ^^ to obtain the fluxes as linear transformations of the forces 

r = "^V^Xj. (2.47) 

j 

Here the coefficients L*-' are measurable constants which are defined in terms of correlation 
functions 1491. 



Substituting the linearised fluxes (2.47) into the entropy production (2.43) gives us the 
form 

a = ^V^XiXj. (2.48) 

«J 

The positivity of entropy production constrains the coefficients Ljj to be a quadratic form. 
Thus, the simplest constitutive equations which guarantee the semi-positiveness of the rate 
of entropy production are given by Onsager's reciprocity relations 

L'^ = V\ (2.49) 

The Onsager relations are the central result in the linear thermodynamics of irreversible 
processes programme. Their strength resides in their generality; they showed for the first 
time that non-equilibrium thermodynamics leads to general results independent of the 
molecular model used [M] . 

Onsager's reasoning was inspired by William Thomson's idea of a variational principle 
for dissipative systems, namely, a principle of least dissipation of energy. There is no 
evidence of such a principle for a generic non-equilibrium system. However, for relatively 
small deviations from thermal equilibrium one can argue in favour of a principle of least 
entropy production |23j . To illustrate this point, let us consider the case of heat conduction 



^Microscopic reversibility eissumption. 
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alone. In this case we are left only with the single pair 



dT 

Ji = -K^, (2.50a) 

ox 

d /I 



where J'^ is the one-dimensional heat flow, T is the absolute temperature and k is the 



thermal resistivity. The volume integral (2.42) becomes 

dS f K fdT^ ^ 



_ \ \ ^V. (2.51) 

dt J T^ \dx ' ^ ' 

If in the whole system the temperature deviates locally from some average Tq 

T = Tq + 6T, (2.52) 

neglecting higher order terms, the perturbed entropy production is 

ds f K (ostW^^ 

K fd5TV 



Since the temperature is a function of the position only, we write the temperature perturb- 
ation in terms of the displacement 5x as 

dT 
6T = -—6x. (2.54) 



Hence, since 5T /Tq <^ 1, the variation of entropy production (2.53) becomes 



''it-khi^^'^'''- <"=) 



Note that the minimum, of (2.55) corresponds to the stationary state away from equi- 
librium 

d'^T 
.^ = 0. (2.56) 

Prigogine referred to this result as an example of a kind of inertial property of non- 
equilibrium systems 

. . . When given boundary conditions prevent the system from reaching thermo- 
dynamic equilibrium (that is, zero entropy production) the system settles down 
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to the state of least dissipation . . 

It should be stressed that the theorem of minimum entropy production is formulated 
under the assumption of small deviations from equilibrium. It can easily be shown that 
systems far from thermodynamic equilibrium could behave in quite a different manner, 
in fact, directly opposite from the one expected from the minimum entropy production 
theorem. 

It was also noticed in [38] that this approach has some other limitations, the most 
important for us being that of the parabolic nature of Fourier's law. This implies that 
thermal disturbances propagate information with boundless speed. Cattaneo proposed 
a solution to this problem by taking into account the microscopic time-scale in which 
such disturbances would propagate information in a material medium, thus deducing a 
hyperbolic equation for heat. 



Hyperbolic heat conduction: a road to Extended Thermodynamics 



The one dimensional version of Fourier's law (2.50a) together with an energy balance of 
the form 



dt 



T{t,x) 



d_ 
dx 



J'{t,x), 



leads to the heat equation 



d_ 
di 



T{t,x) 



K^T(t,x) 



0. 



(2.57) 



(2.58) 



We may write the solution of (2.58) for an initial value problem in an unbounded region 
of space in the form 



r(t,x) 



1 



(47rKt)3/2 



r(0, y) exp 



(y 



AKt 



dy. 



(2.59) 



This solution implies that T{t, x) is non-zero zero for all x and t > 0, even when r(0, x) 
may have support on a finite interval only. 

The above result is far from satisfactory. One would expect that temperature disturb- 
ances have a bound on their propagation speeds given by the material internal structure. 
In this interpretation, heat conduction can be understood as the energy exchange during 
collisions of the molecules forming the material. The chance of being hit depends on the 
length a molecule travels before hitting another one. This length is the mean free path. 
Cattaneo proposed a modification to Fourier's law to incorporate this view 



Ji 



dT 

dx 



d_ fdT\ 



(2.60) 
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where r corresponds to a real parameter with units of time. He called (2.60) 'non- 
stationary' Fourier's law. This leads to 



di 



T{t,x) 



52 d f d^ 



(2.611 



Still, this is not a hyperbolic differential equation, this attempt did not remedied the 



problem. Thus, Cattaneo proposed a further modification of (2.60). He assumed that 



the the operator T{d/dt) was 'sufficiently' small so the difference between J'^ and J'^ is 
negligible. Thus, he could approximate the non-stationary Fourier's law by 



J«(t,x)+f^J''(t,x) 



-K-r(t,x), 



(2.62) 



where f = t / k. This result is called the Cattaneo equation. This will be the reference 
point in the coming chapters, when we deal with hyperbolic theories of dissipation. The 
reason for this is that it leads to a telegraph equation for thermal disturbances, i.e. 



^£^(^'^)+l^(*'^) 






0. 



(2.63) 



Equation (2.63) is hyperbolic if the parameter f > 0. In such a case, it predicts the 



propagation of heat pulses at the finite speed 



(2.64) 



Thus, we see that for arbitrary small values of the 'relaxation' time f the characteristic 



speeds become unbounded. Figure |2.1| shows a comparison of the evolution of a pulse by 
considering the heat and the telegraph equations, respectively. 



Note that in order to obtain the hyperbolic heat equation (2.63) two assumptions were 



made. The first is that the 'modified' heat flux J'^ is no longer related to a 'force gradient' 
as in the ordinary Fourier's law. The second is the assumption that the quantity 



d fdT 



dt V dx 



(2.65) 



is small, thus setting a dissipation time-scale which is long compared to the system's 
dynamical evolution. This kind of behaviour is better seen in liquids, as in the case of the 
second sound in superfluid helium or the dissipation to solidification time scales of molten 
polymeric fluids ^^ In these cases, heat acts as a 'slow' variable of the process. Thus, one 



^'*An example of a polymeric fluid is 'gel-wax', a fluid typically used in candle making. One can try the 
following experiment: put some molten gel wax in a Petri dish to form a layer of about 1mm. Heat from 
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could treat the heat flux as if it were a state variable. This constitutes the motivation to 
introduce Extended Irreversible Thermodynamics. 

2.2.3 Extended Irreversible Thermodynamics 

Prom the discussion in the previous section, the reader may have some feeling of discomfort 
about the applicability of the heat equation. For many 'practical' applications, Fourier's 
law is a very good approximation. However, when one is interested in sufficiently rapid 
phenomena, the LIT description is not good enough. Moreover, when the underlying 



principles of the dynamical theorj ^^ prevent signalling faster than a critical speed, then 
the traditional approach is simply unacceptable. Therefore, there is a clear conflict between 
Linear Irreversible Thermodynamics and relativity. 

The general philosophy behind a thermodynamic description of macroscopic phenomena 
is that of averaging. In such a process, one reduces the (very large) number of degrees 
of freedom, to a manageable set of 'thermodynamic variables' which are directly related 
to conserved quantities, i.e. mass, energy and momentum. All other variables decay 
very rapidly in the averaging process and thus they are not conserved, e.g. heat. If, 
however, the time scale of certain phenomena is 'fast enough', the decay of some of these 
variables may become slow compared to the total duration of the phenomenon. In this case, 
the otherwise not conserved fast variable, can be treated as a slow, conserved one. The 
Extended Irreversible Thermodynamics programme (FIT) was designed to systematically 
approach thermodynamic phenomena whose set of conserved slow variables is enlarged by 
a set of 'dissipative' yet slow ones which complement the thermodynamic description of a 
system evolving away from thermal equilibrium. The interested reader is referred to the 
work by Jou [38j . 

Continuing our previous fluid example (c./. Section 2.2.1), in FIT one considers the 
energy of the system to be a function of the equilibrium (conserved) quantities and the 
dissipative (slow) fluxes. In our case, we have 

E = E{N,S,V,J''), (2.66) 

or, equivalently, the energy density is written as 

p = p{n,s,q). (2.67) 



below in a uniform manner trying to preserve a steep temperature gradient between the bottom and the 
surface of the fluid. With enough luck hexagonal convection patterns (Benard cells) will form. Remove the 
Petri dish from the heat source. If the cooling down is sufficiently uniform, the wax will stiffen before the 
patterns are completely dissipated. 

^^The underlying Lorentzian geometry encoding the 'principle of causality' through general covariance. 
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Figure 2.1: This figure shows a comparison of the evolution of a thermal disturbance 
modelled with the heat (left) and telegrapher (right) equations. The vertical axis represents 
time, while the horizontal corresponds to one-dimensional space. The initial condition is 
a delta function at the origin of the spatial direction. The evolution occurs towards the 
top of the page. We can see that, for the heat equation, the pulse spreads over the entire 
region without any lag. On the contrary, the evolution provided by the telegraph equation 
propagates a thermal signal which dissipates away to merge the long term solution of the 
heat equation (the transition between 'fast' and 'slow' heat flux). The bottom graph shows 
a how the two distinct evolutions become indistinguishable in the 'long' run. Here, the 
'z'-axis corresponds to temperature, while the spatial and temporal directions should be 
obvious from the description of the top images. 
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Thus, we can write an 'extended' Gibbs relation of the form 

dp = fidn + Tds + fiqdq, (2.68) 

where fig stands for the 'heat chemical potential' given by 

N^^- (2-69) 

It is one of the aims of this work to show how an extended Gibbs relation leads naturally 
to a hyperbolic description of heat conduction in a relativistic setting. The details of this 
will be given in Chapter 5. 

2.3 Irreversible processes in general relativity 

In this chapter, we have introduced the two main components of this thesis. On the one 
hand, we consider the general theory of relativity as our fundamental dynamical theory 
(c./. Section 2.1, footnote 3). The purpose of presenting the full variational construction 
from the Einstein-Hilbert action was two-fold: 

1. it allowed us to spell out unambiguously the meaning of general covariance in the 
context of this work, and 

2. to highlight that the connection of the mathematical formulation of a dynamical 
theory with the empirical world is manifest through the Lagrangian function and its 
symmetries. 

It is hard to overestimate the importance if these two points, specially to correctly interpret 
the outcome of chapters 5 and 6 - the core of this manuscript. 

On the other hand, we acknowledge the fact that many dynamical processes allowed by 
general relativity, have never been observed in nature ^^ and therefore, it may seem the they 
are not allowed by some physical law as a consequence of its 'time-symmetric' naturq^ 
As discussed in chapter 1, here we adopt the view that thermodynamics is a collection of 
self-consistent additional hypothesis - the laws of thermodynamics - that we must add to 



our dynamical theory to obtain a stronger description about nature ^® 



^^This kind of induction does not imply that such processes wiU never be observed, but we may have to 
wait a 'very' long time. 



The choice of the direction of time is a matter of convention. 



17 

^^The addition of the laws of thermodynamics to any dynamical theory 'forbids' a larger class of processes 
to occur than the dynamics alone would. In this sense it is 'easier', at least in principle, to falsify. Thus, 
we obtain a stronger description of nature [62] . 
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If we agree that the principles of general relativity and the laws of thermodynamics 
form the empirical basis of our macroscopic description of nature, how do we assess our 
individual models of physical processes? In the forthcoming discussion we will show that 
in essence, the physical information contained in the Lagrangian function is equivalent to 
the one given by an equation of state. Thus, in a more complete view of natural processes, 
we can describe a system's thermodynamics in a dynamical way or, equivalently, we can 
give a thermodynamic description of its dynamics. 

With this, we close our discussion about the framework of this thesis. 
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Irreversible processes in general relativity represent an outstanding challenge for our ef- 
forts to understand the dynamics of macroscopic systems. The equations governing such 
dynamics were laid down in the previous chapter. In this chapter we deal with the matter 
content of those equations. The aim is to give a sound microscopic motivation for our later 
developments of relativistic non-equilibrium thermodynamics. 

Although most of the material in chapter is not new, we have tried to carefully spell out 
subtleties often omitted in more specialised literature. In particular, on arguments typically 
invoked in 'proofs' of the fundamental theorem of irreversible processes; Boltzmann's %- 
theorem. The following is an account of the results of non-equilibrium relativistic kinetic 
theory which are relevant to keep this work self-contained. Central discussions relating 
to the Einstein-Boltzmann or Einstein- Vlasov systems, have been entirely omitted. The 
interested reader is referred to the classic work of Synge [71], Tauber & Weinberg [73] . 
Stewart [69], Ehlers [26], Treciokas & Ehis f73], Elhs, Matr avers & Treciokas f28] and 
references therein for a range of different aspects of relativistic kinetic theory. 

The spirit of the following discussion, as a preamble to our treatment of relativistic 
dissipative fluids, is the one highlighted by Treciokas & Ellis [3]- Therefore, we will use 
kinetic theory to assist us to obtain the form of the macroscopic transport equations. 
Hence, rather than use it to examine the nature of the fluid theory, it will serve as a 
motivation for it. 

3.1 The symplectic structure of the relativistic phase-space 

This section is the starting point for our exploration of processes in nature which exhibit 
a 'time-asymmetric' evolution. It was stated in the previous chapter (see footnote 3) 
that, within this work, we consider a theory as physical only if it can accommodate the 
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features of our observed world, moreover, we say that a theory is fundamental if it can 
be derived from a variational principle. One can immediately raise an objection - with 
well justified claims - that, under such consideration, we would have to exclude dissipative 
systems altogether. The reason for this is that, in the strict sense, variational techniques 
only apply to conservative systems. The view sustained in this thesis is more relaxed. 
Therefore, the rest of the present work will be devoted to argue that, through the insights 
stemming from the tenets of relativity, a variational construction for transport phenomena 
is possible. 

The central aim of this section is to present a covariant description of the underlying 
microscopic theory of relativistic hydrodynamics. The reason for this exercise is two-fold. 
On the one hand, it will exhibit the role played by general covariance in the context 
of classical mechanics when formulated in terms of states on the relativistic phase-space. 
Furthermore, it will remove the notion of 'universal time' from the formulation of the 
second law of thermodynamics. 

The work presented in this section was motivated by the relativistic treatment of kinetic 
theory as given by Synge [7T], Stewart |69j and Ehlers [26]. The mathematical account is 
well known for non-relativistic systems and can be found in the classical text by Arnold 
[8]. A relativistic version of the same symplectic formalism has been given by Rovelli 
|66j . Hence, this brief introduction to relativistic kinetic theory is not new. However, 
the discussion regarding the covariant formulation of the second law of thermodynamics is 
presented from a novel point of view, allowing us to introduce the notion of 'entropy as a 
fluid'. 

Let us commence our discussion by considering a simple gas, i.e. a collection of a large 
number A^ of particles of a single species. As the particles move through spacetime collisions 
might occur. We will assume for the moment that they are completely elastic so that 
momentum, mass and number of particles remain conserved. When the particles are not 
colliding, the principle of equivalence asserts that they move along the free fall trajectories 
of the spacetime metric g, i.e. the non-intersecting curves, 7 : / C M — > M, which satisfy 
the equations of motion obtained from the Lagrangian function, A : TAi — > M, 

A(x^x") = ^5afe(x^)x'^i^ (3.1) 

Here x° denotes the components of an element of the tangent space T^Ad at a point x 
on the spacetime manifold M. This defines the set of 'natural' coordinates - the eight 
functions (x",x") - for the tangent bundle of spacetime TM. Through the rest of this 
work, a dot represents differentiation with respect to the afHnqJ parameter A. 



^ A parameter is called affine if it is related to proper time by an affine transformation: r — )■ A = ar + 6 
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The spacetime metric provides us with an isomorphism between the tangent and the 
co-tangent bundles of spacetime. Therefore, we can alternatively describe the same free- 
fall trajectories on T*Ai. Let us define the canonical momentum conjugate to x" by the 
spacetime one- form p £ T*A4 whose components are given by 

dL , , 

In this case, the free-fall dynamics of a system of test particles is encoded in the free 
Hamiltonian function, H : T*A4 — > M, defined by 

H{x'',Pa) = ^g'''{x')PaPb. (3.3) 

Here, the set of natural coordinates for T*M is given by the eight coordinate functions 
{x"',Pi,). As the reader may anticipate, the dynamics is independent of such a particular 
choice of coordinates. However, from an operational point of view, the natural coordin- 
ates are indeed 'preferred'. In the context of classical mechanics they are referred to as 
'canonical coordinates'. 

Notice the use of the covariant and contravariant metric representations in the Lag- 
rangian and Hamiltonian descriptions respectively. Their appearance highlights the role 
of the metric components as potentials of the gravitational field in a non-relativistic situ- 
atiorR 



Since the Hamiltonian (3.3) does not depend explicitly on the parametrisation of the 



where a and b are constants and r denotes proper time. 

■^Consider the weak field limit of general relativity. In this case we can decompose the spacetime metric 
into a Minkowski part and a perturbation 

gab ^Vab + hab, where r?a6 = diag(-l, 1, 1, 1) and l^atl < 1- 

It follows from the definition of the inverse metric 

Sa = gacg"'' = {rjac + hac) (rj"^ + h"'' ] = Sj" + Tjach"'' + hacf)"^ + 0{h?). 

Thus, keeping only the linear terms in h, it follows that 

-qaJi" + hac-n" = 0, 
and we have that the inverse metric is simply given by 

ah ab j ab 

g =11 -h . 

To obtain the Newtonian limit, we use the equivalence principle, i.e. from the geodesic equation and the 
equivalence between 'inertial' and 'gravitational' masses of test particles, to make the identification 

ftoo = -2$, 
where $ is the Newtonian gravitational potential. Thus, we can compute both the Lagrangian and Hamilto- 



nian [(3.11 and (3.31 respectively], for a test particle of mass m moving with some velocity u to obtain the 
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curve, its numerical value is constant on the hypersurfaces defined by the pseudo-sphere 

9''\x)paPb = -rr?. (3.4) 

If ?Ti > 0, we call m the mass of a test particle. Two curves passing through the same 
spacetime event with the same direction differ only in their parametrisation. Thus we 
see that their momenta define different "shells". We say that the particles have different 



masses and consequently we call the pseudo-sphere (3.4) the mass shell for particles with 
mass m. 

In non-relativistic classical mechanics, the phase-space of a system of N particles is 
defined by the 6A^-dimensional manifold whose elements are the collection of spatial loca- 
tions and momenta of each particle, i.e. the state of the system. The natural generalisation 
to a relativistic case would be the 8A^-dimensional manifold with the restriction imposed by 
the underlying Lorentzian symmetry of spacetime. The relativistic single-particle phase- 
space is then the collection of points x°' in spacetime together with all the allowed momenta 
Pa through each point. This implies that pa should either be timelike or null and future 
directed. Thus we say that the relativistic single-particle phase-space V is the future non- 
spacelike part of T*M. whose boundary dV corresponds to states of massless particles. 
Alternatively, we can define the relativistic phase space in terms of the mass of the test 



particles. The mass shells (3.4) form 7-dimensional subsets of T*M.. Let us denote by 
Vm each subset for different values of m. We define the single particle phase-space as the 
collection of all the individual mass-shells 



V = [jVm. (3.5) 



The dynamics of a system described by a Hamiltonian function follows from an action 
principle. However, the role played by the Hamiltonian is different from the one played 
by the Lagrangian in the variational principle discussed in Chapter 2. This is a direct 
consequence of the symplectic structure of V. 

3.1.1 Forms and motion 

It was previously said that the 'natural' choice of coordinates for T*M was the 'pair' 
{x"',Pa), where the x"" are the coordinates representing a point x €z Ai and the pa are the 



familiar expressions of classical mechanics (modulo an infinite additive constant!) 

A — -mv — m^ and H — —mv + m^. 
2 2 
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components of the one-form 

a = padx°', where aGT*M.. (3-6) 

By definition, the one-forms that we have encountered so far are elements of the dual to 
the tangent space at each point of spacetime, i.e. they are functions which assign a real 
number to each tangent vector of spacetime. In particular 

a : T^M — > M. (3.7) 

The natural projection of the co-tangent bundle, vr : T*M. — > M., induces the push- 
forward map vr* : TV — > TAi. Therefore, given a tangent vector X £ TmV, we observe 



that the action of the one-form (3.7) on the pushed-forward vector vr^X G TxA4, whose 



base point is a; = '/r(m), defines the phase-space one-form 

Q : TmV — > M. (3.8) 

We can express this by 

a [X] (m) = a [ir^X] (x), where a G T^V. (3.9) 



Equation (3.9) is the definition of the canonical one-form or symplectic potential. Written 
in this form it is clear that a is a one-form on phase-space, not on spacetime. The situation 
is summarised by the commutative diagram; 

T„^V > T^M 

(3.10) 



This rather abstract treatment pays immediate rewards. It is not difficult to show 
that the symplectic structure of the co-tangent bundle (in particular of the phase-space) 
is simply given by the closed two-form 

oj = da = dpa A dx"" . (3.11) 




Here the coordinate expression follows directly from (3.7). Notice that the warning sign is 



in place since the coordinate representation is defined from dd and, therefore, it belongs to 
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a different space! It is tfirougli (3.9) that we should understand the last equality in (3.11) 
without any ambiguity. 

The two-form oj is obviously skew-symmetric, and it follows from its coordinate ex- 
pression that its matrix representation is non-singulaij^ Hence, to defines the symplectic 
structure of phase-space. Since it is a bi-linear non-degenerate form, it provides us with 
an isomorphism between TV and T*V. In particular, given any smooth function H on V, 
the symplectic form uj takes the one-form 

-dH = -- — dpa - -7^ — dx" 

Opa OX°- 

dH d ^ . dH d ^ . 

= T^lT^'^Pb ^ d^ - 7^7^'^Pb ^ dx" 

OPa OX°- OXa Op°^ 

'•{dpb^dx^) = uJ{XH), (3.12) 



^ dpa dx°- dx"- dpa 
where J denotes all contractions, into the Hamiltonian vector field defined by 

Xh^——-— — . (3.13) 

dpa dx°- dx"- dpa ' 



The vector field (3.13) generates a one-parameter family of diffeomorphisms of phase- 
space. The dynamics of a system whose Hamiltonian function is H is given by the integral 
curves of Xh ■ 

A fluid example: vorticity and structure 

Let us make a pause in our present argument and note how the previous construction is 
immediately related to our forthcoming discussion of hydrodynamics [c.f. Chapter pi. 

Consider a single constituent fluid whose elements are described by the conserved 
particle flux n"". Its conjugate momentum Ha can be related to the chemical potential 
(J. through the Euler relation 

P + P = -nVa, (3.14) 

where p and P are the energy density and pressure of the fluid. 

We can express the particle number density current in terms of a family of co-moving 
observers n° = nti° whose normalised four- velocity u'^Ua = — 1 deflnes a congruence para- 
metrised by an affine parameter A. Let us denote the four velocity as the rate of change of 
the coordinate functions of spacetime along these curves, n" = x", where the dot denotes 



''a matrix is non-singular if its determinant is not zero. 
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change with respect to the parameter. The canonical equations of the system are [T^ 

(3.15) 



n 



dnx"' dH dfia u 

— TT- = 7^ — and — — = ^ia-bn 
a\ o^a dA 



dH 

dx"'' 



where H is the 'free-fall' Hamiltonian (3.3). 



Let us contract these equations for the fluid with the particle number density current 



n 



dH dH 



—n 
dH dH 



Ha-hn n 



dx°- dfia dfia dx'^ 



-fia;bn n 



(3.16a) 
(3.16b) 



We see that (3.16a) simply gives the information about the number density of the fluid 



particles, while (3.16b) tells us that the contraction of the momentum equation with the 



particle flux is identically zero. Thus, it follows that the 'forces' acting on the fluid are 
orthogonal to its four-velocity. Hence, the dynamics of this simple system is given by the 
momentum equation 

2^[a;6p'' = -Ha. (3.17) 



The LHS of equation (3.17) allows us to identify the fluid's vorticity with the symplectic 
two-form 

fia;bdx^ A dx" = dfia A dx" . (3.18) 

In this example the only way in which the fluid particles can change their energy is by 
changing their masses (during collisions say). If we do not allow for particle creation or 



mass shifts, the dynamics is constrained to one of the surfaces Vm of (3.5) defined by the 
equation H = m. This is normally known as the Hamiltonian constraint and is a crucial 
ingredient in Hamilton's variational principle. 



Notice that the vanishing in the contraction of (3.17) with the flux n" is equivalent to 



the (trivial) vanishing of uj{Xh,Xh)- We shall return to this point in Chapter^ when 
the constraint surface will have more degrees of freedom. 

Motion in phase-space 

The example above gives us some insight into how the classical mechanics of conservative 
systems is formulated. Without much ado, Hamilton's variational principle can be stated 
in the following manneiQ 



*Up to this point we have only considered states defined on the relativistic single-particle phase-space. 
In the case of 'larger' systems, the geometric construction is completely analogous. We just need to 'count' 
the correct number of degrees of freedom. In such a setting, the spacetime manifold A4 is simply replaced 



42 



3.1. Synaplectic structure of phase space Relativistic kinetic theory 

The single-particle variational principle. Let a be the canonical one-form 
of the relativistic single-particle phase-space and let 7 C Vm be a curve whose 
restriction 7 C 7W connects two events, xo,a;i G M^. We say 7 is a physical 
motion if 7 extremises the action 



Sb]= [a. 



This principle is an empirical fact and constitutes the physical basis of classical mech- 
anics. Equivalently, we can use Stokes' theorem to state the same principle in terms of the 
symplectic form in the following way. Consider one of the hyper-surfaces Vm- The restric- 
tion of the canonical one- form allows us to define a two- form (Jo. Since Vm is a hyper-surface 
of V of co-dimension one, Cj is necessarily degenerate, i.e. the bi-linear form uj has a null 
eigenvectoijjX. The empirical fact in this case states that the physical motions correspond 
to the projected orbits of X into spacetime. 

Let us denote by Qh the one dimensional sub-group of the diffeomorphisms of phase- 
space generated by Xh- If the mass of every individual test particle is conserved, Xh is 
tangent to each Vm- We can think of the 6-dimensional quotient manifold Vm/Gn as the 
non-relativistic phase space of classical mechanics for each inertial observer. 

In order to construct a differential operator which encodes the evolution of a dynamical 



system, we can substitute the definition of the free-fall Hamiltonian (3.3) into (3.13) to 
obtain 

L ^ ui.,XH) = {-,//} = A.^ - /V^O.Pe J^. (3.19) 



Here L? are the Christoffel symbols associated with the metric g . We call (3.19) the 



Liouville operator. This differential operator acts on functions defined on V. Thus, if 
G : V ^M. and L[G] = 0, then G is a constant of the motion. 

Notice that, equivalently, the Liouville operator is simply a one form on T*V defined 
by the canonical two- form lo and the Hamiltonian vector field Xjj- This highlights the 
leading role played by u in the geometric description of a dynamical system. Furthermore, 
as we will show, the symplectic form acts as a building block for the phase-space volume 
measure. 



by the space of kinematic variables; the configuration manifold C whose coordinates are typically denoted 
by q. The relativistic phase-space is then the future, non-spacelike part of T*C with coordinates {q°',Pa). 
Here the indices run over AN values, i.e. we have {q^,q^, ..., g*^,pi,p2, ...,P4jv). 

^A null eigenvector of a bilinear form is a non-vanishing vector with zero eigenvalue. 
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3.1.2 Volume elements 

In this sub-section, we try to keep a completely covariant, coordinate- free language. How- 
ever, we will make a clear connection with the coordinate language used in most relativistic 
fluid dynamics work. 

We need to construct a set of measures. One for the base manifold and one for the fibre 
through each point of the base manifold and also a prescription to transform the volume 
form as we move along spacetime. 

The volume measure in Ai is 

1] = ^/^dx^ A dx^ A dx^ A da;^ = -eabcddx''dx''dx''dx'^ (3.20) 



4! 
The cotangent space of Ai at x, TxAi, has its own volume given by 



TT 



-g)^d/ Adp^ Adp2 Ad/, (3.21) 



this is the volume element in 4- momentum space at each point. 

We can now construct a volume form for V in terms of products of the two-form oj 

n = oj Auj Aoj Auj 

= [dpo A dx°] A • • • A [dp3 A dx^] 

= - [dx° A • • • A dx^] A [dpo A • • • A dps] 

= - [dxO A • • • A dx^'] A [d(5oaoP"") A • • • d{gsa,y)] 

= -goao9ia^92a293a-, [dx° A • • • A dx^] A [dp"" A • • • dp"^] 

= ?7A7r. (3.22) 

Now we can show a result of central importance in the development of the second law of 
thermodynamics. We can ask ourselves what happens to the volume form we just defined 
as it is dragged along the flow of Xh- From the Cartan identity for Lie differentiation of 
r-forms we see 

£xh^ = Xh\ doj + d {Xh\ w) = (3.23) 

and as a consequence 

Cxh^ = (3.24) 

Since 8— forms in V are a one dimensional vector space and Cx^ is an 8— form for some 
vector field X (not necessarily Hamiltonian), it must be proportional to O. The divergence 
of X is defined in such manner as the constant of proportionality between the Lie-dragged 
volume form and itself. Thus we see that any Hamiltonian vector field is divergence-free, 
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that is, the phase flow defined by such field is incompressible! This is precisely what the 
Liouville's theorem states. 

Before we move on, we need an extra measure. Consider the 7— form given by 

n^ =XH\^ = -4:dH AloAojAuj. (3.25) 

The hyper-surfaces of the phase space V which will be of special interest are those for 



which (3.25) reduces to 

n^ = p^oa A vr (3.26) 



where 



oa = ^eafccddxMx^dx'^. (3.27) 



We denote such hyper-surfaces by S. One can verify that the 7— form (3.25) is closed 
and it is also dragged along the phase flow 

Cxh^^ = 0. (3.28) 

We can consider Q!^ as the volume measure on the space of phase orbits, since Liouville's 
theorem guarantees that every cross section of a tube of phase orbit have the same volume, 
thus it can be regarded as a measure of the number of orbits in each tube. Physically the 
integration of J7'' over a hyper-surface S counts the number of states intersecting E [26] . 

3.1.3 Individual mass-shell volume measure 

Finally, we will be interested in cases where the mass distribution is discrete, i.e. we will 
be interested in different mass shells of definite mass m, Vm- We observe that another 



way of interpreting (3.24) is to say that Q. is an L— invariant measure on V and from that 
we see that h{p'^)U, is too, for any function h. To get rid of the contribution from all the 
irrelevant mass shells consider the function 

/i(p2) = 2He{p)5{p'^ + m^) (3.29) 

where He is the Heaviside function which is equal to one if p is future directed and zero 
otherwise. Note that upon integration with respect to Q this will only take into account 
the contributions from the future part of the mass shell Vm- Using coordinates {x^,p'^,p^) 
on V and 'integrating' h{p'^)Q over p^ we obtain 

^m = ?? A VTm (3.30) 
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-^dp^ A dp2 A d/ (3.31) 



on Vm, where 

3.2 Boltzmann's equation 

In the rest of this chapter we consider a spacetime region and assume that it is filled with 
a system of many particles. Following [26], we suppose that the interaction between the 
particles can be divided into long range and short range forces. Then it seems reasonable 
to account for the long range forces in terms of a mean field, in which each particle moves 
like a test particle except when they undergo a collision, which correspond to the short 
range interactions. The short range forces can be accounted for in terms of idealised point 
collisions, whose probability of occurrence is governed by cross-sections taken from special- 
relativistic scattering theorjj^ We also assume here that the collisions only occur between 
pairs of particles, thus, we are considering the case of a dilute gas. 

In this work, the only long range interaction we consider is gravity. If we take the 
spacetime metric as a mean gravitational field, then the associated phase flow describes all 
possible free motions of particles between collisions. 

The assumptions of kinetic theory are physically plausible for dilute gases in which the 
mean free path between successive short range interactions is much larger than the range 
of those interactions. In what follows, we generalise in a covariant manner Boltzmann's 
programme of a one-particle phase space. 

A gas state is represented as a collection of phase orbits corresponding to the free fall 
states occupied by particles between collisions. A collision corresponds to jumps of the 
participating particles from old to new phase orbits. During a collision, particle has no 
phase orbit assigned to it. The state of a gas in phase space is fully characterised by 
the occupation number, N = N[T,], which assigns to a given compact oriented spacetime 
hyper-surface S the number of occupied orbits intersecting it. An individual gas state is 
too complicated to explore in detail. In practice, one uses a Gibbs-ensemble, which is a 
collection of such individual states representing a macro-state of the gas, and studies the 
average distribution of states. We assume that the average distribution of occupied states 
can be described by a continuous density function with respect to the volume measure fi'' 

(iV[S]) = / hn\ (3.32) 

We call /e the one particle distribution function of the gas. It is a Lorentz invariant 



^This is the so-called 'molecular chaos' hypothesis. It is the only time asymmetric assumption of kinetic 
theory. A deeper understanding of this point is required to explain the origin of the asymmetry of the 
second law, in particular the problem of the 'arrow of time'. 
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representing the density of particles in a spacetime region of 3- volume d^x and momenta in 
the range d^p with energies in dE. Throughout the rest of this work we will work with the 
average distribution of occupied states and will drop the angle brackets. Also, for small 
enough regions, the distribution function at a point is independent of the hyper-surface we 
use and, therefore, we will drop the subscript for it. 

Consider a compact region of phase-space D with boundary dD. During a collision 
inside D particles exchange momentum. If we count creations in the positive sense and 
annihilations in the negative, then the total number of collisions inside D can be expressed 
by an analogous distribution 



N[dD] = / gn (3.33) 

Jd 

where g is interpreted as a collision density in D. 

Assuming that / is continuously differentiable, the net number of collisions in D is 

given by 



N[dD] = / fn^ (3.34) 

JdD 

and using Stokes' theorem, the product rule for exterior differentiation and the definition 
of r^^ we obtain 

N[dD] = [ fn^= [ d{fn^) = [ dfAn^= /" d/ a (Xnin) = [ L{m (3.35) 

JdD Jd Jd Jd Jd 

Thus, we have obtained the covariant version of Boltzmann's equation 

Lif) = g. (3.36) 

Specific forms of the collision term g depend on the model under consideration. In partic- 
ular, it is through the mathematical form of this term that one can materialise the time 
asymmetric molecular chaos and the dilute gas assumptions of kinetic theory, i.e. that two 
un-correlated particles before a collision will be correlated for all times afterwards, even if 
they never meet again. We will not go further with a discussion relating to the solutions of 
Boltzmann's equation. We just need to keep in mind that the 'irreversibility hypotheses' 



for kinetic theory are manifest through (3.36). 



3.3 Entropy 

The central notion in any discussion about irreversibility is that of entropy. Unlike most 
quantities that we have encountered so far, for which we have a clear and definite intuition, 
entropy will be the most elaborate and abstract concept we will work with during the 
course of the present work. It is not the intention here to repeat the philosophical meaning 
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of the concept - beautifully addressed by Penrose [60] and Price |63] - but to obtain an 
operational definition to support our future developments on the relativistic dynamics of 
non-equilibrium processes. The material presented in this section is mostly referential. 
The construction follows that of Synge |71j and Ehlers [26] . 

Consider a box divided in A^ equal cells and a set of u marbles. Suppose we accom- 
modate the marbles in the cells in some random manner and then ask how random our 
arrangement is. This may sound terribly vague so let us first compare two possible ar- 
rangements and ask which one is more likely to be the outcome of a random assignment. 
To do this, let us denote each arrangement by an string of numbers 

A = [z^i,J^2,--- ,'^n], (3.37a) 

A' = [u[,u'^,--- ,<]■ (3.37b) 

Here Vi and v'- represent the number of marbles in the ith box for each distribution. 
Obviously, the total number of marbles does not depend on any particular arrangement 



E^^ = E^*' = ^- (3.38) 



i=l i=l 



In a truly random assignment, each box has an equal probability of 1/n of receiving 
a marble. It is not difficult to check that the probability to obtain any of the outcomes 



(3.37) is given by an expression of the form 



P[A] = ^^ -J-Y ■ (3.39) 

If the number of marbles and cells is not very large, the answer follows from the direct 



computation of (3.39) for each array. Indeed, if the numbers are "small enough" one can 



answer such question by mere inspection. However, as the numbers get larger, equation 



(3.39) is not very useful to compare the likelihood of the outcome A given A' . To make 



this expression more manageable, we work with its logarithm 

n 

log P[A\ = log z/! — log JT f j ! — i^ log n 

i=l 
n 

= log!/! — NJlogi^i! — z^logn. (3.40) 



1=1 



This will be very useful in the case where the number of cells, the number of marbles and 
the number of marbles in each cell are very large. In such situation, we can approximate 
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the logarithms in (3.40) by Stirhng's formula 



logo! ~ a log a — a. 



(3.41) 



Thus, using (3.38), we obtain 

n 

log P[A] = V log z^ — / Vi log Vi — V log 



n. 



(3.42) 



%=\ 



Now it is straightforward to compare the probabilities to obtain each outcome just by 
taking the difference of their logarithms 



logP[^] -logP[^'] = log 



P{A!\ 



c-c, 



(3.43) 



where we have made the abbreviation 

n n 

( = -^Uilogi^i and C' = -"^i^i^ogu-. 



(3.44) 



=1 



i=l 



Therefore, the desired relation for both probabilities follows from (3.43) and (3.44). That 
is 



oc 


=^ 


P[A] > P[Al 


c = c' 


=^ 


P[A] = P[Al 


c<c' 


=^ 


P[A] < P[A']. 



Thus we see that the quantity C measures the 'randomness' of a given array A. 

In the example above we were just interested in the 'likelihood' of finding different 
static configurations of very large samples of balls grouped in cells. We can extend the 
analysis to the case where we are also interested in the motions of the balls inside each cell 
(provided there is a mechanism that governs such motion). In that case, the distribution 
of balls will not just depend on their positions within the cells, but also on their velocities. 
In the limit when the marbles and the cells are arbitrarily small, we are essentially looking 
at a gas in a bo:j|j Now, we have seen in the previous section that the total number of 



particles at a given time is given by the occupation functional (3.32). Here an instantaneous 



configuration is only defined with respect to a class of observers; those sharing the same 



^At this level we are not yet concerned about the particular details of the gas. For our purposes it suffices 
that each marble/particle is indistinguishable from one another. This corresponds to the "hard-spheres" 
approximation. 
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simultaneity hyper-surface. As before, let S be a 7-dimensional hyper-surface of phase 
spacq^ The total number of particles and their respective momentum is given by 



Pain = f InPa^l 



(3.45a) 
(3.45b) 



In analogy with the total number of marbles (3.38), equations (3.45) do not depend on 



the particular distribution function of the particles; their specific arrangement within the 
volume Qn- Therefore, to be able to compare the degree of randomness among different 



distributions, we follow (3.44) and define 



^[S,/n] 



(/n log /n) 



(3.46) 



Notice that, unlike the total number of particles and the momentum on a given hyper- 
surface, the quantity H does depend on the particular distribution function. Thus, the 
quantity "^[2, /n] can be regarded as a measure of the "randomness" of the distribution /n 
at a specific moment characterised by the hyper-surface S. 

Consider a space-like hyper-surface Sq. If we want to know what is the most random 
distribution of a very large number of particles N with total momentum P" at a given 



"time" represented by Sq, we just need to extremise (3.46) with respect to all the possible 



distributions. A crucial thing to be noticed is that, without any restriction, such question 
does not make sense. It would be like asking what is the most random distribution of 
marbles in a box without saying how many marbles there are in the first place! In our 
case the relevant question is which is the most random distribution over Sq of an array 



of A^ particles (3.45a) with momentum (3.45b). To answer this question, we introduce the 



Lagrange multipliers «„ and /3a and work instead with 



% = % + a,,N + /3"P„, 



(3.47) 



where fS"" is assumed to be a timelike vector field. Then, we look for the distribution 
function that makes % an extremum 



5n= [ [-(log /n + 1) + an + /3V] 5fn^i = 0- 

Jt.0 



(3.48) 



Note that the one-particle phase space is an S-dimensional manifold. In this sense S leaves one degree 
of freedom. If the projection of E down to spacetime is a spacelike hyper-surface G, then ( |3.45a[ ) effectively 
counts the number of particles at a given time. 
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We can immediately read off from the integrand the most random distribution of N 
identical particles with averaged momentum P" over the hyper-surface Sq to be 

/o = exp [an(x) + /3"(x)pa - 1] . (3.49) 

We say that a system is thermalised at a certain moment characterised by Sq if it 



is described by the distribution function (3.49). Note that the choice of hyper-surface 



is crucial through the whole argument, as it represent a surface of simultaneity for a 
family of observers. However, just as in our previous example of marbles in a box, if each 
box has an equal probability of getting a ball, the result is independent of the particular 
geometry. Therefore, the next question one could ask is how to compare distributions 
between different hyper-surfaces, provided there is a manner of passing from one to another. 

3.3.1 H-theorem 

Now that we know the distribution function which extremises the randomness functional 
^[So,/n] of a gas described by /n at a given time Sq, we are interested to know how to 
relate such measure from one hyper-surface to another. A precise and objective answer 
to this problem is far from trivial. In a rather speculative sense, this is the definition of 
irreversibility; the fundamental notion of the arrow of time. Again, it is not our aim to 
pursue such quest in any depth, but to obtain a clear characterisation of the irreversibility 
of relativistic processes. 

Consider the flow generated by the vector field X (not necessarily Hamiltonian) . The 
change of Ti along the flow of X, parametrised by e, is 

d.,.„ _ d 



de Js(e) 

= -J Cxifnioguni] 

= -/[(! + log fn)Cxfn]ni - [ U log fnCxnl (3.50) 

Our first task is to find the conditions under which the quantity H is conserved along 
the flow of X. That is, when two hyper-surfaces have an equally random distribution of 
particles and momenta. To do this, we need to be careful; we need to guarantee that 
we are comparing the right regions of phase space. This fact becomes evident by looking 



at the second term in the RHS of (3.50); we need to take into account evolutions of 
the volume form Q^ along the flow. To keep the discussion at the simplest level, we only 
consider those flows which are volume preserving. Thus, from the volume dragging formula. 
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equation (3.28), we see that if the vector field X is Hamiltonian, the volume element Q,^ is 



preserved. In this case, 7i is constant along the flow of Xh if the distribution function /n 
is also dragged by the flow, i.e. if /„ is constant along the orbits of Xh- 



From the discussion in the previous section we can ask if the distribution function (3.49) 



is dragged by the phase flow generated by H. Consider the Hamiltonian vector field Xh 



and the distribution function (3.49). From the definition of the Liouville operator (3.19) 
we have 

CxJo = XHfo = L[fo]. (3.51) 



It is straightforward to evaluate this expression to obtain 

L[fo] = exp (an + /^V) tn'>a + P'^"''^ Pa. 



iPb 



Substituting (3.52) back into (3.50) we obtain 



<m 



H 



{(l + log/o)L[/o]}o5,- f hloghCxA 

On + /3>a) exp (an + P^Pa) \afPa + /S^^'^'VaPfe] ^i 



0, 



(3.52) 



(3.53) 



where we have made use of the volume preserving condition (3.28). Therefore, the distri 



bution function /o is dragged by the phase flow only if the Lagrange multipliers satisfy the 
conditions 



a. 



and /? 



(a;6) 



0. 



(3.54) 



Note that the choice of a Hamiltonian observer is crucial to make the right comparison of 
the distribution functions that extremise Ti. over two distinct hyper-surfaces. To see this, 
let us recall the marbles in cells problem. A choice of a flow which does not preserve the 
volume element ri^ , would be equivalent to comparing the distributions of marbles in boxes 
of the same volume but with different number of cells. 



Equations (3.54) represent one of the central results of relativistic kinetic theory. They 



are the macroscopic conditions of global thermodynamic equilibrium [22]. In the next 
section we will express these relations in terms of fluid variables. For the time being, let 
us note that in order for the functional % to be dragged by a Hamiltonian observer, a 
necessary condition is that the spacetime admits a timelike Killing vector field, i.e. the 
spacetime must be stationary. 



We have arrived at the conclusion that the macroscopic conditions (3.54) are equivalent 
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to the micro-physics condition expressed by the colHsion-less Boltzmann equation 



L[M = 0. 



(3.55) 



Thus, we see that for any Haniiltonian observer, a thermahsed system in a given hyper- 



surface So [described by the distribution function (3.49)] will be thermalised in a neigh- 
bourhood of So and we say that the system is in equilibrium. A distribution function 



satisfying (3.49), but not (3.54), is called a local equilibrium distribution [2j. Are there 



any other distribution functions which are dragged by the phase flow? The answer to that 
rests at the core of the 'H-theorem. 



"H-theorem. Let /„ be a distribution function describing the instantaneous 
configuration of a system of particles on a hyper-surface S with volume measure 



On defined in terms of a Hamiltonian vector field Xh- If /o satisfies (3.49) over 
a hyper-surface So and /o 7^ /n, then S is in the the past of So- 

This theorem has been proved a number of times in the literature under two funda- 
mental assumptions: when the system is a dilute gas and molecular chaos is present. The 
following argument will show how we need physical input to realise a proof of the above 
mathematical statement. 

The first thing to note is that, implicit in such a statement, there is a convention. 
We have chosen the direction of the future to correspond to "positive" displacements 
of the parameter of the fiow e. Now, consider the change of H along Xh between two 
neighbouring hyper-surfaces S and S' of a system described by the distribution functions 
/n and /n respectively 



n^'jL] 



'H[j:,fn]+^n 

ae 



e. 



H 



Prom (3.50) we can immediately see that if we can guarantee that the inequality 



L[U] > 0, 



(3.56) 



(3.57) 



is satisfied, then the last term in the RHS of (3.56) 



j[i + \ogU]L[U]ni<Q. 



(3.58) 



Therefore, it follows that the value of % decreases in the positive direction of the flow of 
Xh- That is, if e > 

^[S',/;]<7^[S,/n]. (3.59) 
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In particular, the equilibrium distribution function satisfies 



^[So,/o]<^[S,/„] 



(3.60) 



for any distribution function f^. Thus, the equilibrium value TIq = 7i[T,, Jq] corresponds to 
a minimum. The entropy of a system is defined as the negative of the functional Ti. Hence, 



we see that (3.59) implies that entropy increases towards the future of a given observer. 



Notice that the conclusion (3.59) depends on the validity of the inequality (3.57). This is 



a particular version of the fundamental problem of the arrow of time. In most cases this 
result is simply postulated and in the present work we do as well. However, we should 



keep in mind through the rest of this work that (3.57) is as far as we can go within the 
realm of the classical theory of irreversible processes 



3.4 Macroscopic fluid equations 

To close this chapter, we obtain the macroscopic fluid description of the gas. As noted in 
|74j , a fluid model is appropriate if the collision dominance is assumed to imply that the 
gas is sufficiently close to equilibrium to allow the definition of a unique 4- velocity and use 
of a conventional thermodynamic formalism. 

So far we have been considering the averaged occupation number to describe the number 
of particles crossing a given hyper-surface. Now we focus on those hyper-surfaces for which 



we can use the individual mass-shell volume element (3.30) where the integration over 



momentum space is simply (3.31). In this case the occupation functional (3.32) can be 
written as 



mn 



On 



G 



fuP^T^n 



SB 



(3.61) 



where S" is the part of the mass-shell which is contained in S. Written in this way, the 
term in brackets represents the number density current for the n species, 



n 



/nP^VTn, 



s« 



(3.62) 



which integrated over the oriented volume element aa gives the total number of particles 
passing through G with momentum in S". In what follows, if the region of integration 
of phase space T,* is omitted, the integration (3.62) should be over the entire mass-shelM 
For any smooth distribution function /„ > 0, the particle density current n" is a timelike, 
future directed vector field in spacetime. 



In order to guarantee that all the integrals of the type (3.621 exist, we shall always assume that all 



distribution functions vanish sufficiently fast for large energies at infinity on the mass-shell. 
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In a similar manner, the current density associated with the four-momentum distribu- 
tion of the system defines the kinetic stress-energy tensor of particles of the n species 



-lab 



/nPVvTn. 



(3.63) 



These current densities are special cases of moments of the distribution function in 
momentum space. 



To obtain conservation laws, we need to compute the covariant divergence of (3.62) 



and (3.63). This is easily done as an application of Liouville's theorem. Thus, we write 



(see 1261 



n 



/nTTn 



L[fn]TTn 



-lab 



y/npVvrn =yL[/nK7r„. 



(3.64) 
(3.65) 



Notice that for a simple gas, consisting of a single particle species, Einstein's equations 
together with the Bianchi identities imply that the divergence of the energy-momentum 
tensor vanishes, and in this case the distribution function /n satisfies 



L[fn] = 



(3.66) 



and therefore, (3.64) is a continuity equation stating the particle number current n" is 
conserved. 

For a mixture of particles, we need to consider different distribution functions integrated 
over different mass-shells. In this case there might not be a conservation of energy for each 
individual constituent, but again, Bianchi identities impose a balance law 



-lafe 



En 



ah _ V^ 

'A ;6 - 2^ 

A A 



f LA[fAV7rA = 



(3.67) 



where the subscript A denotes each individual species. 

We can extend the reasoning above and construct an entropy current density from 



(3.46) in an analogous way to (3.62) 



S[T.] = I aaS" 
IG 



(3.68) 
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where again G is a spacelike hyper-surface and 

s" = -y(/nlog/nK^n. (3.69) 

This provides us with another route to the second law of thermodynamics, as encoded 



in the ?^-theorem stated in the previous section [see (3.59)]. In this case, if we consider 
two spacelike hyper-surfaces G and G' as part of the boundary of a 'cylindrical' region of 
spacetime and using Stokes's theorem we have 

S'-S= [ s'^aa - I s'^aa = [ s^cja = [ s\ri (3.70) 

Jg' Jg JdD Jd 



where rj is the spacetime four- volume form (3.20). 



It is clear then, that the total entropy of an adiabatically isolated system will not 
decrease in time if and only if 

s% > (3.71) 

This requirement is equivalent to the condition on the distribution function ( |3.57| ). 
There is no direct way to prove this without further assumptions about the dynamics of 
the gas. The asymmetry in the direction of the arrow of time can thus be reinterpreted as 
an asymmetry in the boundary conditions. 

3.4.1 Thermal equilibrium and the Tolman-Ehrenfest effect 

This subsection is entirely based on the work by Ehlers |26j and has been included for 
completeness and future reference. 

For an observer moving with the normalised four velocity given by the particle flux 

u"" = n^/n, (3.72) 



the equilibrium distribution function (3.49) depends only on the projection of the four 
momentum into its own frame, E = —u"'pa. In this case, one can write the densitised 
thermodynamic quantities p, n, s and p, representing energy, particle number, entropy 
densities respectively together with the pressure p as 

p= NEdE, (3.73a) 

Jm 
/•oo 

n= NdE, (3.73b) 

Jm 

s = /3{p + p) — an, (3.73c) 
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where 

P 



1 /"°° 

-/ NiE"^ -m'^)E-^dE. (3.74) 

"J Jm 

These quantities are defined on a surface of equihbrium states whose coordinates are 



Lagrange multiphers (a,/3). One can show that, in addition to (3.73c), they satisfy the 
equation 

dp = (3-^ds + af3-'^dn. (3.75) 



Thus, comparing with the Gibbs relation (2.39) of section 2.2.1, we can interpret the 



Lagrange multiphers as the temperature and relativistic chemical potential by 

T = /3-\ (3.76a) 

/i = a/3"\ (3.76b) 



There is a subtle point to note in the above discussion. The identification (3.76) comes 
from the identification of the equilibrium thermodynamic variables with the Lagrange 
multipliers defined to obtain the equilibrium distribution function. These, however, retain 
their meaning only when the distribution itself is dragged by the vector field defining the 
evolution of the system, i.e. Xh- What happens if we choose some other vector field which 
does not commute with Xh'^ Does a moving body look cold? A step towards the answer 
was found by Tolman. Equilibrium in a region of spacetime demands that the Lagrange 
multipliers satisfy the evolution equation L[f) = 0, i.e. that the combination a + (S^-pa 



should be a constant of the motion. In particular, the second equilibrium condition (3.54) 
implies the existence of a timelike Killing vector in spacetime, that is, thermal equilibrium 
can be attained only in stationary spacetimes. This is quite a remarkable constraint for 
a relativistic description of thermodynamic systems. In particular, this implies that for 
relativistic thermodynamics the temperature T is not spatially constant. To see this, we 
use coordinates adapted to the Killing vector field /3" so that 

gab-o = 0, (3.77) 



and, from (3.76a), the temperature satisfies 

5afe/S"/3' = goo = -/?' = -T-\ (3.78) 

Therefore, we can see that in a stationary gravitational field the equilibrium temperature 
varies according to the law 

T=(-5oo)-'/'. (3.79) 
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This is called the Tolman-Ehrenfest effect. In consequence, a gravitationally red-shifted 
photon preserves the energy with which it was emitted, that is, there is no energy transfer 
while freely propagating in a gravitational field. 

We shall come back to this in more detail in Chapter 5. 

3.4.2 Deviations from thermal equilibrium 

The above discussion gives us a gist of what is yet to come. The Tolman-Ehrenfest effect 
is a purely relativistic effect for systems in thermal equilibrium which shows that the 
Newtonian idea that spatial variations of temperature alone imply a non-vanishing heat 
fiux cannot be extended to the relativistic realm. Therefore, in relativistic settings, a 
non-uniform temperature distribution is not necessarily associated with a non-equilibrium 
situation. 

Within the framework of kinetic theory, we defined thermal equilibrium in terms of a 
distribution function /n and an observeip^ Xh by requiring that, in the spacelike hypersur- 



face locally orthogonal to Xh, the distribution function should be given by (3.49), where 



the Lagrange multipliers, a and /3", satisfy the conditions (3.54). Thus, in the macroscopic 
limit we could associate them with the relativistic chemical potential and the equilibrium 
temperature, respectively. How do we obtain the macroscopic limit of a non-equilibrium 
situation? or, rather, how do we define entropy (and therefore temperature) in a situation 
away from equilibrium? 

We have a partial answer to the first question, and we face a serious problem with the 
second. On the one hand, we can consider an expansion order by order in deviations from 
thermal equilibrium. This is the well known Grad's approach described in detail for the 
relativistic case by Stewart [70] . Such an approach leads to the set of constitutive equations 
of the Israel & Stewart model of relativistic dissipation discussed in the next chapter. On 
the other hand, following such a path, we will inevitably face the problem of interpreting 
the Lagrange multipliers of the equilibrium function /q. This does not represent an issue 
as long as the macroscopic model is one for which the local equilibrium assumption holds. 
That is, in a regime where LIT programme should operate. This is a surprising result 
because the Israel and Stewart approach models systems with finite relaxation times, as 
one would expect from an EIT point of view. 

To stress the physical relevance of the ambiguities arising from any attempt to define 
a non-equilibrium temperature we consider the following thought experiment proposed 
by Casas- Vazquez and Jou [17J. The zeroth-law of thermodynamics, as introduced in 

^''Or a set of inertial observers; those defined through a timeUke vector field Xe which commutes with 
Xh. 
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Chapter 1, can be used to define an ideal thermometer as the objects which measure 
the same temperature throughout a system which is in thermal equilibrium. However, 
away from equilibrium this is no longer the case. Imagine a system with two components; 
matter and radiation, and suppose we have two thermometers; one with perfectly reflecting 
walls and one whose walls are black. The first one will be unaware of the radiation and 
indicate the temperature of the matter only, whilst the latter will be sensitive to both, 
radiation and matter. In an equilibrium situation, both thermometers will produce the 
same reading. However, when we depart from equilibrium, matter and entropy may have 
different temperatures, thus the thermometers would provide us with different readings. 
We will address this particular problem in section 5.3.1 from a dynamical point of view. 
For the time being, the purpose of this discussion is to make the reader aware that there 
is, indeed, a physical issue when dealing with systems away from thermal equilibrium. 



59 



Relativistic theories of dissipation 



This chapter presents a review of various theories of relativistic dissipation. The aim is to 
make evident how the distinct models of relativistic dissipation share the same philosophy 
in their construction, that is, by providing an energy momentum tensor written in a par- 
ticular frame together with a suitable definition for the entropy flow in order to impose 
the second law of thermodynamics. 

We begin this review with a discussion of the so-called 'first-order' theories of dissipa- 
tion. In particular, we work through the construction of Eckart's model - the first extension 
of non-equilibrium thermodynamics to a relativistic framework. In this case, the simple 
definition of the entropy flow leads to Eckart's hypothesis, which is nothing more than the 
relativistic analogue of Fourier's law. Therefore, a thermal disturbance will propagate with 
an unbounded speed, leading to internal inconsistencies in a relativistic setting. For the 
sake of completeness, we mention the Landau & Lifshitz model, which is not fundamentally 
different from Eckart's original proposal. 

Owing to the inconsistencies arising from the first-order theories, a class of 'second- 
order' theories emerged. Here we will introduce two distinct approaches of this kind which 
differ essentially in their thermodynamic assumptions. Firstly, we will discuss the model 
proposed by Israel & Stewart [37]. They generalise the definition of the entropy flow to 
include 'second order' corrections from thermal equilibrium. This is normally interpreted 
as a truncated series on deviations from equilibrium. The second-order terms of such an 
expansion are given by all the possible covariant combinations of the scalars, vectors and 
tensors available in the theory. These additional couplings need to be obtained by external 
means, i.e. by indirect measurement or through a direct kinetic calculation. Interestingly, 
it has been argued by Geroch and Lindblom |32j that a direct measurement of these terms 
may not be possible. However, experiments showing the existence of second-sound in solids 
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or, alternatively, in superfluid helium seem to contradict their result. The imposition of 
the second law in the second-order theories is made in an analogous manner to the Eckart 
'first-order' model. However, in this case it leads to a generalised version of the Cattaneo 



equation [c.f. equation (2.62)] and, hence, to a relaxed propagation of thermal signals. 
There is a subtle point to note here. In the Israel and Stewart programme, one assumes 
the validity of the local equilibrium hypothesis and yet one correctly obtains a finite speed 
for heat propagation. We shall discuss these matters towards the end of section 4.2.1. 

Secondly, we include a discussion of a different point of view sustained by Carter. 
Although Carter's original proposal aimed at a simple way of doing 'off the peg' calcu- 
lations, the discussion presented here will serve two purposes. On the one hand it will 
shed light onto a highly obscure reference on relativistic thermodynamics. On the other, 
such approach constitutes the original motivation for the multifluid approach to relativ- 
istic dissipation presented in this thesis. For completeness and accessibility, we include a 
transcript of Carter's original views in Appendix A. 

Let us emphasise that the various theories presented in this chapter does not exhaust 
all the possibilities for the first and second order theories available. It is simply a collection 
of the most representative and widely used. The inclusion of a wider selection of models, 
we believe, would not contribute in a significant manner to our discussion. However, it 
would be remiss not to mention some other roads that have been taken. Such is the case 
of the work by Geroch and Lindblom on relativistic theories of dissipation of divergence 
type [32j , or the recent re-derivation of a second-order theory from kinetic theory by Koide 
|21j . Finally, the interested reader may wish to consult the review article by Herrera and 



Pavon [21] for a broader discussion about hyperbolic theories of dissipation. 

4.1 First-order theories 

Let us begin this review with the simplest generalization to describe dissipative process in 
a relativistic context. In 1940, Eckart published a collection of articles entitled 'Thermody- 
namics of Irreversible Processes', the third of which was devoted to irreversible processes 
in relativity. This publication set the basic strategy followed by later developments in 
relativistic theories of dissipation. A generic theory of relativistic dissipation consists of at 
least one fluid, whose particle number density is conserved, together with a local energy 
balance written in a manner that allows us to impose the second law of thermodynamics. 
This section is based on Eckart's original work |25j . 
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4.1.1 Eckart's model 

Let n" represent the particle number density flux of a single species fluid system. In the 
previous chapter we saw that we can express the conservation law for n"" by the expression 

n% = 0. (4.1) 

In Eckart's model, one uses this fact to choose a frame in which to carry out all physical 
measurements. Thus, the Eckart frame is deflned by the family of observers moving with 
the normalised four-velocity parallel to the matter fluid 

u'' = — (4.2) 

n 

where n^ = —gabn^n- 

Given an arbitrary vector field F"", we can form the pair 

/ = -gabu'^F', (4.3a) 

r = h\F\ (4.3b) 

where we have introduced the orthogonal projector to the observer's four-velocity 

h% = 6% + u^Ub. (4.4) 

One can easily verify that, indeed, hf^i^u = 0. Therefore, we can decompose any vector 
field F"" at each point of spacetime into its parallel and orthogonal components to uP" 

F" = fu"" + r. (4.5) 

This is a very useful decomposition to describe vector quantities from the point of view 
of an observer moving together with a fiuid. In a similar manner, an arbitrary tensor F"^ 
can be decomposed in terms of its components 

cj) = UaUbF"'', (4.6a) 

r = h^UcF'^, (4.6b) 

</,"'' = h''Ji\F^'^. (4.6c) 

Instead of proving a variational principle from which the conservation of energy fol- 
lows as an identity, most relativistic theories of dissipation postulate that there exists a 
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symmetric energy momentum tensor which satisfies the conservation law 



T"''.,, = 0. 



Then, using the decomposition (4.6), they write it in its general form 



T'^'' = pu^'u^ + 2g('^u'') + ^'^^ 



(4.7) 



{U 



where 



q^Ua = and Tr'^'ua = 



(4.9) 



In this sense, any divergence-free symmetric second rank tensor serves to define a matter 
model if we interpret (with the right units) p as the energy density, q'^ as the transverse 
momentum and tt as the anisotropic stress tensor as measured by an observer moving 
with the particle flux. With this physical interpretation, one can define the internal energy 
e through the relation |33) 

n{e + a) = p (4.10) 

where a is an arbitrary constant. 



The dynamical relations contained in our matter model imposed by the constraint (4.7) 



are easily obtained by contracting the energy-momentum tensor T""^ with the observer's 
four velocity u"" to obtain 



UaT 



ab 



Ua-.hT'"' = 0. 



(4.11) 



Let us begin with the first term in the above expression. It follows from the decomposition 

(4.12) 



(4.8) that 



UaT 



ab 



-pu -q , 



and, from the definition of the internal energy (4.10), it becomes 






lab 



nle + a)u 



+ <l\ 



(4.13) 



Finally, using the conservation of the particle number density fiux, equation (4.1 ), it follows 
that 



and equation (4.13) reduces to 



an" I = 0, 

■.b 



-UaT 



lab 



ne + q ... 



(4.14) 



(4.15) 
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The second term in the LHS of (4.11) can be written simply as 



(4.16) 



Thus, the local form of the energy balance (4.11 ) from the fluid observer's point of view is 
given by 

(4.17) 



ne + g% + q^iia + Ua.bT^"^ = 0. 



Before giving the physical significance of (4.17), a further simplification can be made. 



Note that when the energy-momentum tensor takes the form 



T"'' = {p + p)u''u'' + pg' 



,ab 



(4.18) 



we can define the hydrostatic pressure of the fluid through the trace of the stress tensor 



1 



Therefore, in the general case, this allows us to define the viscous stress tensor 



-)ah 



-7r'^^+p/l'^^ 



(4.19) 



(4.20) 



which is assumed to be a linear function of Ua-fi, in terms of the hydrostatic pressure. It 
follows that the contractions of this viscous stress with the observer's four velocity vanish 



n.P'^^ = 0, PI = 0, 



(4.21) 



and we can obtain a general expression for P 



ab 



^b = X 


9 
nuacubd„, ua.cucd„, 

2h h U(^^.,i) - -h h, Uc-d 


= x 


^{a;b) ^ 2^{a^b) _ \yC ^^ 
o 


= A 


'a"^ + 2n("'u^)" 


; 



(4.22) 



where A is the viscosity coefficient and we have introduced a as the traceless shear tensor 



1 



-.U^'.Mab- 



(4.23) 



C^ab — U{a;b) o ■ 

Let us define the the invariant specific volume v as the inverse of the particle number 
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density n and use it to express the purely spatial part of Ua-b as 

h'^ Ua;b = u .f, = ni). 



(4.24) 



Thus, substituting the above expressions for the viscous stress, equation (4.20), and the 



spatial projection (4.24) back into the local energy balance (4.17) we obtain 



n{e + pv) + q\ + q^iia - Ua.bP = 



(4.25) 



In this form, equation (4.25 ) is completely analogous to the non-relativistic energy balancaM 
It is only the term containing the four-acceleration tt" which has no Newtonian counterpart. 
Formally, it results from the fact that infinitesimal 3-spaces orthogonal to the observer's 
worldline are not parallel to each other, but relatively tipped because of the curvature of 
such a line [26j. One can interpret this as a contribution due to the inertia of heat. 

For a simple fluid, where the internal energy is only a function of the pressure p and 
the specific volume v, there are always two functions, 6 and s say, such that 



de 

dv 



+ P = e 



as 

dv' 



(4.26) 



Hence, we can write the quantity inside the bracket in the first term of the energy balance 



(4.25) in terms of these two functions as 



e-\- pv = 9s. 



(4.27) 



This is simply the Gibbs relation we found in section 2.2.1 [c.f. equation (2.39)]. Thus, 
the two functions 9 and s correspond to the equilibrium temperature and entropy density. 



In terms of these variables, we can re-express (4.25) to get 



ns+{9-'q%a 



abn — l 



9^ 



q''[9-a + 9Ua]+Ua-bP'''9 



(4.28) 



Written in this form, we can impose the second law of thermodynamics on the dynamics 
of the matter model described by the general energy momentum tensor T"* in terms of the 
thermodynamic quantities n, p and s relative to observer's frame in the following manner. 
Consider a vector field s" representing the entropy flux within the fluid [c.f. section 3.4]. 



^We can write the Newtonian energy balance in as 

de 



n— +V-q-(^-V)-v = 



where e is the internal energy, q is the heat flow, !3^ is the total stress and v is the fluid's three-velocity. 
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Prom (4.28) one can infer that 



sn" + -q" 



(4.29) 



Note that this is far from being the only manner in which we could have written the 
entropy density current. It is, however, the simplest for it is 'linear' in departures from 
thermal equilibrium, i.e it is of 'first-order' in the heat flow. In this case, the second law 



of thermodynamics takes the local form [c.f. equation (3.71)] 



ns 



{0~\n;a > 0. 



(4.30) 



As we have argued in the previous chapter, the inequality has to be imposed from the 
outset; it cannot be proven within the dynamical content of the matter model. 



Eckart's hypothesis 

In order to make the divergence of the entropy flux positive definite, the simplest assump- 



tion is that the two terms in the RHS of (4.28) are independent and non-negative 



On the one hand, the first term in the RHS of (4.28) lead Eckart to propose the 
relativistic analogue of Fourier's law, namely 



-Kh''''\9.h + 9uh\. 



(4.31) 



As before, the proportionality scalar k represents the thermal conductivity of the fiuid. 



It is also clear that equation (4.31) is orthogonal to the fluid's four- velocity. The spatial 



projection of the covariant derivative of the function 9 corresponds to the relativistic tem- 
perature gradient. As we have already discussed, there is no Newtonian analogue of the 
acceleration term —h°'^9ui,. Such a term implies an isothermal flow of heat in accelerated 
matter in the direction opposite to the acceleration [25] . We will provide further details on 
this point in the next chapter. Soon it will become clear, that every reasonable relativistic 
generalization of Fourier's law contains this acceleration term as a consequence of the effect 
of the thermal inertia. 



On the other hand, from the definition of the viscous-stress tensor (4.22) and the 
decomposition of Ua-h into a trace, shear, vorticity and acceleration parts 



1 c ^ 

Ua;b = ^U .^hab + (Tab + ^ab ~ UaUb, 



using a symmetry argument and the orthogonality condition (4.21), it follows that 



na:bP' 



•ab 



CJabP' 



•ab 



\a\ 



(4.32) 



(4.33) 
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Hence, we can write the divergence of the entropy flux expHcitly as a quadratic function 
of its sources 

This constitutes the first attempt to generahse Fourier's law to the relativistic regime. 
However, it can be readily seen that this result contains a covariant analogue of the heat 
equation, without altering its parabolic character. Indeed, it was shown later by Hiscock 
and Lindblom that the Eckart model not only suffers from non-causal behaviour, but also 
it possesses unstable modes for thermal perturbations of ordinary matter models [35] . 

4.1.2 Landau & Lifshitz theory 

The theory of relativistic dissipation proposed by Landau and Lifshitz as an alternative to 
Eckart 's model shares essentially the same features as in the original proposal. Therefore, 
there is no need to extend the discussion in this particular direction beyond a few main 
results to keep the discussion complete and self-contained. The assumptions and calcu- 
lations are completely analogous to the ones in the previous section and can be found in 

The main difference this approach has with respect to the one followed by Eckart is the 
choice of the observer's frame. In the case of Eckart's model we used co-moving observers 
to describe the dynamics of the fluid whilst Landau & Lifshitz used a timelike eigenvector 
of the energy momentum tensor of the matter. In both cases, the equations of motion 
are determined by the conservation of the particle number density flow and the vanishing 
of the divergence of the energy-momentum tensor. Thus, the dynamics are completely 
equivalent and, although the specific form of the local energy balance will take a different 
form, the physical content of both approaches is necessarily the same. Crucially, they lead 
to the same ansatz to ensure the positivity of the entropy fiux. Therefore, they also share 
the same pathologies as found by Hiscock and Lindblom. 

4.2 Second-order theories 

This two-part section contains an account of the well known Israel & Stewart second-order 
theory together with the original version of the pioneering work of Carter known as the 
'regular model'. One of the aims of this section is to highlight the main ontological differ- 
ence between the two 'second-order' approaches, namely, the local equilibrium assumption 
for small deviations from thermal equilibrium in the former and, the EIT approach to 
non-equilibrium thermodynamics of relativistic systems of the latter. There is no claim 
of originality of the material presented in this section. However, it has been written in a 
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manner that, we hope, wiU facihtate an understanding of the details in the original sources. 



4.2.1 The Israel &z Stewart second-order theory 

In view of the unsatisfactory results of Eckart's theory, Israel |37j and Stewart ^70j de- 
veloped a new strategy to solve the inconsistencies of the relativistic first-order theories 
of dissipation. Their approach, firmly grounded in relativistic kinetic theory, was named 
the transient relativistic thermodynamics, but it is better known as the second-order the- 
ory of relativistic dissipation. This model follows the same dynamical construction as in 
Eckart's model up to the definition of the entropy density flux, where they included a 
full set of second-order corrections in the entropy sources. This truncated expansion leads 



to a relativistic generalization of Cattaneo's equation (2.62). The following is a detailed 
derivation of such an extension. 



Let us begin by writing the decomposition of the stress-energy tensor [see (4.8)] in a 
slightly different manner 



T'^^ = pu'^u'' + {p + T)h''^ + 2u^''q^^ + T' 



ab 



(4.35) 



The difference from (4.8) lies in the explicit inclusion of a viscous pressure term r. Again 
q"" is the heat flow orthogonal to the matter's four velocity , while r"^ and r are the stresses 



caused by viscosity in the fluid. Here the tensor t satisfies the relations 



= u"-Tab = t\ = T[ab] ■ 



(4.36) 



Thus, it relates to Eckart's definition through 

r°'' = -7r'^'' + (p + r)/^'^^ 



(4.37) 



As before, the equations of motion are given by the conservation laws (4.1) and (4.11) 



and the imposition of the second law is completely analogous to the previous section. 



Prom the modified form of the energy-momentum tensor (4.35) and properties (4.36) of 



the viscous-stress tensor, the divergence of the entropy density current, equation (4.28), 
can be written as 



1 



q''[0;a + OUa] + {Ua;b)T' 



"^^-^ - TU" 



\a^ 



(4.38) 



Here s"" is the entropy current defined in Eckart's model, equation (4.29), the brackets () 
represent the symmetric traceless part of a given second rank tensor, i.e. for a given tensor 
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A 



ab 



(A 



ab/ 



2 " 



2A 



(ab) 



Krih^fA 



ef 



(4.39) 



and, crucially, 6 is the equilibrium temperature of the fluid. It is obtained through the 
assumption that the local equilibrium hypothesis holds for small deviations from thermal 
equilibrium, i.e. that the energy density remains a function of the particle number and 
entropy densities when the system is driven 'slightly' away from thermal equilibrium. This 
means that, implicitly, the heat flux is still considered to be a fast variable in the EIT 
sense (c./. section 2.2.3) and, therefore, it cannot be treated as a state variable. 

In order to satisfy the requirement of the second law, in addition to Eckart's hypothesis 



(4.31), we require that 



r 

ab 






(4.40a) 
(4.40b) 



where ^ and rj are the bulk and shear viscosity coefficients, respectively. Note that in the 
case of the Eckart model we only had the coefficient A. This is because here we have made 
an explicit distinction between the trace and shear parts of the congruence. 



We see from (4.22) that t is proportional to P . Hence, the quadratic form for the 



ab 



divergence of the entropy density current, equation (4.34), becomes 



'" " Ce ^ ke^ ^ 2r]9 - 



(4.41) 



Motivated by the non-causal and unstable behaviour of Eckart's proposal, Israel and 



Stewart proposed to generalise the definition of the entropy density flux (4.29) by including 
a complete set of second order terms, namely 



sn"- + 



1 



Por^ + I3iq'qb + P2Tbcr'''] ^ + ao^ + "i"'"''^ 



(4.42) 



Here the coefficients /3o, /3i, /32, ao ^-nd ai correspond to the different couplings for the 
second-order terms. In particular, these quantities need to be provided by some other 
means, i.e. by direct measurement or directly through kinetic theory, but cannot be 
obtained from an equilibrium equation of state. 



To compute the entropy production s"'.^6, we can split (4.42) into first and second order 



pieces. The divergence of the linear part is given by (4.38) 

. „ 1 



^I ;a' 



:q''[e;a + OUa] + {Ua;b)T 



"^ TU\ 



(4.43) 
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Therefore we only have to compute the divergence of the second-order part 



■S// ;a^ 



2 
+ 



(3or^u^\ ^ / /?igfegV \ ^ (P2Tbc 



^bc^a 



\a 



V 



+ 



aiqbT 



ab 



(4.44) 



Expanding term by term, we have the following factorisations for the first square bracket 



2,, a 



pQT^U' 



Piqbq^u" 



^2Tbcr'>'u- 



(^, ... 



r-aU^Po 



;a 



ab 






P2U''\ , Tab;cU''l32 
Tab + 2 



(4.45a) 
(4.45b) 

(4.45c) 



Whilst the second bracket in (4.44) is given by 



V 



J -.a^r + Ta^^ + \ — \ Q T, 



aiqbT 



ab 



ablbai.fai-] ^ab,„ ^°''«i 



(4.46a) 
(4.46b) 



In order to factorise the last two equations, Israel (see j37j) introduced two extra thermo- 
dynamic coefficients 70 and 71 such that 



70 + 71 = 1- 



(4.47) 



Thus, it follows that (4.46a) and (4.46b) become 



faoTq" 

faiqbT^ 
V 



^\a-^ + ( ^r ) 'Til 

\a 



+ T 






ab 



Qb;a-^ + 



%!! 



+ '7° 



q 70 



^-'"T^y-jJb^^^'. 



(4.48a) 
(4.48b) 
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Finally, collecting all the terms, we obtain the expression for the entropy production; 



s''-J = -r 



-q 



<a + Vn»/3o + ^f^] re-,%ao ^"° 



2\ e 



'Ojo 






laO - T-„aQ 



"0 



t6ji 



1 /P2U' 



Ua:b + 



2\ e 



TabO + Tab;cU''P2 " gfeja"! 



ai 



QbOli 



(4.49) 



In order to satisfy the second law constraint (4.41 ), the simplest - though not the most 



general - choice is to make each term positive. Thus, the Israel and Stewart theory makes 
the following identifications 



T = -i 



u% + T.^aU^o + I (^) re - q%a, - (^) ^ q^Ojo 



kOh 



ab 



^4 + ^" + ^^). «^^--;^«o-(t).-^ti 



+ QbicU^Pi - T\^ai - y-j) r^e-fo 



(4.50a) 



(4.50b) 



r-^ = -2^h'^-h"'lu,.,d + \(^] T,de + T,a-eu'h-qd;cai-{^) qaeii). (4.50c) 



Heat conduction 

For the purely heat conducting case we artificially turn off the viscous terms, i.e. we set 
r = t"'^ = 0. In this case, the energy-momentum tensor reduces to 



r"'' = (/) + p)uV+p(7"'' + 2u 



O-b 4- 'y^Mrfi) 



(4.51) 



and the entropy current is simply given by 



1 ii°- 

s - sn + -q - -^Piqbq ■ 



(4.52) 
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The simplified version of the divergence of the entropy density flux is 



3^1- I ^ f(^^^ \ a > bo 

?;a^ + na + - I — ^ ) QaO + Qa-bU Pi 



(4.53) 



and we see that the quadratic ansatz (4.50b) becomes 



-kOh' 



ab 



, 1 ■ 1 

;67^ +Ub + - 



^xW- 



QbO + Qa^cU^Pi 



(4.54) 



This is the relativistic generalization of the Cattaneo equation given in section 2.2.2. It 
represents the central result of the purely heat conducting Israel & Stewart model and it 
will serve as a point of comparison for the results in the next chapter. Note that the last 



term in (4.54) is essentially a time derivative of the heat transport four- vector. Therefore, 



its coefficient /3i can be interpreted as a 'relaxation time' for the propagation of thermal 
disturbances. The subtle point we mentioned earlier lies in the fact that, although we 
assumed that the local equilibrium hypothesis holdq^ we have obtained a critical time- 
scale on which dissipation operates, a result which follows from an EIT point-of-view. 
Thus, the reader may wonder where did such a time scale come from? The answer is 
simple, we have removed the thermodynamic information of heat as a state variable at the 
expense of a complete set of coupling coefficients for the second order terms in the entropy 
flux. The remarkable feature of this approach is that it correctly reproduces an EIT result 
from the thermodynamic principles of LIT. 

4.2.2 Carter's theory of dissipation 

The final section of this chapter is devoted to the original work by Carter on a relativistic 
heat conduction. Although it was not made explicit at the time, soon we will notice that 
Carter's ideas are rooted in those which gave rise to the the EIT programme. Furthermore, 



he presented an argument in which the Cattaneo equation (2.62) follows directly from a 



Gibbs relation where the heat flux is included as a state variable. As before, we have tried 
to make every calculation explicit and consistent with the method discussed in the previous 
sections. 

Let us now show that Carter's thermodynamic programme belongs to the class of EIT 
theorieq^ Following the same methodology as in the previous sections, let us consider a 



■^The energy density is a function of the particle number and entropy densities alone, p = p{n, s). 

^We can readily see this from his original arguments on these matters. For completeness, we have 
included a transcript of a highly unknown reference where this fact was made explicit first. The reader will 
find a recollection of Carter's original ideas in appendix A. 
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single species fluid. Recalling the decomposition of a general energy momentum tensor 
with respect to the four- velocity of an observers moving with the particles of the fluid, we 



write [see equation (4.8)] 



T"" = pn^u" + 2q'^''u''> + ■k""" . (4.55) 

In addition to the particle number density flow n", let us introduce an entropy density 
current s" - in general not parallel with n". We can decompose such a current by (c./. 
section 3.4) 

s^ = s*u'' + f, (4.56) 

where s* is the entropy density measured in the Eckart (matter) frame, and j" is the 
component of the entropy current transverse to the matter flow. 



As before, the equations of motion are given by the conservation laws (4.1 ) and (4.7) and 
the second law of thermodynamics is expressed by the positivity of the entropy production. 
Thus, we can express the equations of motion together with the second law in terms of the 
projections 

n".^ = n + u^.^n = 0, (4.57a) 

s\=s + u\s + f,>Q. (4.57b) 



Using the decomposition of Ua-fi, equation (4.32), we can define the symmetric part of 

Ua-b as 

1 

C^ab = -T^U^'.Jlab + Cfab- (4.58) 

Thus, the equations of motion imply a local energy balance equivalent to Eckart's expres- 



sion (4.17) 



r^.ftu'^ = p + u'.^t^p + g% + q'iib + TT^'dab = 0. (4.59) 

Up to this point, the construction is identical to Eckart's model. It is here, where we 
need to introduce a thermodynamic assumption to extract the second law from the energy 
balance that Carter's model differs from Eckart's and Israel & Stewart. 

Let us consider the possible functional dependence of the energy on the state variables. 
In the special situation, when the medium is in strict thermal equilibrium, the energy 
density is a function of the particle number and entropy densities alone. However, in a 
general situation, we could also expect p to be a function of the non-equilibrium magnitude 
of the entropy transfer vector, i.e. 

p = p{n,s,j). (4.60) 
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In an ordinary situation, we would consider states which are close to equilibrium. In 
such a case, j'^ is very small quantity. Moreover, local isotropy requires that the dependence 
on j^ of the energy density should be of second order in the neighbourhood of equilibrium, 
that is 

p = p{n,s,0) + ff{f). (4.61) 

Therefore, a general variation of the energy density takes the form 

dp = pdn + 9ds + figdj, (4.62) 

which we recognize as the extended Gibbs relation of section 2.2.3. Thus, Carter's theory of 
dissipation adopts an EIT point of view before imposing the second law of thermodynamics. 
The specific details of how the inclusion of the heat flux as a state variable leads to a 
relativistic version of Cattaneo equation are presented in the next chapter. Carter's regular 
model is a particular case of the general multifluid formalism for relativistic dissipation. 
For completeness, we finish our discussion as in the original version of Carter's manuscript. 
However, this incomplete view does not present us with the final outcome of the theory. 
The reader has been warned. 

In order to take this last step, let us consider the Legendre-transformed energy density 

P = P- Mqj- (4.63) 

Differentiating with respect to the affine parameter of the observer's worldline gives 

^ = H'h + eS+ j''Cu[Pq\b - TabP'^j''- (4.64) 



where we have made the substitution 



p%f = £u [p\f\ - Ua-bP^^f. (4.65) 

Let us further assume that the heat fiow is simply proportional to the transverse entropy 
current j"", so we can write 

g» = ej^, (4.66) 

where we call 9 the thermodynamic temperature. Notice that we cannot tell if this is the 



temperature that appears in the extended Gibbs relation (4.62). 

In analogy with Eckart's model, collecting the previous result and substituting into the 
local energy balance, one can show that the divergence of the entropy density current can 
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be written as the sum of three terms 



s\0 = -j 



+ 



b + £u [Oub + /ij' 
h"" {fin + es-p)+ i^y - n 



(^ab 



+ f.,b[0-e). (4.67) 

Once again, the simplest way to impose the second law of thermodynamics is to require 
each term to be positive. Thus, we obtain a heat conduction equation of the form 

e.b+CuTT^ = -Yabq''. (4.68) 

Here, the quantity vr^ can be interpreted as momentum conjugate to the entropy density 
flux given by 

7r2 = eua + fi'^, (4.69) 

and the transverse tensor Yab is the thermal resistivity. 

Since the other two terms cannot be made strictly positive, it is required that they 
both vanish. This imposes the further restriction for the energy momentum tensor 

vr'^^ = /i«*(^n + 9s-p) + /x^j^ (4.70) 

and allows us to conclude that the proportionality coefficient in our assumption for the 
heat flow is indeed the thermostatic temperature of the fluid 

e = e. (4.71) 

This last equality is also a consequence of the regularity ansatz. It is to severe a constraint 
for the relaxation time of thermal phenomena described by this early model. Indeed, this 
was a crucial part of the argument in the criticism of Hiscock and Olson of the regular 
model [57]. In the next chapter, we construct from scratch a multifluid approach to heat 
conduction based on the ideas presented in this subsection. We will show that the regularity 
ansatz removes the information contained in the equation of state from the relaxation time. 
Therefore, finding a model which would violate local causality was not very difficult. Thus, 
as expected by Hiscock and Olson, if the regularity ansatz is relaxed, we stand a better 
chance to obtain a consistent relativistic theory of heat. 
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In this chapter we construct a theory of two interacting fluids in General Relativity by 
means of an action principle. Such approach was pioneered by Taub [72] and later by 
Friedman &; Schutz [29j and championed by Carter [13] in the so called convective vari- 
ational approach. Here the central role is played by a single scalar function, the master 
function, which we can understand as a generalisation of the Lagrangian for a purely 
conservative model. 

Variational principles for perfect fluids in General Relativity have been widely studied 
in the literature, see for instance [iOJ and references therein. We start by constructing the 
action from a Lagrangian density as a function of all the possible covariant combinations 
of the kinematical variables, namely the flows n"", s"" and the metric gab- By considering 
an unconstrained variation of the action we obtain a general expression for the conjugate 
momenta. From this we realise that in the generic case the momenta are not aligned with 
the flows - we call this effect entrainment. Later we will show that this is an essential 
ingredient to extend the analysis to dissipative systems. 

To obtain the purely conservative equations of motion, we introduce the idea of a 3- 
dimensional material space [131 dS] HSJ |9], whose points represent idealised particles of 
the medium. We must associate a different material space for each particle species at 
each point of space-time. This allows us to impose the current conservation constraints 
by "pulling" the closed volume 3-forms "back" into space-time, whose duals will be the 
conserved particle number density flow and the "conserved piece" of the entropy flow. By 
imposing these constraints we obtain the equations of motion in the form of a force balance 
equation corresponding to the local version of the energy conservation law. 

We proceed with a thermodynamical analysis. There is an obvious choice of reference 
frame to perform such analysis. However, by defining the relevant quantities as space-time 
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scalars, the choice of observers is irrelevant for the theory - but is of central importance 
for the correct interpretation. 

Finally by imposing the second law constraint on the equations of motion we obtain 
an entropy production rate which is in perfect agreement with the phenomenological one 
produced in the previous section. 

5.1 Oldroyd-Carter rheometric (convective) construction 

We start our discussion on canonical heat conduction by constructing the appropriate 
physical quantities. The main idea throughout this presentation is to follow the variational 
formalism as discussed in Chapter 2 while keeping track of the symmetry underlying the 



equations of motion (2.22). 



The variational approach to relativistic thermodynamics, as presented by Carter [T6] . 
builds on the notion of a lower dimensional matter space which represents the local config- 
uration of an extended body in spacetime |26j . This provides an instructive way to identify 
physical quantities which are invariant under the group of local transformations acting on 
spacetime which characterise the evolution of such a 'body'. 

The set of invariant quantities of the motion represent the state of the medium - its 
intrinsic properties. Thus, the matter space description can suitably be applied to the 
study of equations of state, as discussed in [lOJ. We use the matter space to construct 
the notion of Lagrangian variation in order to obtain the equations of motion of a system 
whose state is represented by conserved particle currents. 

Consider a timelike vector field U whose integral curves generate locally the one- 
parameter family of diffeomorphisms ips £ Diff(A^). For a small enough region 0, of 
spacetime M., these curves can be thought of as the particles of a medium, rigid or elastic, 
or even a fluid. Letting Gjj = Diff (A^)|[/ denote the one-dimensional sub-group of Diff (A^) 
generated by U, we define the matter space as the quotient manifold 

^{n,u) = n/gu, (5.1) 

whose points are the equivalence classes defined by the orbits of Qu. Since the region Q 
is four-dimensional and the orbits of Gjj are the one-dimensional integral curves of f/, it 
follows that the matter space ,^ is a three dimensional manifold |58j . Following Carter 
|10j . we call Q the rheometric group. 

Notice that we can construct as many matter spaces as we have timelike vector fields. 
Thus, anticipating the forthcoming discussion, a two-fluid system generates rheometrical 
groups, Gn and Qs say, acting over the same region Q, and therefore two independent matter 
spaces ^ri and ^s [Sj. 
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Although the matter space is a manifold of its own, its definition implies a projection 
map 

f : x" ^ X^ix"), A = 1,2, 3, (5.2) 

which takes an event x" G il into its equivalence class X G ^/Gu defined by the unique 
worldline passing though it. We can think of X as local coordinates for the points in the 
matter space. 

Quotient manifolds may have a very limited range of geometric structures defined on 
them. For instance, we cannot expect to have an unambiguously defined metric over ^, in 
fact, only in the case of a rigid congruencaj is it possible to construct such a metric [26]. 
However, our main physical ansatz is to assume that a canonical measure is well defined. 
Let us denote by n such a measure which, in the above coordinate representation, can be 
written as 

n^ = NABcdX^ A dX^ A dX^ . (5.3) 



The interpretation of this volume form is obtained by means of the projection map ( |5.2[ ). 
The map / naturally induces the pull-back /* from each cotangent space Tft^\^ to Tx 
and the push-forward /* from the tangent space at T^^il onto Ttu)^- Two remarks are in 
place [26]: 

1. At each point x G 0, the restriction of /=„ to the subspace of Tx^ orthogonal to U is 
a bijection. 

2. The range of /* consists of the 1-forms which are annihilated by U . 

The above remarks suggest that there is natural correspondence between [/-orthogonal 
objects - tensors whose contraction with u"" vanish - and tensors in the matter space. This 
is indeed the case, as shown by Oldroyd [56] and Carter [TSj [l5l [TO] . However, the lack of a 
unique metric tensor field in SS allows the map to distinguish between raised and lowered 
indiceo 



The pull-back of the canonical volume form (5.3) is the spacetime 3- form 



/*n^ = NABcX\X^i,X^.^^dx'' A dx^ A dx^ = Uabcdx" A dx^ A dx^ = n. (5.4) 

^We call a congruence rigid if its expansion and shear are identically zero. 

^For a general non-rigid congruence, we can use the projection map to push forward the spacetime metric 
(?"'' at each spacelike hypersurface orthogonal to U to obtain a family of tensor fields on B. For an elastic 
body, each individual tensor of the family defines a 'state of strain'. We can extend this notion to the case 
of a fiuid. However, we do not have a priori a unique manner to identify one-forms and vector fields on B. 
In the case of an elastic material, the identification is made, conventionally, through the pushed-forward 
metric at the 'instant' the material was formed. For a detailed discussion about the state of strain of an 
elastic body see |26| . 
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By virtue of the definition of ^ and the compatibihty between pull-backs and exterior 



derivatives, d/*n = /*dn , we see that (5.4) defines a closed form 



dn = /*dn^ = ^y^ n^ X^aX\X^f,X^^ dx'^ A dx"" A dx^ A dx^ = 0. (5.5) 

oX 1 1 I I 

It is easy to see that this is equivalent to some conservation law obtained by applying 



the dual operator given by the spacetime metric to (5.5), 



* dn = ^n"''^''^€„feed = o n[abc;d] = 0, (5-6) 



which follows from the same symmetry reasoning. Thus, equation (5.6) defines the con- 
served current 

n% = where n" = -Ubcde"'""^ = nu", and u^Ua = -1. (5.7) 

We can interpret the canonical measure n on the matter space ^ as the number 
density n of conserved particles whose world-lines are a subset of the flow generated by 
the normalised vector field f/ in fi. We can say, in a less formal manner, that U generates 
the conserved particle flow. 

It may seem that this construction obscures a fact that we could have taken for granted 
as a starting point. However, the true value of using a matter space to describe the motion 
of an elastic or fluid medium will shortly become evident when we use the relevant field 
variations in the action principle. 

Notice that functions (0-forms) defined on ^, by construction, are invariants of the 
flow generated by U. Conversely, any function on spacetime which is invariant under the 
action of Qjj has a natural counterpart in ,^. Thus, there is a one-to-one correspondence 
between ^jy-invariant functions in 0, (functions on spacetime which are preserved by the 
flow) and arbitrary functions in ^ [58j. In particular, the coordinate functions X are 



dragged by the flow generated by the current (5.7). We can easily see this, by the explicit 
calculation 

CnX^ = X^^n- = X^.^ribcde''''^ 

= ^7VscD6"''^'^X^,X^,X^,X^, = 0, (5.8) 

where the last equality follows from the commutative property of partial derivatives and 
the anti-symmetry of e. The study of more general rheometric invariant objects is described 
in detail in ^^. Our interest here is to use the matter space to assist in the construction 
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of an action principle for thermodynamics following the reasoning described in chapter 2. 

5.1.1 Lagrangian and Eulerian variations 

In the construction of the variational principle for general relativity, we used a one- 
parameter family of diffeomorphisms ips locally generated by an arbitrary vector field 
^, acting on spacetime fields - including the metric - at fixed points in spacetime whose 
coordinate representation plays no role. In the forthcoming exposition of a variational 
principle for relativistic thermodynamics, however, we want to follow Carter's work |15j 
and introduce two different types of variation. 

The first type of variation we wish to consider corresponds to those whose effect on 
the fields in Q is described in terms of a coordinate system which is itself dragged by the 
displacement vector ^° = Sx"". We call such variations Lagrangian and denote them by 5i. 
In particular, as we will explain, Lagrangian variations leave the projection map / fixed 

unms]. 

The second kind arise whenever we have a definite prescription for identifying the 
points of spacetime and we consider a variation of the fields at a fixed spacetime point. 
Such a prescription is, necessarily, not generally covariant since, as discussed in Section 



2.1.2, it gives a physical meaning to the spacetime point at which the variation is taken. In 
this case, a general perturbation cannot be required to preserve the projection map and, 
therefore, we need to make a distinction between this Eulerian type of variations from 
their Lagrangian counterpart. Let us denote this type of variations by 6e. 

For any field quantity ^ and Lagrangian displacement ^ generating the one-parameter 



family of diffeomorphisms (fs, as in (2.5), we have that the difference between Lagrangian 



and Eulerian variations is the Lie-derivative \15\ ITOl [T6l [6] 

(de - 6e)^ = C^^. (5.9) 

The explanation for this expression - found in many treatments of relativistic elasticity 
- is sketched below by making a distinction between active and passive diffeomorphisms. 



The spirit of the argument is to facilitate the interpretation of (5.9) and not to provide a 
rigorous proof for it. 

Active and passive diffeomorphisms 

Although in the course of this work we have been using extensively the idea of diffeomorph- 
isms of spacetime, we wish to address a subtle but relevant point; the difference between 
those which we refer as active from the ones we call passive [66j. 
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An active diffeomorphism ip is a differentiable map with differentiable inverse from Ai 
into itself. This is to be distinguished from a passive diffeomorphism (j), which is a similar 
map from a coordinate system for an open region of spacetime into another coordinate 
system. A coordinate system for an open subset fi C 7W is also an invertible differentiable 
map assigning to each point in the interior of Q its coordinate representation in M . Thus, 
we have the following mappings: 

1. Active diffeomorphism: (p : q >-^ p = (p{q), where p,q £ Q. 

2. Coordinate system: ^ : p \-^ x = ^{p), where p £ Q, and x G M . 

3. Passive diffeomorphism: (p : x' ^-^ x = </>(x'), where x',x G M . 

For simplicity reasons, we consider the case of a scalar field only. A scalar field ^ is 
a function on spacetime. An active diffeomorphism defines a new scalar field at the same 
spacetime point 

^{q) = ^{ip{q)). (5.10) 

Also, a scalar field ^ on spacetime determines a function t/j on M defined by 



^(x) = ^{i9-^{x)), 



(5.11) 



which through a passive diffeomorphism cp defines another scalar field in a new coordinate 
representation given by 

^\x') = ip (0(x')) . (5.12) 

This can be summarised by the commutative diagram 



Active 



Lagrangian 




(5.13) 



Passive 



Eulerian 
Let us now consider the one-parameter families of active and passive diffeomorphisms 



locally generated by ^ and C, respectively. As in section 2.1.2 let (pQ = idJvi. and 0o = idM"^. 
There are two ways we can compare the different numerical values of the fields ^ and ^. 
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First, we can use the upper part of the diagram to visuahse the setting we discussed in 



section 2.1.2 There, we defined the variation of the field ^ [see equation (2.6)] in terms of 



the one-parameter family ips- Following the top cycle of (5.13) we construct the difference 



of two distinct fields, ^ and ^, at the same spacetime point. 



A^ = ^{q) - ^(q) = ^{ipsiq)) - ^(q). 



(5.14) 



We can divide this variation by the change in the parameter which, when evaluated at 
s = 0, is indeed the Lie-derivative of the field ^ with respect to the vector field S, at the 



point q. Therefore, (5.14) reproduces the RHS of (5.9) 



The other way of comparing the two fields is by taking the lower route of (5.13). 



Starting at some point q € Q, we can look at the difference between the value of the field 
^{q) and its associated value in terms of a coordinate representation ip'[x') after being 
affected by some passive diffeomorphism (pto 



^{q) - V(x) = ilj {{cPt o i^)iq)) - V^(^(g)). 



(5.15) 



Then we can take another point p £ Q and proceed in the opposite direction 



^{p) - ^'{x') = V'((./>-t o ^){p)) - v^'(^(p)). 



(5.16) 



The comparison of the two different paths gives the result 



A^ 



^(q)-V(x)] -[^{p)-i;'{x')] 
^{q)-^{p)] -[i;'{x')-i^{x)]. 



(5.17) 



Notice the similarity between the first term in (5.17) and the intermediate expression 



in (5.14). Despite the suggestive look, they are not the same; they are evaluated at 



different spacetime points. In the limit when the parameter t tends to zero, this defines 
the Lagrangian variation of ^ 



5^^ = ^ ((i?-i o(Pto ^){q)) - ^{p). 



(5.18) 



Now we can explicitly see how the "coordinates are dragged" by the Lagrangian displace- 



ment (^. Looking at the first term of (5.18), we realise that each diffeomorphism in the 



family acts directly on the coordinate system i} associated with the point q. 



The second term in (5.17) is easier to interpret. It is simply the difference between two 

l"^, and, therefore, at the same spacetime 



representations of the field at the same point in 
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point. Thus, this term corresponds to the Eulerian variation 

5e* = il^'ix') - ip{x) = i^iMx')) - tp{x') = ^ {{'&~^ o ^t){x')) - *(t9-1(x')). (5.19) 



Note that in this case the exphcit choice of a coordinate system prevents (5.19) from being 



equal to the general covariant expression (5.14) 



The final point to notice is that both expressions, (5.18) and (5.19), compare values 



of the field ^ at spacetime locations, independent of the coordinate choice, thus (5.17) is 



a generally covariant expression. However, in order to evaluate (5.18) in the limit when 



the parameter t is zero, we need the condition -i? o ^g ° "*? = idTW independently of 
the coordinate system we choose. This bounds the displacement C to the vector field ^ 



generating the family ips- Therefore, in the infinitesimal limit, (5.14) coincides with (5.17) 



and, using the definitions for the Lagrangian and Eulerian variations, (5.18) and (5.19), 



we obtain the expression for the Lie derivative (5.9). 



5.1.2 Eulerian variation of conserved currents 

Prom the previous discussion on the rheometric invariance of the coordinate functions of 
^, an arbitrary variation of X is of purely Eulerian type and is simply given by 



d,X^ = [5e - C^] X' 



-Li-X^ 



-Hp.X' 



-^^a^l- 



(5.20) 



Thus, in general, changes of rheometric invariant objects due to an arbitrary coordinate 
displacement ^ are solely given in terms of displacements orthogonal to the flow. In this 



case, we can freely remove the _L sub-script from (5.20), knowing that the coordinate 



variations parallel to the flow will not contribute to the Eulerian perturbation. 

Let us consider the Lagrangian and Eulerian variations of the objects we have defined 



so far. A Lagrangian variation of the closed 3-form (5.4) is given by a displacement tangent 
to the flow of U 

biVL = C\un = XCun, (5.21) 

where A is some parameter of the flow. From the Cartan identity for Lie-differentiation of 



r-forms [see equation (3.23)] we have 



Cun = d(C/Jn) = d{^''-^\nabcdx'' A dx^ A dx'^). 



The a-contraction on the RHS of (5.22) is 



[Sen]a = diu^-Uabcdx" A dx") = {u'^n'^€abcd);edx'' A dx" A dx 



(5.22) 



(5.23) 
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This follows from our ansatz ( 5.3 ) , which implies the conserved current in the form ( 5.7 ) and 



a symmetry argument. In the same manner, the contractions [(J^n]^ and [(5^n]c are identically 



zero. Thus, the Lagrangian perturbation of the pulled-back volume form (5.21) vanishes, 



as expected from the definition of Lagrangian variations and rheometric invariance. 



The Eulerian variation of n follows from the above discussion and the definition (5.9) 



4n = -C^n = -d(Cjn). (5.24) 



Using the same contraction decomposition as in (5.23), the a contraction is 



[den]a = -{enabc);d<ix'' A dx" A dx" = -{i^nabc + nabc;dC)^^ ^ ^^ ^ dx'^- (5-25) 

It is not difficult to see that by adding up all the contractions, and shuffling the indices in 
a suitable order, we obtain the well known coordinate expression [HI |6] 



^enabc = - \^i ■,[a^\d\bc] + ?^afec;dC j , (5-26) 

where the index between the bars is excluded from the antisymmetrisation bracket. Thus, 



using the definition (5.7) and contracting with the volume form eabcd, whose Eulerian 
variation is [6] 

^ef^abcd = -^(^abcdg'^'^^egrs, (5-27) 

we obtain a key result for our forthcoming variational principle; the Eulerian variations of 
the conserved current 

5^n^ = -L^n'^ - rf (t-,^ + ^^''^eSfec) , (5.28) 

cannot be varied independently of the metric. This is, however, not surprising. An arbit- 
rary variation of the coordinates will not preserve the projection map in general. 

5.2 The two-fluids action principle 

We start this section by considering an unconstrained variation of the action for a system 
of consisting of two fluids. This will exhibit the general formalism contained in a multi- 
fluid variational principle. It will show that, in order to obtain the correct dynamics, we 
will need to impose a set of constraints on the world-line variation of the currents. It will 
also serve as further motivation to use the mathematically more elaborate, but physically 
deeper, assumption of a well defined volume form for two independent matter spaces in 
the case of the study of conserved currents. 
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5.2.1 The unconstrained variation 

Let us define the relativistic invariant Lagrangian-type density A as the integrand for the 



matter sector of the Einstein-Hilbert action (2.4). We refer to this density as the master 



function. The symmetry group of general relativity requires the action to be generally 
covariant, therefore, assuming that the system is isotropic, the master function can only 
depend on the scalars we can form from the two fluxes and the metric. As customary, 
following the notation previously introduced, we use n" and s" to denote the particle 
number and entropy density currents respectively. Thus, we have the three scalars 



n' 


= 


- gabn^n^ 


s' 


= 


- gabs'" s'', 


f 


= 


- gabn^'s^, 



(5.29a) 
(5.29b) 
(5.29c) 



and the matter action [see (2.4)] 



Sm{^) = / A(n,s,j)^/^dO. 
Jo. 



(5.30) 



An unconstrained variation leads to 



,, dK^ dA^ dA^. 
6A = -—6n + -—6s + -—6j. 
on OS Oj 



(5.31) 
Using equations (|5.29|), we can change the passive density variations for dynamical 



variations of the world-lines generated by the fluxes and the metric. That is, 

5n = - —[2gabn''5n'' + n''n''5gab], 



1 



5s = --[2g,bs''5s'' + s''s'6gab], 
Is 

5] = - hgah{n%s'' + s''5n^) + n'^sHgabi 
^3 



together with the fact 



1 (9A _ aA 1 aA _ OA 1 9A _ aA 

2n dn dn^ ' 2s ds ds'^ ' 2j dj dj'^ ' 



(5.32a) 
(5.32b) 
(5.32c) 



(5.33) 
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we find that the variation (|5.3l|) becomes 
SA = 



dA 



„5A 



612" + 
+ 



.9A 

■^ r^ 9 'Sa ri .9 ^O 



dA 

'p' 

dA 



(5s" 



dA ^ ^ dA ^ ^ dA ^ fj 

jjrj^Z fJgZ jjjZ 



Sg, 



ab- 



(5.34) 



From (5.34), we can read off' the conjugate momenta associated with each of the fluxes 



dA 



lia = ^,=gab{B^n' + A-^s') 



^ ~ ds- ~ 
where we have introduced the coefficients 

dA 



QabiB^s' + A'^'n' 



i3° 



drfi 



B^ = -2— A""^ = - — 

Qg2 ' Qj2 



(5.35a) 
(5.35b) 

(5.36) 



The conjugate variables (5.35a) and (5.35b) demonstrate the fundamental role of the 



master function [ 16j . The distinct role of the ffuxes and their conjugate momenta are 
mostly omitted in the ffuids literature. A key advantage of the variational approach is 
that these quantities are immediately determined by the form of the master function. 
Furthermore, it is clear that the momenta are not aligned with the respective currents, 



owing to the fact that the master function depends on the relative flow (5.29c). This is a 



physically relevant effect. Fundamentally, there is no argument to rule out the dependence 
of the master function on the relative ffow. In fact, this coupling is associated with the 
entrainment effect that is known to be of central importance in multi-fluid systems [6]. 

In the case of a super-fluid neutron star core, the entrainment arises owing to the strong 
interaction and couples the neutron and proton ffuxes. In the heat problem we want to 
address, we will show that the entrainment between matter and entropy is associated with 
the thermal relaxation of the medium. 

It remains to identify the coefficient for the metric variation. From the discussion in 



chapter 2, it is related to the energy momentum tensor through (2.17). From the symmetry 



of 5gab and the definitions (5.35a) and (5.35b), we can re-write the variation (5.34) in terms 



of the fiuxes and their momenta to obtain 



5A = fxM" + dcds" + -5 (/^cn" + ecs'')5gab. 



(5.37) 



It is not difficult to see that the metric contraction with each individual quantity inside 
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the brackets in (5.37), iicn"' and 6cS°', is symmetric only when the entrainment coefficient 



^"^ is identically zero 

However, it is evident that these terms satisfy the Noether identity 



(5.38a) 
(5.38b) 

(5.39) 



regardless of entrainment. This should indeed be the case if we expect to produce a stress- 
energy tensor from the variation of the master function with respect to the metric. 



From (5.37) and the variational definition (2.17), we can read off the stress-energy 
tensor 



-tab 



,o,„b 



QO. b 



/i^n" + r s" - A(n, s, g)5 



ab 



(5.40) 



The equations of motion are given by a variational identity [see (2.22)] as the local con 



servation law of the stress-energy tensor (5.40). If the variation of the fluxes, n" and s". 



were left un-constrained, such variational identity would require the momenta (5.35a) and 



(5.35b) to vanish in all cases. To show this, we explicitly evaluate the local conservation 



law, and then note that this implies the orthogonality between the currents and the 'forces' 
given by 

TJ'.^^,n'' = fia-^bn^n" + Oa^bs'^n" + 2n^.(^/x,)n'^ + 2s^.(f,e,)n'^ = 0, (5.41) 

with an analogous expression for Tj'.f^s"'. Even in the simplest case, when both currents 
are aligned, the only possibility for these equations to hold is that both momenta are 
identically zero. In particular, this is an indication that the particle number density flux 
has only three independent components and, therefore, the only compatible dynamics with 
an unconstrained variation is trivial. Thus, the only hope to obtain non-trivial dynamics 



from this type of Lagrangian master function, where equation (5.41) follows as a vari 



ational identity, is to impose some kind of constraints on the variations of the currents. 
In particular, if we assume from the onset that at least one of the fluxes coming into the 
action principle is conserved, we can fully exploit the matter space formalism introduced 
in the previous section. Furthermore, we can use it to explore the rheometric symmetry of 
the system. 



87 



5.2. The two-Ruids action principle 



Variational thermal dynamics 



5.2.2 The constrained variation 

Let us assume that the particle number density and entropy currents are conserved. This 
corresponds to the existence of two different matter spaces together with their respective 
volume forms. This imposes the constrained Eulerian variation of each current in the form 



(5.28). Thus, assuming two different Lagrangian displacements, ^n and ^s, we have 



.^abi 



5.n' = -C^X - n' Ca + :,9'"'5e9ab , 



5,s' = -C.y-sUe, + -g-%gab 



(5.42a) 
(5.42b) 



Some feeling of discomfort may arise from the use of equations (5.42) directly into the 



variation of the master function. The current displacements in (5.37) are of the general form 



(2.6), whereas equations (5.42) correspond to the Eulerian part of the general variation 



(5.9). However, as far as the variation principle is concerned, there is no difference in 



choosing one kind of variation from the other. The reason for this is given by Carter [15j 
in the following argument. Consider the Lagrangian and Eulerian variations of the action 



(5.30) produced by a general displacement ^. Then, using formula (5.9), their difference is 

(5.43) 



given by 



{5, - 6e) s{n) 



I C, [A. 
Jn 



-g dn. 



The integrand is straightforward to evaluate 



A[A7=5] = A;a^^/=5 + Ca^V^ = (Ar );a%/=5, 



(5.44) 



where we have made use of the formula C^y 
theorem, we have 

{6e - 5e) S{n) = 



an 



= -y/^.^"^. Therefore, using Stokes 
Aed^a, (5.45) 



where dS^ denotes the surface volume element of the boundary d^. Thus, from our regular 
assumption, the displacement ^ vanishes on the boundary of the region 0, and, therefore, 
the action principle does not distinguish between Eulerian and Lagrangian variations. 

Before we continue with our variation analysis, let us make one further assumption; we 
consider the metric as a passive field. By this we mean that, to linear order, a Lagrangian 
variation of the metric is zero. Thus, the metric acts as a predetermined background. This 
does not contravene general covariance in the sense that the equations of motion will still 
follow as a Noether identity, but their symmetry may not be as broad as Diff (M) but a 
sub-group of it. 
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Substituting the Eulerian variations (5.42) into equation (5.37), tlie variation of the 
master function becomes 



5eA 



-n 



e^ + 2eiV])^, + (-.Vs + 2e;yi)0, 



+ 



^-n' + e''s'-(^licn' + 9cs')g'"'] d^gab- (5.46) 



It is straightforward to show that the first term in each bracket of (5.46), contracted with 
the corresponding momentum, can be written as 



-n-.^^C^a = - 2/X[,^,]n'' + n%M, C " 2C^'Va + 2en''Va ^ , 



■I'^r,^], 



^[«T,^, 



a i-b/ 



~s\esOa 



'^Ola-Ms' + S\^,9a) C - 2d^"^a + (2^ts''^ea 



(5.47a) 
(5.47b) 



Thus, the constrained variation (5.46) can be written in terms of the two different dis- 



placement vectors for the fluxes, ^n and ^s, together with a general displacement for the 
Eulerian variation of the metric and a divergence term which does not contribute to the 
variational principle. 



Using equation (5.9) and the assumption of a passive metric we have 



Se9ab = -2C(a;b), 



(5.48) 



and after the usual integration by parts, ignoring all the boundary terms, the Eulerian 



variation of the action (5.30) is [compare with (2.8)] 






where we introduce the notation for the 'forces' 



-gdn, 



(5.49) 



fa = '2fJ-[a;b]n^ + "-^/^a, 
fa = 26'[a;fc]S + S .fj9a, 



(5.50a) 
(5.50b) 



and the energy-momentum tensor is found directly from (5.46) to be 



T,' = flan'' + Oas' + ^5„^ 



(5.51) 



where we define the generalised pressure, ^, as a function of the conjugate momenta by 



^ = A-fian^-eaS''. 



(5.52) 
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Using (5.51) and (5.52) only, it is an easy task to show that the divergence of the 



energy-momentum tensor can be written as a "sum of forces" 



J a I J a' 



(5.53) 



Furthermore, as stated in chapter 2, if the matter fields satisfy their own equations of 



motion, that is, if each of the forces (5.50) are identically zero, the local conservation law 



for the energy-momentum tensor follows as an identity [see (2.22)]. In particular, this 
implies the force balance 

ra = -fa- (5.54) 



Thus, the equations of motion are given by the system (5.53) together with with ( |5.54[ ). 
It is worth mentioning that the reasoning followed to obtain those equations also applies 



to the unconstrained variation (5.37). However, the only solution to those is trivial. That 



is because we needed to impose the orthogonality condition with respect to the particle 



number density flux. In the case of the constrained equations of motion (5.53), the ortho' 



gonality follows as an identity consequence of the symmetry of (5.50a) 



f>'' 



2/X[a;fep^^'' 



(5.55) 



where we use the conservation law (5.7) to eliminate the second term of (5.50a). Therefore, 



by virtue of the force balance, equation (5.54), the particle flux is also orthogonal to /^ 
and hence 

(5.56) 






holds in all cases. 

In the strict variational sense, we should have an analogous identity for the entropy 
flux s". This would imply that the particle number density current and the entropy flux 



are parallel, as can be seen by the symmetry of both equations (5.50) when contracted 



with each of the fluxes provided both currents satisfy a conservation law of the type (5.7). 



Furthermore, it means that both currents are dual to two, in general distinct, spacetime 
volume 3-forms pulled back from their respective matter spaces. In this case only, both 
matter spaces coincide and the canonical volume form becomes degenerate. However, as 
far as the variational principle is concerned, the dynamics of the system are implied by the 



vanishing of (5.50) regardless of any production term. It is the identity (5.55) that implies 



the non-trivial behaviour of the system, which follows from the conservation law of one of 
the currents only. This is the key point we analyse in detail in the coming sections. 

Let us make a final remark about the symmetry group of the equations of motion 



implied by the identity (5.56). The orthogonality between the particle number density 



90 



5.3. Thermodynamic interpretation Variational thermal dynamics 

current and the divergence of the stress-energy tensor shows that both forces, /° and f^, 
are [/-orthogonal tensors and, therefore, then evolution has a suitable description in the 
matter space corresponding to the particles. 

5.3 Thermodynamic interpretation 

The action principle and the equations of motion we have derived are those of an arbitrary 
two-fluid system. So far the labels n and s, denoting particle number density current and 
entropy flux, have just been ad-hoc tags to distinguish between two generic fluid species. 
This section is devoted to show that, indeed, these tags can be given the thermodynamic 
meaning they express. 

As noted in the last remark in the previous section, when the particle number density 
current satisfies the conservation law (5.7), the forces (5.50) are [/-orthogonaPl We will 



assume from the outset that this is the only conserved current. This suggests that it is 



natural to choose observers moving with the matter flow [see (5.7)]. This coincides with 
the choice made by Eckart, as we have seen in the previous chapter. In the following 
discussion we use the terms matter and Eckart frame to express the same observer choice. 
More complex settings, e.g. when dealing with more than one conserved current (excluding 
entropy) or when reactions are present, make the choice of frame less obvious. We will 
extend this discussion further at the end of the section. 

From the previous discussion, it is clear that when the entropy flux is not conserved, 
it cannot be parallel to the matter flow. Thus, it can be decomposed into co-moving and 
[/-orthogonal components. The entropy flux is then expressed as 

s'' = s*{u'' + w''), s* = -s''ua, (5.57) 

where s* is the projection of the entropy current into the matter frame and w"" is the 
relative velocity between the two frames which, by definition, satisfies the orthogonality 
condition w"'Ua = 0. 

Let us denote the normalised entropy four velocity by u^. In this frame the entropy 



fiux is written simply as s°' = su^, where s is given by (5.29c). From this and (5.57), we 
see that 

s* = s«na) = s7, where -f = {u" + w^l = {1 - w'^)''^/'^, (5.58) 

is the redshift associated with the relative motion of the two frames. 



Analogously, we can decompose the momenta (5.35) into 'co-moving' and 'orthogonal' 



''Here U corresponds to the normalised four- velocity associated with the particle number density flux. 
See discussion in Section 5.1. 
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pieces. Substituting (5.57) into the momenta (5.35) one obtains 



Oa = {B^S* + ^^^n) Ua + B''s*Wa = 0*Ua + 6^ 



(5.59a) 
(5.59b) 



Here we denote with a 'star' the components of the momenta (5.35a) and (5.35b) ahgned 



with U and with a 'natural' symbol the pieces which are orthogonal. Thus, we have the 
definitions 



//* EE -n>a = ^''n + ^°^s* and ^^=^^'3*, 
e* = -u'^Oa = A^'^n + B^s* and 6^ = B's*. 



(5.60a) 
(5.60b) 



It is the aim of this subsection to show that the 'starred' quantities in (5.60) are the 
ones susceptible to a thermodynamic interpretation. 

When dealing with a generic multi-fluid system, the components of each momentum 
measured in its associated frame - the one defined by the current it is conjugate to - 
correspond to the chemical potential of that species. In the case of the entropy flux, we 



use the tag thermal for its conjugate momentum, equation (5.35b). It is reasonable that 



the projection of the thermal momentum (5.60b) into the matter frame should correspond 



to the temperature an observer moving with the flow would measure. In the following 
subsection, we argue that the quantity 6* indeed satisfies the definition of temperature 
in the Gibbs sense; as the rate of change of the energy with respect to the entropy while 
keeping all other thermodynamic variables fixed. 

5.3.1 The temperature problem 

Thermodynamic properties such as temperature or pressure are uniquely defined only in 
the case when the system under consideration is in equilibrium. Intuitively this makes 
sense since, in order to carry out a measurement, the measuring device must have time 
to reach equilibrium with the system. Such measurement is only meaningful as long as 
the time scale required to obtain a result is shorter than the evolution time of the system. 
This, however, does not prevent a generalization of the various thermodynamic concepts. 
The procedure may not be unique, but one should at least require the generalized concepts 
to be internally consistent within the chosen extended thermodynamics model. 

In our present discussion we are facing a larger problem. We must find an internally 
consistent set of thermodynamic variables to describe the state of the system when it is out 
of equilibrium. This approach should allow us to infer the thermodynamic temperature 
from an extended Gibbs relation in terms of this set of variables. Furthermore, since 
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there are no well established transformation laws for the temperature (see [45j). we need 
to show that when the analysis is carried out in the matter frame, the thermodynamic 



temperature corresponds exactly to the dynamical temperature (5.60b). Therefore, in 



the spirit of physical interpretation, unless stated otherwise, we shall consider all physical 
measurements to be taken in - and with respect to - the matter frame defined by the vector 
field U. We attach a star symbol to distinguish any quantity defined in this manner. 



Let us start by re- writing the variationally defined two-fluid stress-energy tensor (5.51 ) 



into its more common dissipative form [see definitions (4.8) and (4.35)]. Following the 



general decomposition (4.6), the energy-momentum tensor (5.51) takes the form 



Tj> = p*UaU^ + QaU^ + Uaq^ + Pj" . (5.61) 

Here, the energy density p* , the heat flux q°' and the purely spatial stresses P""^ can be 



obtained by direct computation [see (4.6a)-(4.6c)] 



p* = ^*n + e*s* - 



^, 



(5.62a) 
(5.62b) 
(5.62c) 



where h^ is the projection into the space orthogonal to U [see definition (4.4)] and where 
we have defined the quantities 



a'' = 8*10", 



Pa 



"Wa 



(5.63a) 
(5.63b) 



which are simply the spacelike components of the entropy flux and the thermal momentum 
respectively. 



It is, however, physically more interesting to note that the energy density (5.62a) could 



equally have been obtained as a Legendre-type transform on the master function in terms 



of the quantities (5.63). To see this, we also note that the definition of the generalised 



pressure, equation (5.52), is also a Legendre-type transform on the master function 



-^{p,0*,o) = A + pn + 6*s* -pa 



(5.64) 



where a = s* {waw"') ' and p = B^a. Thus, substituting (5.64) into (5.62a), the energy 
density becomes 

p*{n,s\p) = -K + paa''. (5.65) 
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This result is in complete analogy with (4.63). In fact, this relation informs us that the 



definitions (5.63) are indeed key variables in our thermodynamic model. Therefore, by 



virtue of (5.65) and (5.64), the variation of the energy density 



dp* = fi*dn + 9*ds* + a^dpa, 



(5.66) 



shows that indeed, the dynamical temperature (5.60b) is the thermodynamic temperature 
as measured in the matter frame U. 

Another point worth noting at this stage is that the variational analysis leads naturally 



to the presence of shear terms in the energy momentum tensor, (5.62c). Such terms are 



usually associated with viscous stresses and it is interesting to note that they arise even 
when we consider the pure heat conduction problem. 

To close our discussion on the thermodynamic interpretation of the model, let us note 
that in the limit when the two currents are parallel, cr"" = Pa = 0, the energy density (5.65) 
becomes p = —A and the starred notation becomes redundant. In this case, the quantities 



(5.60a) and (5.60b) become the usual chemical potential and thermostatic temperature 



respectively 



/^ 



dp _ 

dn 


dA 
dn 


2n ^\ = nB"", 
dn^ 


dp 
ds 


dA 
ds 


-.t^-'^' 



Thus, equation (5.62a) becomes the familiar equilibrium Euler relation 



p + ^ = pn + 6s, 



(5.67) 
(5.68) 

(5.69) 



which gives a consistent and unambiguous thermodynamic interpretation of the two-fluid 
model, at least in the case when the relative flow vanishes. Furthermore, this is an in- 
dication that the non-equilibrium case is effectively characterised by the misalignment of 
the two fluxes. However, a correct understanding of the nature of equilibrium can only be 
achieved by means of the second law of thermodynamics; through the entropy production. 



5.4 Rheometric covariant thermodynamics 

Now we are about to present the central result of this chapter. Its simplicity will be a direct 



consequence of the dynamical multi- fluid programme introduced in section 5.2, combined 
with the ability to draw thermodynamic conclusions from such approach. 

We have seen that - independently of any observer - the conservation of the particle 



number density current (5.7) implies that each of the forces (5.50) are [/-orthogonal tensors. 
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Also, the result from the previous section shows that the multi-fluid formalism has a natural 
thermodynamic interpretation when the physical measurements are taken in the matter 
frame U. 

Let us make a pertinent remark about the role of general covariance in a generic theory 



of dissipation. It is clear that (5.61) is a mere decomposition of the manifestly general 



covariant tensor (5.51) relative to a particular observer; the matter flow U in our casan 



For a manifestly general covariant theory of matter - defined as in section 2.1.2 - the 
internal variables, e.g. heat flow or shear stresses, play a passive role in the dynamics; 
they are implicit in the equations of motion. It is the specific choice of observer that gives 
them a definite physical interpretation; their meaning is bound to such a choice. 

It may seem that the above remark represents a backwards step in our pursuit a cov- 



ariant theory of thermodynamics. Quite the opposite, the symmetry of the action (5.30) 
is the complete group of diffeomorphisms of spacetime. However, we have shown that in 
order to obtain non-trivial dynamics in the physical case where one has a conservation 



law of the form (5.7), the relevant symmetry is that of the rheometric group Qrj. Thus, 



although the dynamics is clearly independent of any possible choice of observer [see (5.50), 



(5.53) and (5.54)], by construction, the most suitable family of observers from a physical 
and geometric point of view are those moving with the matter flow. In this sense, local 
deviations away from thermodynamic equilibrium have a natural description in the matter 
space ^{Q,U). 

The study of departures from equilibrium is objectively characterised by a non- vanishing 
entropy production. As discussed earlier in this chapter, we will use the explicit freedom 



present in the 'force' equation (5.50b) to explore the dynamical consequences of a positive 
entropy production term. In the rest of this section, we will divide our thermodynamic 
analysis in two parts. First we deal with the internal 'rheometric' dynamics, taking place 
in the matter space 3^. Then, in the second part, we shall discuss the implications of 
thermodynamic departures from equilibrium in the 'objective' covariant dynamics of the 
system. 



■'This remark trivialises a point often made by Garcia-Colin and collaborators [30] in their efforts to 
merge the principles of linear irreversible thermodynamics with those of relativistic fluids. They argue that 

...this [the inclusion of the heat flow in the stress-energy tensor] is against the tenets of the 
general theory of relativity, the stress energy tensor includes only all forms of mechanical 
energy... Heat cannot be incorporated into its structure... |30| 



From our discussion, it is clear that the presence of the transverse energy flow (5.62b I in a particular 
expression for the energy-momentum tensor is a consequence of the choice of observer. We could have 
equally chosen another observer whose four-velocity is aligned with a time-like eigenvector of the stress 



energy tensor at each spacetime point. In such case, the decomposition (5.611 would be free from terms 
containing g". However, such an observer will see (in a general case) a net particle flux in its rest frame. 
Historically, this is the difference between the Eckart (matter) frame and Landau descriptions. 
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5.4.1 Internal dynamics: the relativistic Cattaneo equation 

The aim of this subsection is to show that the internal thermodynamic evolution of a 
material system - represented by the particle number density current - leads in a natural 
wajjjto a relativistic analogue of the Cattaneo equation [see (2.62)]. 

Let us first note that by making use of the definition of the transverse energy flow 



(5.62b), we can re-express the entropy flux (5.57) as 



s*u'' + 



(5.70) 



The conjugate momentum (5.35b) can also be decomposed in terms of the matter four 



velocity U and the heat flow q. Although the calculation is straightforward, we wish to alert 
that such decomposition highlights a feature which is of prime importance to guarantee 
the stability and causality of the theory; entrainment. As a reminder, entrainment is a 
generic effect in multi-fluids systems which tilts the momenta relative to their respective 
currents, i.e. in a generally coupled multi-fluids system, currents and momenta are not 
'aligned'. One way of characterising such effect is given by the scalar quantity 



-A'^'ns*. 



(5.71) 



Thus, using once again the expression for the heat flow (5.62b), it is easy to see that the 



decomposition of the thermal momenta (5.59b) can be written as 



1 e^-e* 



(5.72) 



Finally, in order to preserve the clarity of the analysis, let us introduce the shorthand 
notation for the quantity inside the brackets in the above equation 



/3 



1 9^-1: 
s* s*6* 



s*e* ■ 



(5.73) 



It is worth noting that entrainment makes the decomposition (5.72) non-trivial. If we 



ignore the entrainment coupling between the currents in the Lagrangian density by setting 
A'^^ = 0, it follows from the definitions (5.60b) that 6^ = 6*. This particular case was 



considered by Carter [llj and, from the analysis made by Hiscock and Olson |57j . we know 
that it leads to a model that exhibits stability and causality problems. We will shortly 



^By now it should be clear that, when dealing with non-equilibrium thermodynamics, we understand 
'natural' to mean the simplicity motivated quadratic ansatz for the entropy production. It cannot be 
overemphasised that such ansatz is the simplest - not the most general - way to ensure that the second law 
holds. 
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return to this point. 

In terms of the decompositions for the entropy flux and thermal momenta, equations 



( 5.70 ) and ( 5.72 ) , and by virtue of the orthogonality of the force (5.50b ) with the normalised 



vector field U tangent to the particle world-lines, the entropy creation rate is found to 
satisfy 



a b 



,o.„b 



s%r = 2eia-b]u''s'' = 20[,.,]n"(?' 



1 



(5.74) 



Expanding the RHS of (5.74) immediately leads to 



s .„ 






(5.75) 



As we have been doing in the previous chapter, the simplest assumption one can make 



to secure the positivity of (5.74) is to introduce the constitutive equation for the heat flow 



q" = -Ke*h 



ab 



,* 1 . U^P . P ^ U^ 



which we can immediately re-arrange to obtain a Cattaneo-like relation 
Here we have introduced the effective thermal conductivity 



and the thermal relaxation time 



l + K/3 

kP 



(5.76) 



(5.77) 



(5.78) 



(5.79) 



l + Kp 

It is worth noting that, if /3 varies on a time scale Tp which is long compared with the 
relaxation time, then 

kR t 

(5.80) 



■ kP T ^ 

up ~ — ~ — *^ Ij 

Tp Ti3 



and in this case one would simply have 



K ~ K and r ~ np. 



(5.81) 



Thus we see that the dynamics obtained from the constrained variation of the two-fluid 
Lagrangian as described by an observer moving together with the matter flow, not only 
has a natural thermodynamic interpretation, but it provides us with the key constitutive 
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equation to obtain a hyperbolic theory of heat conduction. This was already the case of 
Carter's regular model, however, it is now clear that the omission of entrainment would 



make the relaxation time, equation (5.79) together with (5.73), independent of the matter 



description as given by the master function or through an equation of state. This disad- 
vantage was the main tool used by Hiscock and Olson [57J to show that, despite being 
hyperbolic, linear perturbations of realistic matter models (a non relativistic Boltzmann 
gas and a strongly degenerate Fermi gas) violate causality when described in terms of the 
regular model. Hence, the variational approach can be used to construct a hyperbolic 
theory of dissipation, but entrainment is the key feature to guarantee causality. The linear 
stability and causality analysis of a general two-fluid system will be developed in the next 
chapter. 



Thermal equilibrium and the Tolman-Ehrenfest effect 

So far we have been interested on the internal dynamics of our two-fluid construction. 



The central idea has been to relax the conservation law (|5.7|) for the entropy current and 

0, given by the 



analyse the variational dynamics governed by the Noether identity T"''' 



■,b 



pair of equations (5.53) and (5.54). We have seen in Chapter 2 that thermal equilibrium is 



defined to be the state for which the entropy production vanishes. Therefore, as we have 
previously argued, the misalignment of the two currents n" and s" corresponds to a generic 
non-equilibrium situation. 



In the present case, the only source of entropy is heat [c.f. equations (5.75) and (5.76)]. 
Therefore, thermal equilibrium is equivalent to the vanishing of the heat flow. In this case, 



the entropy is carried along with the matter and we obtain the usual Euler relation (5.69). 



By adding the force equations (5.50a) and (5.50b), we find that 



(n^ + s9) iia + hf'a ifJ-^bn + 9;bs) = 0. 



(5.82) 



It follows directly from the Euler relation (5.69) and the Gibbs relation for the pressure 



(analogous to (5.66) setting the last term to zero) that 

h\ [P.,b + {P + p)ub] = 



from which we immediately obtain 



iP + p)Ua = -h\P;b. 



(5.83) 



(5.84) 



As expected, we have the usual relation between the acceleration and the pressure gradient. 
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Finally, the entropy momentum equation can be cast in the form 

h\ [e-M + oiib] = 0. 



or, equivalently 



Oua 



-h\e,,. 



Comparing equations (5.84) and (5.86), we see that 

1 



1 



P + p 



h\P;b = SaO;h. 



(5.85) 



(5.86) 



(5.87) 



Thus, the equilibrium dynamics are characterised by a balance law between temperature 
and pressure gradients. It is worth noting that, even in the case of thermal equilibrium, 
the acceleration term plays a central role in the dynamics. 

Let us recall a result from section 3.4.1. In our previous discussion about thermal equi- 



librium [c.f. equations (3.54)] we conclude that in order to attain such states, a necessary 



condition is the existence of a timelike Killing vector field. In this case, the spacetime is 
said to be stationary. Therefore, in coordinates adapted to such a Killing field {gab,o = 0), 
the four-acceleration can be expressed in terms of the connection coefficients as 



-Too, with Too = -^goo,bg 



ab 



(5.88) 



Recalling that the acceleration is orthogonal to the matter flow, the entropy momentum 



equation (5.86) can be written as 



^,kg' 



1 



:9oo,k9 



(5.89) 



In the Newtonian limit, 500 is associated with the gravitational potential ^g by the expres- 
sion 



500 



1 + 



2$, 



(5.90) 



where we have restored the value of the speed of light to c. Substituting (5.90) back in 



(5.89) one immediately obtains the Tolman-Ehrenfest effect 



V^ 



where g is the Newtonian gravitational acceleration. 



(5.91) 
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5.4.2 Further remarks on the internal thermal dynamics 

Before moving forward in our analysis of the second law, it will be useful to consider a 
slightly different point of view. 

As a first step, we recall that the entropy force /^ is [/-orthogonal. Thus, contracting 



/° from equation (5.50a) with u"' and using the result in the force /^, equation (5.50b), we 
arrive at the rather elegant expression 



^*/a = -2'U''S {9[c0a]-b + (^b;[cGa]) ■ 



(5.92) 



This expression is manifestly [/-orthogonal. Moreover, it emphasizes the relevance of the 
entropy conjugate momentum Oa- However, if we want to gain insight into the key factors 



that contribute to the force, then we need to expand (5.92). To do this, we contract the 



entropy force (5.50b) with the entropy flux (5.70) to obtain 



{s^ea)T, = s'^fl = --q-fl. 



(5.93) 



This shows that the entropy production Fg = s"".^ only depends on the component of the 
entropy force that is parallel to the heat flux. In general, given that the entropy force has 
no components in the direction of the particle flow, we can make the decomposition 



fa = fha + t. 



(5.94) 



From (5.93), it is clear that f^ cannot contribute to the entropy production. Thus, there 



are two degrees of freedom in the entropy dynamics which are not constrained by the 
second law. This is an important point since there is no obvious way to distinguish the 
viability of models with different forms of f^ . 

It is straightforward to compute the values of /" and f^. The former can simply be 



read out from (5.93) as 



/" 



.3*9 



*2 






(5.95) 



This component should vanish when there is no heat transport. Again, there are a number 
of ways in which we could impose this condition. The simplest corresponds to our previous 



ansatz (5.76), leading to our generalization of Cattaneo equation (5.77). This is not a 



surprising result. However, it should be noticed that in this case we were never concerned 



about the second law at all! We arrived at this result solely by requiring that (5.95) 



vanishes in the absence of heat. As we have seen before, such an ansatz guarantees the 
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positivity of the entropy production Fs in (5.93). 



To compute the perpendicular component f^, let us introduce the projector orthogonal 
to q"" 

(5.96) 



a ra 

- b=(> 



1 



b - -^Q Qb- 



Then, using ansatz (5.76), it follows that 



J a 



-2s*e*h^„ V 



P;bq^ + ^q]b^ - qb;dq''(3 



(5.97) 



We see that the variational approach leads to the presence of terms that, even though 
they involve the heat flux, are not associated with entropy production. As far as we are 
aware, the dynamical role of these terms has not been discussed in detail in the literature 
even though similar terms are (as we will soon see) also included in the Israel & Stewart 
formalism. The variational model leads to these terms taking a specific form. In particular. 



those in equation (5.97) are all quadratic in q", the deviation from equilibrium. At this 



order, the most general case would allow a force of the form 



fit = h\ l.\ A,q^ + {q%,B + q^.^q'^C 



(5.98) 



with Aa, B and C unspecified coefficients. There may also, in principle, be first order terms. 



Clearly, equation (5.97) represents a particular case where all the coefficients follow from 



/3. Hence, the variational model is a particular example of the general class of permissible 
theories. The fact that all these models satisfy the second law of thermodynamics means 
that we cannot express a preference at this point. A very interesting question concerns 
whether there are situations where ft has a distinguishable effect on the dynamics of the 
system. 

Comments on the four-acceleration 

The key result of this section is the relativistic extension of Cattaneo equation. We obtained 
this by examining the dynamics of the two- fluid model in the Eckart frame, where the 
system has a natural thermodynamic interpretation and the state of the system is well 
described by the matter space. 

To facilitate a comparison with earlier results, let us recall the form of the relativ- 
istic version of Fourier's law as derived by Eckart and the constitutive equation for heat 
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transport obtained by Cattaneo [c.f. equation (4.31) and (2.62)] 



nh"'' [e,b + oiib] , 

d .: 



-K^—T 



dt- 



1,2,3. 



The structure of the variational constitutive relation, equation (5.77) combines the fea- 



tures present in both of the above equations, including the acceleration term which in the 
variational context arises as a consequence of the equations of motion. 

The role of the four acceleration can naturally be interpreted as being due to an effective 
'mass' per unit entropy given precisely by the temperature. As we will see later, this term 
has no counterpart in the Newtonian problem. Formally, this term originates from the 



local energy balance, equation (4.17). As we have discussed in the previous chapter [c.f. 



the discussion after equation (4.25)], it results from the fact that the infinitesimal 3-spaces 
orthogonal to the matter world lines are not parallel, but relatively tipped over because 
of the curvature of the world line. This leads to the interpretation of the four-acceleration 
contribution in terms of the effective inertia of heat. 



5.4.3 Equivalence with Israel & Stewart 

Since the two approaches are based on different strategies, any comparison between the 
variational model and the Israel &: Stewart results must be done carefully. Notably, the 
Israel & Stewart model is based on an expansion including terms up to second order in the 
deviation from equilibrium in the entropy flux. 

Meanwhile, the variational analysis did not involve such an expansion. As a result. 



the final equation for the heat fiux, equation (5.77), contains higher order terms. The 



corresponding Israel &: Stewart version [c.f. equation (4.50b)] reduces to 



q" = kTh 



ab 



T,b + Tiib + TpiQb +(^u") T\b 



(5.99) 



when we neglect the viscous terms r and t . We observe that equation is manifestly linear 
in q"". 

In order to compare the results, we focus on the linear deviation from equilibrium of 
the temperature 9* . Thus, we have 



r =T + 0{q^), 



(5.100) 
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and it follows that the generalised Cattaneo equation (5.77) can be approximated by 



(5.101) 



Here it is worth noting that in the second term in the first bracket, we could use the 
standard decomposition of Ua-b in terms of the expansion, shear, vorticity and acceleration. 
This would lead to terms that were explicitly excluded from the Israel & Stewart model 
at the point where we focussed on the case with r = r = 0. Basically, the variational 
analysis leads to the presence of terms that couple the heat flux to the shear and expansion 



of the flow. As these were artificially excluded from the analysis that led to equation ( 5.99 ) , 



we cannot count this as a difference between the two models. In fact, the full comparison 
carried out by Prion show that these terms agree in the two descriptions |65j . 

Keeping terms up to second order (treating the shear and the divergence of m° as 



first-order quantities), equation (5.99) can be written as 



rir + q" 



,afe 



kjh'''' {Tb + Tub) , 



(5.102) 



where the Israel &: Stewart thermal relaxation time is given by 



Ti = /3ikT 



1 + 



.kT' 



u 



(5.103) 



and the thermal conductivity is 



HI 



1 + 



.kT^ 



-i -1 



(5.104) 



Now it is clear that the two equations for the heat flux are formally identical, and 
we can identify the parameters in the two models. The upshot of this is that the models 
will only produce different results at higher order deviations from equilibrium. Given that 
this regime is hardly tested at all, we cannot at this stage comment on which of the two 
descriptions (if either) may be the most appropriate. Having said that, it is clear that the 
variational approach is formally elegant and the fact that it also applies far from equilibrium 
(at least in principle) may be relevant. An interesting question concerns whether there are 
situations where the, rather specific, set of higher order terms predicted by the variational 
analysis affect the non- linear dynamics. 
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5.4.4 Covariance and entropy production: an electrodynamic analogue 

To close this section, let us discuss an analogy with electrodynamics. Consider the motion 
of a congruence of test charged particles inside a region J7 of a fixed spacetime background 
under the influence of an electromagnetic field whose Faraday tensor is F"'^. The equations 
of motion for such a particle are given by 

w^-u^ = qF%u\ (5.105) 

where u"" represents the normalised timelike vector field tangent to the particles' worldlines 
and q represents its charge density. 

In classical electrodynamics, a source of the electromagnetic field acts effectively as a 
singularity, i.e. as missing point in the interior ofQ. Thus, for a source- free electromagnetic 
field, O, must be a simply connected region and the homogeneous Maxwell's equations, 

dF = 0, (5.106) 

imply the existence of a a gauge dependent potential one-form A such that 

F = dA. (5.107) 



In this case, we can write the equations of motion (5.105) as the 'force' equation 



/r = 2A[„;,]/, (5.108) 

where /|™ stands for the electromagnetic force and we have defined the charge density 
current by 

f = qu". (5.109) 

Thus, the electromagnetic force has the exact same form as that of the force equation 



of a multifluid system whose individual components are conserved [c.f. equations (5.50) 
with vanishing production terms]. In particular, for our two-fluid case, each of the con- 
jugate momenta //a and 6a, play the role of the potential one- form for electromagnetism. 
Moreover, has we have argued at the end of section 5.2.2, when more than one fluid species 
is present the dynamics allows for one of the species to have non- vanishing production, as 



in the case of (5.50b) 



This is the fact that we use to obtain an expression for the entropy production. There is 
a striking feature in this: What would be the analogous electromagnetic expression of the 
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entropy force (5.50b)? We can write down the answer straight away 

fr = 2A[,.,p' + j\Aa. (5.110) 

Note that, in this case, the gauge potential would have an observable effect in the motion of 



the test particles! Fortunately, we can invoke the variational principle from which (5.105) 
was derived and verify that charge conservation in i} implies that j"^ = 0, removing every 
gauge dependence from the dynamics. 

In electrodynamics, gauge invariance means that the potentials have no observable 
effects; only a potential difference is meaningful. In the same manner, in special relativity 
Lorentz invariance tells us that velocities on their own are meaningless, the physically 
relevant information is contained in relative velocities. Furthermore, in general relativity, 
general covariance removes the physical significance of individual spacetime points, only 
coincidences are important. Thus, let us close this chapter with the following question: 
What is happening here, that allows the thermal momenta to have an effect on the dynamics 
away from equilibrium':|j 



^Perhaps a principle of restricted covariance is forced upon us by the dynamics of thermal systems 
away from equilibrium. The answer to this question may give a physical meaning to the covariance of our 
dynamical laws. For a discussion on the physical significance of general covariance, the reader is referred 
to [171[77]. 
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6 



Adiabatic perturbations of equilibrium 



In the preceding chapter we developed a theory of heat conduction from a generic two-fluid 
variational principle. It would be natural then to expect that some generic features of a 
two-fluid system should appear in the context of our approach to relativistic dissipation. In 
particular, in this chapter we study the onset and development of instabilities arising from 
the interaction between the fluids, either through their chemical coupling or entrainment. 

We address the issues of stability and local causality of adiabatic perturbations, that 
is, we look for a criterion to be imposed on a given matter model (equation of state) in 
order to maintain the sound speeds bounded by that of light. Thus, this chapter presents 
an analysis of the local propagation of plane waves on a given background space-time. 

Compared to the standard single-fluid analysis, a multi- fluid system has more degrees 
of freedom; we need to account for relative flows between the different components. This is 
an essential requirement for the two-stream instability jSj- This is a generic phenomenon 
that does not require particular fine-tuning to be triggered, nor is limited to any specific 
physical system. The only requirements are that there must be a relative (background) 
fiow and a coupling between the fluids. Such instabilities are known to exist for a variety 
of configurations. For example, in shearing motion at an interface, we have the well known 
Kelvin-Helmholtz instability. This chapter presents a covariant generalization of the two- 
stream instability analysis to relativistic situations. 

The stability and causality analysis of relativistic theories of dissipation as presented 
by Hiscock and Lindblom [35] was tailored to suit the variables of the Israel & Stewart 
expansion. One of their main results, as we have previously discussed (c./ chapter 4), 
was that of the generic instabilities of 'first-order' theories. To close this chapter, we 
provide a different view regarding the stability of such theories. Building on the two-fiuid 
construction of the previous chapter, we carry out a 'consistent' first-order analysis and 
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show that one can have stable first-order models. However, when the relaxation time is 
varied over the dissipation time-scale, the system under consideration becomes unstable. 

6.1 Linear perturbations of thermal equilibrium 

In this section we deal with linear adiabatic perturbations of equilibrium states, i.e. those 
occurring without heat transfer. 

From the discussion in the previous chapter, the thermal dynamics of a single heat- 
conducting fiuid can be described in terms of the interacting components of an 'enlarged' 
two-fluid system consisting of particles and entropy. This view was motivated by the 
thermodynamic limit of non-equilibrium relativistic kinetic theory, where we obtained a 
well defined entropy fiux whose main effect is to redistribute the internal energy of the 
system towards the state of thermal equilibriurq^ In this manner, the internal dynamics 
of a dissipative fluid are given by a force-balance relation between a 'material' and a 
'thermal' fiuicQ Thus, the dynamics of thermal equilibrium is given by 

(6.1a) 
(6.1b) 
(6.1c) 



n% = 0, 




s% = 0, 




-'- a \b J a ' J a 


= 



Here T is the canonically defined two-fiuid energy momentum tensor [see equation (5.51 )], 



n"" and s"" are the conserved particle number density and entropy density currents and, using 



the conservation laws (6.1a) and (6.1b), the forces /" and /^ are given by [see equations 



(5.50)] 



/° = 2/i[,.,]n^ (6.2a) 

/a = 20Ms'- (6.2b) 

To obtain the linearised dynamics of the system, we consider only those perturbations 
which have no effect on the background geometry, that is, those for which 6gab = 0. Thus, 
considering the set of linear perturbations 

n-^n^ + Sn-, f-^f- + Sf^, (6.3a) 

s-^s'^ + Ss'^, f^^f^ + 6fl, (6.3b) 



^See the definition at the beginning of section 2.2. 

^See the discussion in section 5.2.2. In particular, the variational definition of the forces (|5.50[), and the 



force balance expressed by (5.531 and (5.541 
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the linearised equations of motion are 

5n% = 0, Sf^ = n''6fi^a-h], (6.4a) 

55% = 0, 6r, = s'69[,.^,], (6.4b) 

which must satisfy the energy conservation constraint expressed by 

Sfa + Sf: = 0. (6.5) 

Let us make a remark about the dynamics of thermal equilibrium. As discussed at the 
end of section 5.2.2, in thermal equilibrium the two matter-spaces used to define each of 
the currents coincide. That is, for a compact region of spacetime fi, the one-parameter 
families of diffeomorphisms generated by each of the particle species' worldlines, are the 
same. Therefore, there is a unique normalised timelike vector field whose integral curves, up 
to reparametrization, represents the two different fiowq^ In this case, where the dynamics 
are completely constrained, the system becomes completely integrable, i.e. 

fJ-la-M = d[a;b] = 0. (6.6) 



Thus, each of the forces (6.2) vanishes independently. Moreover, if the perturbations are 



adiabatic, we need to require, in addition of (6.5), that the perturbed components of the 



force balance (6.5) individually satisfy 



n''Sn[a;b] = 0, (6.7a) 

s'Se^^-M = 0. (6.7b) 

If this where not the case, the non-zero value of the perturbed force acting on each dis- 
placement of the individual currents would generate a relative flow between them, thus 
producing heat. 

In order to express the linearised equations of motion, let us recall the definition of the 



conjugate momenta ^a and Oa given in the previous chapter [equations (5.35a) and (5.35b)] 



where the coefficients B^, B^ and A^^ are given by the partial derivatives of the master 



^See the discussion in section 5.1. 
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function A(n,s,j) [equation (5.36)] 

B-^-2p,, B^^-2^ and A^ ^ -^,. 

Therefore, the perturbations of the momenta which leave the background metric unchanged 
take the form 



S0a = gabis^SB' + B'6s'' + n'^SA'"' + A'^'Sn^). 



(6.8a) 
(6.8b) 



where the variations of the coefficients (5.36) are given by 

d 



SB'' 



d_ 
d-n? 



-2—^6A, SB"" 



2 -^(5 A and (5^"" 

OS'' 



_d_ 



6A. 



(6.9) 



Following the work by Andersson and Comer [6j, we introduce the 'abbreviations' for the 
quantities 



B2b = B'^Qab - 2-^gac9hdrfn'^, 



Bib = B'9. 



dB' 



ah 



•^-^QacghdS S 
OS'' 



A. 



ah — 9ac9bd 



^[n^.^ + nV] + ^"-'^ 



I 9B^ , c J J dA^^ ^ J 



ATb = -9ac9bd (^^[^'^ + ^ "1 + Qj2 



■n n 



dB^ , dA^^ 
A'S, = A^'^Qab - gac9bd ( ^n'^n'' + -wT^s'^s'^ + 



jjjZ jjjA jjjZ 



c d 

s n 



/(.ab 



Xab 



JB^ c d 

-2^r-Y9ac9bdn s , 

OS'' 

^dB' c d 

-2—^gacgbdS n , 
on'' 



(6.10a) 
(6.10b) 
(6.10c) 

(6.10d) 

(6.10e) 

(6.10f) 
(6.10g) 



in order to separate the perturbed momenta in terms of the variations of the currents. 
Thus, we have 



5f,a = {B^ab + A2)Sn' + ixTb + A2)Ss\ 
dOa = {Bib + ATb)6s' + ixTb + A2)6n\ 



(6.11a) 
(6.11b) 



This is a very significant result. In the previous chapter, we defined entrainment as the 
misalignment between one of the fluxes and its conjugate momentum due to the presence of 
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the other fluid. Such an efi^ect is a consequence of the dependence of the master function A 
on the relative flux between the currents, characterised by the coefficient A'^^ above. Here, 



equations (6.11) teh us that, even in the case that the master function is independent of 
the relative flux, the perturbed momenta have non- vanishing components in the direction 
of both perturbed currents. 

To understand better this more complicated setting, we step back for a moment to 
analyse wave propagation in a perfect single fluid. 

6.2 Sound waves in a single fluid 

To set the stage for our forthcoming analysis, we consider flrst the propagation of sound 
waves on a single perfect fluid given by the conserved particle current n". In this section, 
we perform a local analysis of linear perturbations of the fluid propagating on a generic 
background geometry given by the spacetime metric gab- 

Let us make the following plane-wave ansatz for the particle number density current 

5n" = A"exp(-iA;feX^). (6.12) 

Here we assume that both, the amplitude A"' and the wave four- vector A;", have vanishing 
covariant derivatives 

A\ = and ka^b = 0. (6.13) 



For a single perfect fluid, the perturbed momentum (6.11a) reduces to 

5^la = Bl^5n\ (6.14) 

Following the convention set in the previous chapter, we work in the Eckart frame. Thus, 
we can decompose the amplitude and the wave vector into their co-moving and transverse 
components 

k'i^ = kr,{anu"' + kl), (6.15a) 

Al = A'lu'' + Al^. (6.15b) 

Here, we have made the usual definitions 

^nCn = —kau"" and kjz^ = h ^kb- (6.16) 

where fc" is a unit three-vector pointing in the spatial direction of the wave propagation. 
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Similarly, the components of the amplitude are expressed by 

and (7n measures the phase velocity of the waves as seen in the Eckart frame. 



(6.17) 



The perturbed conservation law, equation (6.4a), becomes 



dn% = [A'' exp{-ihx%a 

= [A% + iA^'i-h.^ax' - 6'M] expi-ikbx') 
= iA°-ka exp{-ikdx'^) = 

which implies that the waves are transverse in the space-time sense 

kaA"^ = 0. 



(6.18) 



(6.19) 



It is a straightforward exercise to show that the contraction A'^B^, ,, is identically zero. 



c[a;b] 



Therefore, (6.7a) becomes 



= in'A''k[i,Bl.exp{-ikdx'') = 0, 



and thus we express the linearised force by 



n%bB:.A'' = 0. 



(6.20) 



(6.21) 



The dispersion relation can be easily obtained by contracting equation (6.21) with k . 
Assuming that ricA'^ ^ 0, we may define the speed of sound by 



2 Slogs'^ 

c = 1 H — . 

^ d log n 

Thus, we can rewrite the coefficient B^^^ in terms of the sound speed c^ as 

dB^ 
B^f, = B'^gab - 2-^gacgbdn''n'^ 

2n2 dB"" 

gab + UaUb - UaUb - -T^TT^UaUb 



(6.22) 



B"{hab-c\uaUb). 



hab - UaUb 1 + 



B"^ 9n2 
dlogB""^ 
dlogn 



(6.23) 
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Substituting this result in the first term of the anti-symmetrised bracket of the reduced 



perturbation of the force, equation (6.21), we obtain 

k^n^kaBlAl = {kak'') v'B^ {he - clubUc) Al, (6.24) 

whilst the second term is 

k''n%B^,A'^ = (n^kb) k'^ [S" {hac - cluaU,)] A''. (6.25) 

Therefore, the single fluid dispersion relation takes the simple form 

k-n'k^^Bl^^Al = 0. (6.26) 

Note that this is equivalent to 

al-cl = 0, (6.27) 



where, from the first equation in (6.16), we have 



(Tn = -TrKUa. (6.28) 

Now we can see that Gn corresponds to the phase velocity of the waves as measured in 
the Eckart frame. Now, we just need to show that Cn is effectively the speed of sound. To 
see that this is indeed the case, it will suffice to introduce the pressure and recall that, for 
the single fluid case, the energy density is related to the master function simply by 

p=-A. (6.29) 

In this case, the energy density is a function of the particle number density alone. Thus, 
the Gibbs relation is just 

dp = pdn. (6.30) 

The pressure ^ is defined by the standard Euler relation p + ^ = pn, and in this very 
simple case we can write it as 

^ = -n + n-^ = -p + eV. (6.31) 

an 
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Therefore, 

dp dp dp 

= 1 -I- ri^ 

dp fi" dn 

= 1 + -TT^' 6.32 

dlogn 

and we can see that, indeed, c^ is the speed of sound. 

Let us recapitulate what we have achieved with this exercise. We studied the perturbed 
equations of motion of a single perfect fluid through a plane wave ansatz. We found the 
dispersion relation and defined the speed of sound in terms of variations of the master 
function for the matter model of our choice. This allows us to rule out non-physical 
models, for instance, if the perturbations end up propagating faster than the speed of 
light. In this sense, the relativistic analysis constrains our choices for master functions. 



and hence, of equations of state. Thus, equation (6.32) serves to analyse when the matter 
model we use will give non-physical results. Therefore, let us discuss a bit further the issue 
of local causality. 

On the local causality of the selected matter model 

Consider the wave vector k'^ of the plane wave ansatz. From the above dispersion relation 
we see that 

kak^ = kl{l-cl). (6.33) 

Thus, for causal wave propagation, the speed of sound should be less than the speed of 
light, in natural units 

cl < 1. (6.34) 

Note that the amplitude A'^ of the perturbed current defines the character (timelike, space- 
like or null) of the perturbation. In order to have causal disturbances, A"" must be a timelike 



vector field. Thus, from the transversality of the waves, equation (6.19), k"' is a spacelike 
vector. 



Let us also note that, when the force equation (6.11b) is evaluated in terms of the 
solution to the dispersion relation, we find that, 

Al^ = anAp^. (6.35) 

Thus, the waves are transverse in the ordinary three dimensional sense. 



Equation (6.34) is the first encounter we have with a restriction imposed by causality. 

113 



6.3. Sound waves in the matter and entropy system Adiabatic perturbations 

This implies a constraint on the master function 

— ^ + 1 < 1. (6.36) 

aiogn 

In addition, in order to avoid absolute instabilities, the speed of sound should also be 
bound from below, i.e. 

c\ > 0. (6.37) 

This restriction prevents cases in which cr^ < 0. This can be written in terms of the master 

function by the inequality 

d loe B'' 

^fP^ + 1 > 0. (6.38) 

alogn 

Combining these two results, we obtain a criterion that must be satisfied by any proposed 
matter model that guarantees absolute stability and local causality of the perturbed system 
in terms of the defining master function alone. That is 

d loe B"" 
- 1 < ', < 0. (6.39) 

alogn 

In the following section, we will extend this analysis to obtain a stability and causality 
criterion for the two-fluid system we are concerned with. 

6.3 Sound waves in the matter and entropy system 

Having discussed plane wave propagation for a perfect fluid, we use the same technique 
for the conducting model developed in chapter 5. The aim of this section is to work out 
the dispersion relation for wave propagation in the particle-entropy two-fluid system. To 
make the discussion more general, we relax the thermal equilibrium assumption and work 
out the perturbations in the case where there exists a relative flow between the fluids. 
In this sense, the particular case of thermal equilibrium (aligned flows), will be discussed 
separately in section 6.3.2. 



Using the plane wave ansatz [equation (6.12)], we write the perturbations of the two 
conserved currents as 

dn" = A^ exp{-ikbx'') , (6.40a) 

ds" = A^ exp{-ikbx^). (6.40b) 

Let us note that in these expressions we do not attach a constituent index to the wave 
vector k'^, since waves in a system are such that the constant wave vector is the same for 



all fluids. Thus, from the component decomposition of the wave vector, equation (6.15a), 
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we have 



r = k^ianu'' + K) = Ua.ut + K)- 



(6.41) 



Also, since the dispersion relation is a scalar equation, it will be useful to express the inner 
product Ti;"fe° (where w"" = w°'/w) in terms of the angle /3ns between the two vectors 



w"'k^ = cos /3n 



(6.42) 



Note that there is significant difference from the previous single fluid case. For a 
generic multifluid system, we have to take into account the relative flow w"' between the 



fluids introduced in section 5.3 There, we expressed the entropy flux in the Eckart frame 



by the relation (5.57), that is. 



SMg = Sj{u°- +11)"). 



(6.43) 



Having introduced the entropy current raises a subtlety. The quantities a^ and k^ are 
what would be measured in the matter frame. We could choose the frame defined by 
the entropy flux instead, and make our measurements with respect to the entropy flow. 
Although, from the mathematical point of view there are well-defined transformations 
between the two descriptions, working in the entropy frame would produce unnecessary 
complications due to the fact that thermometers and accelerometers are naturally carried 
by the particles in the fluid; not by the phonons of the system. Moreover, the relative flow 
of the particles with respect to the entropy frame, which we may call w^, will clearly have 
the same magnitude as tu" and is just 



w^ 



■j{w u"- + w"-). 



(6.44) 



Denoting the normalised four velocity associated with the entropy flux by u" , it is useful 
to note that we can write the two four- velocities in terms of their relative flow as 






(6.45) 
(6.46) 



Thus, by contracting each of the four- velocities with the wave vector ka, we obtain a matrix 
equation of the form 



wan - cos /3ns "7 ^ COS /3sn 
-7"^ COS ^ns was - COS /3sn 



fcn 








(6.47) 
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The determinant of the 2x2 matrix should vanish, and this leads to 

crn-WCOS/3ns .„ .„s 

as = cosPsn ^- (6.48) 

wan - COS Pns 

Thus, we have shown that if cj^ < 1 then o"g < 1. This is a natural conclusion since 
causality is a frame independent requirement. 

It was shown in the Newtonian case [5] , that two chemically coupled fluids may become 
two-stream unstable. This result was later generalised in [67] for the relativistic case. In 
what follows, we reproduce the same analysis. 

Let us begin by noting that the perturbed conservation laws, equations (6.4a) and 
(6.4b), give the same transversality condition as in (6.19), that is. 



kaA^ = kaA'^ = 0. (6.49) 

However, the perturbed forces are more complicated than in the previous case. The cov- 



ariant derivative of the perturbed momentum (6.11a) is 



6^la■,c = 5n\^,{Bl, + All) + Ss',,{x:i + ^aD 

= ikM'nKb + AiAlfA^xll + A'^All) exp(-ifc,x'^). (6.50) 

Thus, we encounter our first problem. The perturbed equations of motion satisfy 
a similar force balance as in the previous chapter, but that does not directly imply the 
vanishing of the individual perturbed forces. In a generic multi-fluid system, the discussion 
becomes significantly more elaborate at a very early stage. In our case, however, we can use 
the fact that, in the strict equilibrium sense, we can use the integrability of the individual 
forces to obtain 

= JC:,Ai + JC:iA', = 0, (6.51) 

where we have introduced the 'n-dispersion' tensors 

/C°, = n^(fe[,e^], + fe[,^^P, (6.52) 

jc:i = n'^ik[,x:!t+k[cA:^b)- (6.53) 

Similarly, the perturbed entropy force is 

6f, = ICl,Al + IC2A'n = (6.54) 
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where the 's-dispersion' tensors are 

)Cl, = s'^{k[,Bl^, + k[,A':^,) (6.55) 

^S = «^(A:[cX:^ + M:]^)- (6-56) 



FoHowing the argument in section 6.2 we contract both equations (6.51) and (6.54) 
with the wave vector fc". This wih resuh in a 2 x 2 matrix problem which we address 
by systematically increasing the degree of interaction between the two fluids. In order 
to do so, it is convenient to isolate further the different contributions that appear in the 
dispersion matrices. 



^^6 = n=fc[,e5fe = -^BMKcTnhab + clKub), (6.57) 

'I, = s'k^,B% = -\B^s{k^a,hl, + c^.kaul), (6.58) 



The bulk contribution in the /C^^^ and /C^^ matrices can be reduced to ^67J 

1 

2 
1 

2 
while the entrainment piece becomes 

s^^l^ = rfk^cA'^f, =-7n|B|JjJ(A;n<TnU;a - 2A;a)nf, - kaWh] 

+ '^-A'^s{kn(TnWa - ka){ub + Wfe) k (6.59) 

n ' ) 

<h = S'^[c^a]6 =lfs{^ln[{ks<ysWl - 2ka)ul - kaWl] 

+ 7-^Tn(fesas< - ka){ul + wl)\. (6.60) 

s ' ) 

We have also introduced here the orthogonal projector /i^^ associated with the entropy 
four-velocity u^ and the speed of sound for the entropy perturbations Cg [c.f. equation 



(6.22)] 



Rn = 

'-',ns - 

RS , 

A-as ^ 
•^,ns - 


- 1 / 


' dB^ 
' dB' 

3 

~ 5 


- 1 / 


dA'' 



2 51og^" 

Cg = 1 + ^r^ • (6.61) 

''Slogs ^ ' 

We have introduced the shorthand notation for derivatives of the coefficients B^, B^ and 

(6.62) 
(6.63) 
(6.64) 
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and for the cross-coupling terms in IC^^ 



fXy ., 



ab 



C^n 



ab 



n^cX'a 



S '^[cXa]b 



■Cjn^/B^ka{ub + Wb) 
C-fsy/WB^ka{ul + wl). 



(6.65) 
(6.66) 



Here, 



e 



' ^2„/^° 



1 / dB^ 



Finally, we see that the entrainment part has significantly more presence 



(6.67) 



^ab — ^ l^lc-^a 



{ 



b 
A^'^'knanhab 



n ' \ n ' 



+ ik^an (^T-'Bfns + -^rns) ^aUbj, 



c 



^ '^lc-^a]b 
S 



\{ 



A'^'ha.hl 



■^si^sii-ab 



A''' + 



n 






n. 



+ ^haJj-Bl, + A: 



,ns 



«}■ 



n 



+ ^[^-Bl^ + A 



kaUb 



kaUb 



(6.68) 



(6.69) 



As stated in [67], a multi-fluid system must have non-zero bulk properties, the other 
terms can be absent depending on the equation of state. We now proceed with our sound 
propagation analysis by gradually increasing the coupling between the matter and entropy 
fluids. 



6.3.1 The completely decoupled case 

In this very simple case, the two fluids are completely independent of one another. The 
forces vanish independently and the system is absolutely integrable. Thus, although the 
currents may not be aligned, the dynamics are not affected by any relative flow. In this 
particular case, the only non-zero contributions to the dispersion matrices /C are equations 



(6.57) and (6.58) and we can proceed as in section 6.2. Thus, the contraction of equation 



(6.51) and (6.54) with the wave vector 



k'^df^ = 0, 

k-5r, = 0, 



(6.70) 
(6.71) 
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produces the simple 2x2 matrix problem 



Q Rs^^2 „2^ 



S°(ct- - C„ 






Therefore, the dispersion relation is simply given by 

{'yl-cl){al-cl)=Q 



(6.72) 



(6.73) 



As we have seen before, a^ corresponds to the squared phase three- velocity as measured 
in the Eckart frame. Similarly, a^ is the squared phase velocity of the 'entropy waves' as 
measured in the frame defined by li" , whose transformation to the Eckart frame is given 



by (6.48). Thus we see that, for the uncoupled system, we must require that the equation 



of state satisfies (6.39) for both, particles and entropy. 



This simple case has served the purpose of defining the quantities Cn and Cg. As we 
have argued in the single fluid case, these quantities are identified with the decoupled sound 
speeds. It will be useful to keep this in mind in the following more elaborate discussion. 

6.3.2 Cross-constituent coupling 

The next situation corresponds to that where the system is sensitive to the relative flow 



w"" but not to entrainment. That is, in addition of (6.57) and (6.58), we have contributions 



due to the cross-coupling tensors (6.65) and (6.66). Again, we reduce the problem to a 



2x2 matrix by contracting (6.51) and (6.54) with the wave vector k"". The cross-coupling 



produces the expected off-diagonal terms 



-VWWC B\al-cl) 



ulA'i 



(6.74) 



where C^ is the cross constituent coupling defined by 



C 



^n^s 



2ns 



as° 



ds 



(6.75) 



The dispersion relation is again the determinant of the 2x2 matrix and in this case 
gives 

[al-cl){al-cl)=C\ (6.76) 



We can see from the definition of C , (6.75), that in order for it to be less than zero. 



either B^ or B^ should be negative. For ordinary matter, or entropy, this is not generally 



119 



6.3. Sound waves in the matter and entropy system 



Adiabatic perturbations 



the case. Hence we do not deal with such possibihties in this work. 



Equation (6.76) shows the richer phenomenology of the coupling between matter and 



entropy, or of a two-fluid system in general. This allows for two-stream unstable flows, 
as discussed below (see also [67j and [5]). Once more we are going to take small steps to 
analyse this situation. 



Aligned flow: the stability of thermal equilibrium 



We can now use the angular relationship between o"n and as given by (6.48) to write a 



dispersion relations solely in terms of a^- Addressing the question of absolute stability 
(Cn > 0), we set the relative flow w"' = in (6.48), for which we would have as = Un and 



(6.76) becomes the quadratic expression in a'^ 

I 2n2 / 2 I 2\ 2 I / 2 2 n'i\ 
K) -(Cn + Cs)<7„ + (CnCs "C ) 

whose roots are 



0, 



(6.77) 



cr„ 






r.2^2 



4C2 



(6.78) 



In order to avoid complex u^, the square root in (6.78) has to be real, which follows 
from the fact that C^ > 0. Now, to have absolute stability it is sufficient to require that 



the absolute value of the negative square root in (6.78) be less than the first term, that is 



{ci^ciY>{ci-ci)^^C\ 



that is equivalent to the requirement 



C < Cj^Cg . 



(6.79) 



(6.80) 



This can be written in terms of the bulk quantities we expect to be positive, B'^B^ > 0, 
and so absolute stability constrains the master function to satisfy 



dlogB'' dlogB^ / dlogB"" 



d log s d log n 



dlogn 



1 + 



d log B^ 
dlogs 



(6.81) 



The causality constraint requires a^ < 1. For this we notice that if the absolute stability 
condition (6.80) is satisfied, and C^ > 0, we need only to guarantee that the positive root 
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of (6.78) is causal, since the negative one is always smaller 



^(c^ + a + \/(cS-c2)2+4C2<l 



V(C2 - C2)2 + 4C2 < (1 - d) + (1 - c!) 

C'<{l-cl){l-cl), (6.82) 

which in terms of the master function becomes 

d log s d log n ~ d log n d log s 

Note that causality and absolute stability set upper bounds for C^. In general it can 
be shown [67] that if c^ + Cg < 1 then any absolutely stable master function is also causal. 
Conversely, if c^ + Cg > 1, any causal master function is absolutely stable. 



Equation (6.83) is one of our most important results since it represents the stability of 



thermal equilibrium. Thus, a matter model whose master function (equation of state) fails 



to satisfy the inequality (6.83) can be regarded as non-physical. 

In order to obtain a comparison with the Newtonian problem, we consider the case of 
'small' relative flow. Although this would not strictly correspond to an adiabatic perturba- 
tion, it will allow us to show the similarities with a better known problem. In this case, we 
will be in position to examine if a dynamical two-stream instability is present by solving 



(6.76) for some relative flow and then take the the limit in which the flows are aligned. In 



what follows we will assume that the master function is absolutely stable and causal. 
Small relative flow: 'almost' adiabatic perturbations 



Assuming that w and dn are both much smaller than the speed of light (c = 1) in (6.48), 



the dispersion relation (6.76) becomes 

(f^n - Cn)[(^n - W COS p^^f - C^] = C^. (6.84) 

Introducing the variables 



x = ^, y= ^^°^^- , b'=('-)\ a' = % (6.85) 

Cs Cs \Cs J C„ 



we get 



^^[(x-y2)_i] = i. (6.86) 
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Written in this form, the problem is identical to the Newtonian plane wave propagation 
discussed by Andersson, Comer and Prix [5] where they show the regime in which the two- 
stream instability operates. It is also shown in |67) that the particular example used by 
Andersson, Comer and Prix obeys the causality and absolute stability criteria derived 
above. First, let us note that, due to the presence of a velocity scale given by the speed 
of light, we cannot completely scale out the velocities. Thus, the relativistic analysis of 
stability will, in general, contain an extra parameter compared to the Newtonian case. 
Following |67j, we take that parameter to be the speed of sound for the entropy, Cg, which 
without loss of generality, we can take larger than the Cn- Using these parameter the 



absolute stability criterion (6.80) becomes 
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<1, 



(6.87) 



which is satisfied in the model discussed above. 



We now turn to the relativistic dispersion relation (6.78). Written as an equation for 



(Tn, it constitutes a non-trivial quartic. If we use the same re-scaling (6.85), then (6.78) 

becomes 

(x-y)2 



7-2(1 



c22;2 



) + c^{x-yy 



1 



1. 



(6.88) 



Although the general solution of ( |6.88 ) is available, it is quite a complicated expression 
and offer very little insight for our analysis. Thus, we tackle the problem numerically for a 



set of parameter values that maintain absolute stability and causality. Figure 6.1 provides 
plots of the real and imaginary parts for the solutions of (6.88) in the case when the 



particle and entropy flow are aligned. The solutions for a^ are taken to be functions of the 
relative flow parameter y and the coupling C'^, with the choice Cg = 0.5, b"^ = 1 and = 0. 
Non-zero values for Im((Tn) in the figures indicates the presence of an unstable mode. The 
appearance of unstable modes is reflected in the real parts whenever two frequencies merge. 
This behaviour is typical for this kind of dynamical instability (SjEZj. Note that, although 
in the central figure the instability sets at a 'high' relative flow {w ~ 0.6), unstable modes 
can be produced for arbitrary small values of the coupling C, as can be seen from the third 
column. 

6.3.3 The role of entrainment 

We finally turn our attention to the full problem. Up to this point we have introduced 
three equation of state parameters, i.e. Cn, c^ and C, that are obtained as second derivatives 
of the master function A. Given that we have established the two-stream instability for 
the general cross-constituent coupling and arbitrary background flow, the main reason for 
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Figure 6.1: Plots of the real (top) and imaginary (bottom) parts of the mode frequencies 
(Tn as functions of w^ for Cg = 0.5, 6^ = 1.0, and /3ns = 0. The merger of two frequencies, 
and subsequent non-zero imaginary parts, signal the presence of a two-stream instability. 
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discussing entrainment is because it is an essential ingredient for the correct interpretation 
of relaxation times in the Cattaneo-like equations derived in Chapter 5. We will simplify 
the entrainment case by assuming that the system is close to equilibrium, that is the 
relative velocity is such that w <^ 1, and that the flows are aligned. This does not mean 
that the sound wave speeds or phase velocities have to be small, nor do the flows n" nor 
Ug have to be similarly restricted. 

If we keep deviations from thermal equilibrium to 0{w'^), the master function can be 
approximated as [U |3] 



A = Ao(n2, s^) + Ai(n2, s^){f - V^^), 



(6.89) 



which immediately implies 

^'i, = -2ns 



dXo 9Ai^ ,2 






and 



-2ns 



A" 



dn? dn^ 
dXi 

dXi 



(j2 _ ,/^^) _ Ai 



2n 



as2 



-Ai 



•^,ns ^- 



(6.90) 

(6.91) 
(6.92) 

(6.93) 



In the spirit of understanding a complicated problem, we shall make a further simpli- 
fying approximation by assuming that Ai is constant, so that 






0. 



(6.94) 



This leaves the problem with the simple entrainment parameter Ai. With this approx- 
imation we find 



B^{al - cl) 



-CVbW - A'^'^'^^ial - 1) 
-CVB^ - A-^^{a^ - 1) B^{a^ -cj) 



It follows from (6.33) and (6.41) that 



ulA- 





(6'.95) 









(6.96) 
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and hence the dispersion relations becomes 



{<^l - cl){al - cl) -\C + C^^{l-al){l-<yl) 



(6.97) 



where 



C« 



^" 



^/WW 



and because the inverse of (6.95) must exist, \C\ / 1. 



(6.98) 



We stress that, unhke the simpler cases, the dispersion relation (6.97) does not hold for 



arbitrary propagation direction with respect to the relative velocity. This mirrors the fact 
that entrainment enters the master function in a fundamentally different way: At first- 
order in the relative velocity squared. In the dispersion relation, however, entrainment 
contributes even in the limit of zero relative velocity, as we have seen in the previous 
section, because there are still two sets of interacting sound waves. 

In the limit of aligned flows and small deviations from thermal equilibrium [w <^ 1), 



we see from (6.48) that a^ = a^, and so 



[cr„ 



cDi'jl-c!)- C + C\l-al 



0. 



(6.99) 



Even though when it is difficult to think of a system of fluids which are entrained and 
yet not chemically coupled, it is particularly instructive to consider the entrainment case 
alone, i.e. set C = 0. The corresponding dispersion relation is a quadratic for a'^, and has 
solutions 



C^ + C^ 



2C^^± 



{cl-c!)' + 4C^'{l-cl){l 



(6.100) 



2(1 -C^2) 

For speeds of sound satisfying the causality requirement, Cn, Cs < 1 the discriminant is 
obviously positive and hence the a^ are real. In order to analyse absolute stability and 
causality we need to consider the ranges < C^ < 1 and 1 < C^ separately. We will look 
at C^'^ > 1 first. 



We rewrite (6.100) so that the denominator is positive 



2C^^ 



ci± 



{cl 



r.2\2 



+ AC^\l-cl){l-cl 



2{&' 



1) 



(6.101) 



Since we are imposing Cn, Cs < 1 and C' > 1, the terms outside the square root in the 
numerator are positive, therefore the positive solution is absolutely stable. But, we can 
also show that it cannot be causal. As for the minus solution, it can be shown that it is 
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absolute stable only if C^ < 1 [67J, which is not the case we are considering at the moment. 
Therefore, C^ > 1 does not lead to both absolutely stability and causality, and is therefore 
ruled out. 

In a manner similar to the cross coupling case, the region < C'' < 1 is a different 



story, because the terms outside the radical in the numerator of (6.100) are not of any 



definite sign. This affects the absolute stability analysis more than the determination of 
causality. In fact the causality requirement is sufficiently straightforward in this range for 

In order to assess absolute stability, it is useful to introduce [67] 



(4-ctr + 4C^^{l-4){l-ct) 



(6.102) 



This allows the numerator (6.100) to be written in such a way that the absolute stability 
condition becomes 

t2 - [±2(1 - cl){l - c2)] + [clil - c2) + c2(l - c^)] [2 - (c^ + C^)] < 0, (6.103) 



where the it corresponds to that of (6.100). The final step is to factorise (6.103) to obtain 

{r ± [2 - {cl + c^)] } {r ± [cl{l - c^ + 0^(1 - d))] } < 0, (6.104) 

where if the '+' is taken from the first factor, then the '— ' must be taken in the second, 
and vice versa. In either case, the factor that has the + is positive definite, and so the 
other factor must be less than zero. When the first factor takes the — the inequality leads 
to C^ > 1. The other choice leads to the more restrictive condition C'^ < c^c^- 

To summarise, we have shown that when C^ > 1, there is either no absolute stability 
or causality, which makes this range non-physical. Meanwhile, for < C^ < c^^c^ this 
system is causal and absolutely stable. In terms of our previous definitions, this translates 
to 



(A 



ns\2 



< B'^B'' 1 + 



dlogB'' 



1 + 



d log B^ 



(6.105) 



dlogn J \ dlogs 

as a constraint on the master function, when the relative speed between the two fluids 
is sufficiently small. A numerical analysis shows a similar behaviour to the one shown in 
Figure |6.2[ 



What have we learned? 

This chapter built on the relativistic two-stream instability analysis developed in [5l 167] . 
The central results of this chapter can be summarised by the following two points: 
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1. Although it was not the aim of this chapter to present a two-stream instabihty 
analysis for a system of 'generic' fluids, let us note that there is still work to be done 
in this direction. This is important because, although we have made some progress 
in our understanding of the role of entrainment in linear perturbations, we have not 
been able to produce a complete analysis of the fully coupled two- fluid system. Such 
analysis should become relevant in answering specific questions regarding the onset 
of instability for particular matter models, where the two-fluids under consideration 
are strongly coupled. 

2. In our present thermodynamical discussion we note that in all three cases - decoupled, 
chemically coupled and entrained - equilibrium is always stable against adiabatic 
perturbations. This is, perhaps, not surprising. However, notice that for some matter 
models, if the couplings are strong enough, there might be unstable modes in an 
arbitrarily small neighbourhoods of equilibrium, i.e. in the regime of almost adiabatic 
perturbations. In this sense, the emergence of a two-stream instability for the matter- 
entropy system is a concrete prediction of our multifluid programme for relativistic 
dissipation. 
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Figure 6.2: Plots of the real and imaginary parts of the mode frequencies a^ as functions 
of y and C^ , for Cs = 0.5, and 6^ = 1.0. The lines of intersection and subsequent merger 
of two frequencies signal the presence of a two-stream instability. 



128 



Non-adiabatic perturbations of equilibrium 



In the previous chapter, we have found the criteria that matter models should satisfy 
in order to have a stable state of thermal equilibrium. However, thus far we have only 
considered adiabatic perturbations, i.e. those which do not produce heat. Here, we will 
take a first step towards testing the stability of thermal equilibrium against deviations 
'down hiir the summit of maximum entropy; those where heat has a role to play. 

We have shown in Chapter 5 that the heat conduction problem has two degrees of 
freedom. However, we have only focused on the thermal one, leaving aside the one related 
to the material degree of freedom. Taking both into account, let us consider the case 
where the original state is that of thermal equilibrium, so that the two conservation laws 



for particles and entropy, equations (6.1a) and (6.1b), hold. As discussed in the previous 
chapter, by adiabatic perturbations we mean those which produce heat. Here, we consider 
the case of 'small' non-adiabatic deviations from thermal equilibrium, that is, those which 
are linear in the heat flux. Thus, we linearise the momentum equations associated with 
the transport phenomena given by the Cattaneo-like relations [c.f. section 5.4] to obtain 

aqa + aq" + h^af^;b + fJ'Ua = 0, (7.1a) 

Tqa + qa + K-h^aT;b + I^TUa = 0, (7.1b) 

where the coefficient a is directly related to the entrainment and the equilibrium temper- 
ature T by 

-^''' K - rs _ ^, 

1 Kn 



and we have used the definitions (5.73) and (5.79) to eliminate A'^^ in favour of the material 
properties k and r. 

Note that here we use the local equilibrium hypothesis since s* and 9* differ from their 
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equilibrium values only at second order. Moreover, the first order pressure is obtained from 
the standard equilibrium Gibbs relation 

P;a = A*;an + T;aS. (7.3) 

Therefore, in the case of linear deviations from thermal equilibrium, the fundamental Euler 
relation involves only equilibrium variables, that is 

p + p = fin + Ts. (7-4) 

The key point to note is that the first-order model remains cognizant of its higher order 
origins. Specifically, a and, therefore, r depend on the entrainment of the two-fluid system. 
Thus, although the deviations from equilibrium are linear, the dependence of the master 
function on the relative flux between particles and entropy is a second order effecljj which 
is a crucial ingredient to address the problem of non-adiabatic perturbations. 

In order to analyse the linearised dynamics of the heat problem, let us consider perturb- 
ations away from a uniform background state. We denote such perturbations by 5. Since 
we start at equilibrium, we have g" = n'* = d = 0. Therefore, the perturbed momentum 
equations can be written as 

a6qa + h''J {p,.},) + jidiia = 0, (7.5a) 

T5qa + 6qa + Kh\6 [Tb] + KT6ua = 0. (7.5b) 

Now, we can divide our analysis into two parts. First, we will study the transverse mode 
to obtain the requirements that need to be fulfilled in order to have a stable equilibrium 
state against non-adiabatic perturbations. Secondly, we will look at the dispersion relation 
for the longitudinal mode, containing further information about the stability of the theory 
and where we can impose local causality as described in the previous sections. 

7.1 Transverse modes: first stability criterion 



Noting that the perturbed momenta (7.5) only have spatial components with respect to the 
Eckart frame, we can take the spatial curl in coordinates associated with the background 
flow. Throughout this section, we will denote spatial indices in the traditional manner by 
the Latin letters i, j and k (not to be confused with spacetime indices a, b, c, d ...). Thus, 

^Sgo for instance the original discussion by Carter presented in the last section of chapter 4. 
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acting with the curl operator e'-' Vfc on the perturbed momentum equations we obtain 

^iil' + aQ' = 0, (7.6a) 

tQ' + Q'' + kTW = 0, (7.6b) 

where we have defined the spatial three vectors f7* and Q* by 

U' = e'^''VjUk and Q' = e^^^V jQk, (7.7) 

and e^^^ is the totally-antisymmetric three-volume form whilst Vj is the spatial gradient 
operator. 

Assuming that the perturbations have a dependence on time of the form e*"^* , we arrive 
at the dispersion relation for the transverse perturbations 

ioj [/i (1 + iojr) — iojanT] = 0. (7.8) 

Obviously, the case lj = is a solution. With the aid of the fundamental Euler relation, 



equation (7.4), it is not difficult to show that the other solution of (7.8) can be written as 



u: = ^^^(r-^' . (7.9) 

p + p V p + pj 

This is quite a remarkable result. It shows that the relaxation time r is essential in 
order for the system to be stable. We need the imaginary part of the frequency w to 
be positive. Thus, given that the chemical potential, the particle number density and 
absolute equilibrium temperature are positive quantities, and assuming that the weak 
energy conditiorjj holds, we have 

'^'^ > and — — > 0. (7.10) 



p+p P+P 

Therefore, the stability criterion of the system is given in terms of the relaxation time 

through 

kT 

T> . 7.11 

P + P 

Thus, we arrive at a similar result to the one obtained by Hiscock and Olson [57] regard- 
ing the stability of first-order theories. In particular, the Eckart model - for which the 
relaxation time is identically zero - is inherently unstable. 



^For matter described by the energy momentum tensor Tab, the weak energy condition asserts: that for 
every future directed timelike vector field X, the energy density p — X°'X^Tab > 0. In particular, for a 
perfect fluid, it implies p + p > 0. 
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7.2 Longitudinal modes: causality and the onset of instability 

For the case of longitudinal perturbations, we take the three-divergence of the perturbed 



momenta, equations (7.6), to get 



d d 

ot ot 



0, 






+ VM' + i^S^T + KTVi5u' = 0. 



(7.12a) 
(7.12b) 



Here, we also need to consider the divergence of the perturbed conservation laws. These 
take the form 






1 a^ 

n dt 



5n, 



-nT^s, 



(7.13a) 
(7.13b) 



where we have defined the specific entropy as s = s/n. This problem is clearly more 
involved than the transverse one. To make progress we need to make some decisions. 
First, we need to choose an appropriate set of variables to work with. Secondly, we need 
to provide an equation of state for matter. In the following we will opt to work with the 
perturbed densities 5n and 5s. In this case, we have the set of thermodynamic relations 



5fi 
6T 



dn 


s 9s , 


dT 

dn 


dT 

s ds 



6s, 
6s. 



(7.14a) 
(7.14b) 



It is also useful to keep in mind that the temperature can be interpreted as the 'entropy 
chemical potential', i.e. it is defined by 



T 



dp 
ds 



(7.15) 



where p = p{n, s) represents the equation of state. This immediately implies the other 

thermodynamic relation 

dT _ dp 

dn „ ds 



(7.16) 



This relation allows us to reduce the number of unspecified thermodynamic quantities, 

which is the key reason for the choice to work with the perturbed densities. It is also 

worth noting that 

dp 



dp 
ds 



n 



ds 



(7.17) 
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Again, we will consider the plane wave solutions such that the perturbations behave 
as exp(iwt + kjX^). In a similar way as in the previous sections, we introduce the phase 
velocity of the waves a = uj/k, where fc^ = kik^. Thus, the above relations lead to the 
coupled equations 



n 
T 

n 



2 n dfj. 



a 



H dn 
ndT 
Tiki 



6n + naT 
6n + htT 



a 



1 dT' 

aT dn 


6s = 


1 


K dT' 
tT ds 



5s = 0. 



(7.18a) 
(7.18b) 



The dispersion relation we will obtain from (7.18) will be a quartic on a. Furthermore, 



we immediately notice the appearance of a complex root. To facilitate the interpretation 



of the roots, we complete the square in the last bracket of (7.18b) 

2 



1 



K dT 



a 



kr tT ds 



1 



a 



n 



K dT 



1 



'2TkJ \tT ds 4/i;2r2 
Let us introduce the speed of sound, specific heat and the coefficient b given by 
2_n5/i 1 _ 1 dT ,,_1 



(7.19) 



(7.20) 



together with the auxiliary 'velocities' 



2 _ndT 



4 



1 dT 

aT dn ' 



2 2 

C3 = C, 



4 



K 



4P^2 • 



(7.21a) 

(7.21b) 
(7.21c) 
(7.21d) 



With the above definitions, the system of equations (7.18) takes the more manageable 
form 



- [a^ - 4] Sn + naT [a^ - cj] 6s = 
T 



6s = 0, 



K— \a'^ — c?l 6n + nrT (a — ib) — c\ 
n '- -' L 

where we can read off the dispersion relation as 

^[(a2-c?)(a2-ci)]-[a2-ci][(a-^6)2-cI 



0. 



(7.22) 
(7.23) 

(7.24) 
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Non-adiabatic perturbations 



Here, we have introduced the additional abbreviation 

A = . 



IIT 



Expanding (7.24) gives us the quartic polynomial 

(1 - A)a^ + iBa^ + Ca^ + iDa + E = 0, 
where each of the coefficients are given by 

B = 2b, 

C = b^ + cl + 4-A{cl + 4) , 



D 



-2bc 



3' 



E = Acl4 - clb^ 



2 2 
C3C4. 



(7.25) 



(7.26) 



(7.27a) 
(7.27b) 
(7.27c) 
(7.27d) 



This expression is obviously too complicated for us to be able to make definite statements 
about the solution, without making further assumptions. The most direct strategy would 
be to consider an explicit example equation of state, work out the relevant thermodynamics 



quantities, and then solve (7.26). This would allow us to establish whether the considered 



model is stable and causal. This route is, however, not particularly attractive given the need 
to introduce an explicit model. Nevertheless, we still have hope to make some progress 
by attempting the alternative involution method to find some stability criterion in the 
following manner. Suppose rj, i = 0,1,2,3 are the roots of (7.26). If we now transform 
this set of solutions by an involution into 



•So = 2 ('^o + n +r2 + r3), 
si = 2 ('^o - ri + r2 - rs) , 
•S2 = 2 ('^o + r-i - r2 - rs) , 



S3 



(ro - ri - r2 + ra) 



where we know from (7.26) that 



So 



B 



(7.28a) 
(7.28b) 
(7.28c) 
(7.28d) 

(7.29) 
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Thus, we can evaluate so directly in terms of known thermodynamic quantities 
3 n , • r ,^ -1-1 



k=0 



B 



2(1-^) A- I 



2^ 



akT 



1 



(7.30) 



which, by substituting the definition of a given by equation (7.2), reduces to 



So 



2k 



nfl [{k — TS)T — THfl] 



-1 



(7.31) 



Writing sq as the sum of the real and imaginary parts of the roots of (7.26) 

3 3 



1 



So 



y^ Re(rfe) + i ^ Im(rfc) 



.A;=0 



fc=0 



it follows that 



^Re(rfc) = 



(7.32) 



(7.33) 



fc=0 



Therefore, there would be no hope for the system to be stable if the sum of the imaginary 
parts is negative, 



^lm(rfc) <0. 



(7.34) 



k=0 

However, this condition requires that (k — ts)T — rrifi > 0, that is, 

> T 



P + P 



(7.35) 



for arbitrary small values of the temperature. Thus, the system is absolutely unstable only 
in the case in which the relaxation time vanishes. This is not a new result but it shows 



that our longitudinal analysis is consistent with the transverse one [see equation (7.11)]. 



If we want to continue to consider the problem in (at least to some extent) generality, 
then we need to resort to approximations. As we will see, this is a very instructive route. 
Thus, let us consider the case for which sq vanishes. This corresponds to the limit of very 



large wave number. In this limit, we can see that (7.26) reduces to the bi-quadratic 



{1 - A)a'^ + Ca^ + E = 



whose roots are given by 



4 



1 

2{l^AJ 



\^-C ± [C^ - iil - A)Ef^] . 



(7.36) 



(7.37) 
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Let us examine this expression carefully. First, we evaluate the coefficient 

okT 1 



1-A = 1 



[{p + p)t-kT] >0. 



(7.38) 



Thus, to guarantee the stability of the system (a'^ > 0), the following conditions must hold 



C 



E 



K 


2 K dT / 
, fiT dn \ 


TS\ 

K ) 


>0 


K 

T 


In /dry cl 

T fi\dnj Cv 


<0. 





(7.39a) 
(7.39b) 



This last condition leads immediately to 



^ ^v ^ary 



^T \dn J 



(7.40) 



It is worth noting that this last condition is readily satisfied by all separable equations 
of state of form p = f{n) + g{s). 



Our analysis leads to two further conditions that must be satisfied. Condition (7.39a) 
leads to 



T > K 



Cv on 



dT 2 



(7.41) 



which, provides us with a further lower bound for the relaxation time. Thus, we have a 
second criterion to establish local stability in a given matter model whose equation of state 
is separable. 

Finally, we can approximate the finite wavelength regime through an expansion in 
powers of e = k~^ . Thus, consider the expression 



1 
a = a± + -cJi = a± + eai, 



(7.42) 



where a± solves (7.36). We can evaluate (7.42) up to the fourth power and neglect terms 



of order O(e^). Thus, we have 



a^ = (a± + eai)^ = 4 + Aea^ai + ©(e^), 
a^ = {a± + eaif = o-| + 2ea±ai + 0{e^). 



(7.43a) 
(7.43b) 
(7.43c) 
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Substituting (7.43) into (7.26) while keeping only linear terms in e gives 



4(1 - A)alai + -al + 2Cai - -c^ = 0, (7.44) 

r r 



which leads immediately to 



I cl-al 



-' = 27 2(1 -Vl%C - <^-^=) 

Since all quantities in this expression are constrained to be real and stability requires 
Im((Ti) > 0, we have two cases: 

1. 2(1 - A)al + C = + [B^-4{1- A)E] ^^^>0^cl< 0%. 

2. 2(1 - A)a'i + C = -[B'^-A{1- A)E] ^^^ < ^ cl > a^_. 

Therefore, we can write the conditions for the stability of our complete first-order theory 
as 

a! < c2 < 4, (7.46) 

which is consistent with the notion that mode-mergers signal the onset of instability. 

7.3 What have we learned? 

In this section, we constructed a consistent first-order expansion away from thermal equilib- 
rium. We argued in chapter 4 that first order theories are 'pathological' due to instabilities 
and causality violations. Therefore, two remarks are in order: 

1. On the one hand, the presence of instabilities is not sufficient to rule out a particular 
matter model. There are physical instabilities, as the two-stream case, which provide 
essential information about the dynamics of the system. Moreover, in situations far 
away from thermal equilibrium, thermal instabilities are responsible for structure 
formation, e.g. Benard cells, Taylor instability, Kelvin-Helmholtz, etc... 

2. On the other hand, causality violations may lead to unstable behaviour. However, 
even when no instabilities are present, in any relativistic matter model this is a clear 
violation of our mathematical construction. Therefore, causality is a requirement to 
consider a matter model as physical. 

Here, we have shown that, in contrast with the conclusion of Hiscock and Lindblom 
[35], a first-order expansion away from thermal equilibrium is possible provided we do not 
exclude second-order terms prematurely. This is manifest in our first stability criterion for 



transverse modes, equation (7.11), where we can see that the Eckart model not only fails 
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to be stable, but also causal. That is because from the construction, the relaxation time 
- which is a second-order effect - was absent from the dynamics of the model. Moreover, 



the stability criterion (7.11) gives a non-trivial lower bound for the relaxation time of 



transverse [and longitudinal waves, see (7.35)] modes and, hence, an upper bound for the 



speed of propagation for thermal disturbances. 

To tackle the problem of stability and causality in an analytic manner for arbitrary 
wave numbers is a very difficult task. However, we can obtain approximate criteria from 
an expansion away from the 'geometrical optics' regime, where the dispersion relation 
reduces to a bi-quadratic providing the upper and lower bounds for the phase velocity. In 
this sense, we obtained the onset of instability as a merger of modes in our expansion away 
from thermal equilibrium. Notice that our only true criterion to rule out the instabilities 
as non-physical is well established. That is, if the onset of instability occurs beyond the 
causality regime, every unstable mode will necessarily be non-physical. However, if the 
instability region is within the causality domain, we cannot rule out the instability on 



the grounds of local causality. Moreover, (7.46) gives us the time-scale for the unstable 



behaviour, providing us with the relevant information to observe their evolution given a 
particular matter model, e.g. superfluid states of matter inside a neutron star or, from the 
non-adiabatic pressure perturbations in a cosmological setting [19] . 

This is just the first step towards a better understanding of the dynamical basis of 
dissipative systems in a general relativistic setting. There are still many issues to address. 
Of immediate interest is to apply our stability analysis to specific matter models to test 
and predict the dynamics of heat transport in relativistic systems. 
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Closing remarks 



The central aim of the research presented in this thesis was to elaborate on the geometric 
formalism pioneered by Carter with the aim of understanding heat flow in a large range 
of physical scales and situations. We formulated a theory of relativistic heat conduction 
through a completely generally covariant constrained variational principle for a multi- 
fluid system. In such an approach, we dealt only with the matter sector of the Einstein- 
Hilbert action for a system consisting of two coupled fluids - particles and entropy - whose 
interaction is a crucial ingredient to guarantee stability and causality in situations away 
from thermodynamic equilibrium. 

We started our discussion by introducing the basic empirical principles upon which we 
construct general relativity - the equivalence principle and general covariance - and the 
way they lead to a dynamical theory of geometry and matter. We also discussed that the 
laws of thermodynamics are a collection of self-consistent empirical facts which must be 
added as additional hypotheses to our fundamental theory of motion. In this sense, these 
empirical laws rule out entire classes of phenomena predicted by the dynamics but whose 
occurrence in nature may never be observed, e.g. the flow of heat from colder to warmer 
bodies. This constitutes one of the central topics of fundamental debate on the apparent 
time-asymmetry of nature. In this work we have not entered this debate but it would have 
been remiss not to acknowledge the conceptual difficulties inherent in the inclusion of the 
second law of thermodynamics into any dynamical framework. 

Variational techniques are powerful and elegant mathematical tools. We have formu- 
lated a theory of relativistic heat conduction through a completely generally covariant 
constrained variational principle for a multi-fluid system. In such an approach, we dealt 
only with the matter sector of the Einstein-Hilbert action for a system consisting of two 
coupled fluids - particles and entropy - whose interaction is a crucial ingredient to guarantee 
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stability and causality in situations away from thermodynamic equilibrium. 

The motivation for treating entropy as another 'fluid' coupled to the matter current - 
rather than considering it merely as a component - stemmed from the hydrodynamic limit 
of relativistic kinetic theory. The emergence of an entropy density current, generally not 
aligned with the matter flow, is generic to every theory of relativistic dissipation. The 
novel feature in our programme is that we incorporate the entropy flux into the action 
principle, allowing us to truly analyse the dynamics of the coupled matter and entropy 
system, where the entrainment effect plays a crucial role. 

Compared to the standard single- fluid analysis, a multi- fluid system has more degrees 
of freedom; we need to account for relative flows between them. This is an essential re- 
quirement for the two-stream instability. It is a generic phenomenon that does not require 
particular fine-tuning to be triggered, nor is limited to any specific physical system. The 
only requirements are that there must be a relative (background) flow and a coupling 
between the fluids. Such instabilities are known to exist for a variety of configurations. 
For example, in shearing motion at an interface, this corresponds to the well known Kelvin- 
Helmholtz instability. A key result from our work was a covariant generalization of the 
two-stream instability analysis to relativistic situations. Here, we have shown that the equi- 
librium configurations of our variational construction, as expected, are absolutely stable. 
Moreover, we argued that, for some equations of state, the onset of instability is arbitrarily 
close to the equilibrium regime. It is worth noting that this instability for thermodynamic 
systems is a prediction of our programme for relativistic heat conduction. The dynamical 
role of this instability has not been explored for particular matter models beyond the linear 
regime. 

Finally, we took our first step in non-adiabatic equilibrium perturbations at the linear 
regime. We obtain the criteria for the stability of this type of perturbations. It follows from 
our analysis that the Eckart model is intrinsically unstable because it prematurely ignores 
higher order terms in deviations from equilibrium. In this sense, our results constitute a 
consistent 'first-order' theory of relativistic dissipation. 

8.1 Further developments 

There are a number of future directions which our formalism will allow us to explore. To 
list a few, we have: 

8.1.1 Microscopic foundations of relativistic thermal dynamics 

There are many directions that can be explored regarding the connection between the 
microscopic canonical dynamics of a relativistic system and its averaged thermal dynamics 
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on spacetime. Here, we enlist three steps of a particular strategy: 

1. Kinetic limit of the two-fluid formalism. 

In order to improve our understanding of the variational model in a physical setting, 
and the role of the various non-equilibrium terms, we need to consider our heat 
model in the context of relativistic kinetic theory. The aim should be to develop a 
consistent model for the entrainment and the thermal relaxation of a dissipative fluid. 
A detailed exploration of the constrains on the relativistic single-particle phase-space 
imposed by the two-fluid nature of the macroscopic system is needed. This would 
help us to understand the connection between the variational formalism and the 
geometry of the canonical formulation for the microscopic dynamics of a dissipative 
system. Let us note that the standard approach to relativistic kinetic theory is 
based on two fundamental assumptions: the notions of a dilute gas and the (time- 
asymmetric) molecular chaos hypothesis. It is relevant to ask whether these concepts 
need to be extended in order to match more complex macro-scale physics. After all, 
the variational approach is readily extended to more complicated multi-component 
systems {e.g. with varying composition). It also accounts for non-linear phenomena 
that would require the inclusion of higher moments of the distribution function at 
the microscopic level. 

Thus, our formalism could allow us to address a range of conceptual issues concerning 
the entropy of a relativistic system. Most importantly we may be able to develop a 
link between the two-fluid phenomenology and statistical physics. An understanding 
of this issue is important for a number of relevant applications. The key point 
is that, from a technical point- of-view, the variational model can account for a 
number of 'distinct' entropy components e.g. phonons, rotons, etc... This may be an 
advantage in order to describe specific systems, like relativistic superfluids at finite 
temperatures. However, it is not clear how these models are linked to statistical 
physics. To explore this problem, one can extend existing models for phonons in 
finite temperature superfluids (and the associated phonon hydrodynamics) to more 
complex setting with several 'distinct' entropy-carrying components. 

2. The matter-space dynamics and the geometry of dissipation. 

The notion of individual matter spaces associated with the various fluid components 
played a central role in the variational construction of chapter 5. In this case, the 
matter-space construct essentially boils down to the variation of flow lines associated 
to the individual fluid elements. This problem is relatively simple because the matter 
space has no real physical 'structure'. However, there are situations where it is 
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necessary to generalize this idea. A particular case concerns elastic matter |13l [39] , 
where the elasticity is encoded in the matter-space geometry. This problem is fairly 
well understood. What is less clear at the present time is how we should handle 
evolving system, e.g. when elastic matter undergoes plastic flow or when there are 
composition changes due to reactions in a multi-component system. This problem 
is central to the description of irreversible phenomena. After all, the second law 
requires us to allow for entropy production. In order to do this, as in the heat-flow 
model described above, we need to relax the variational assumptions. 

A natural step forward in this direction would be to develop the variational formalism 
further by considering 'evolving' matter spaces. This idea is technically challenging, 
but we see no reason why one should not be able to make progress on it. We will 
need to add more structure to the description (as in the elastic case) and likely 
account for interaction between individual matter spaces (particularly if there are 
reactions involved), but the mathematics should be sufficiently flexible to deal with 
these issues. The main challenge may be conceptual. At the deeper level we need 
to develop a link between the matter-space and the coarse-graining of phase-space 
in statistical physics. It is, in fact, relevant to ask at what level the statistics should 
be considered. Is it at the spacetime level, or is it in the lower-dimensional matter- 
space? In principle, both answers seem viable but the latter would be an attractive 
(possibly quite revolutionary) solution. In analogy with the description of elastic 
matter one may envisage a model based on the notion of an evolving 'thermal geo- 
metry' (in matter-space) directly linked to the entropy change between hypersurfaces 
in spacetime. The model also requires a dynamical map between the matter space 
and spacetime, in order to link the changes in the local geometric structure of the 
matter configuration (described in terms of normal coordinates) to the macroscopic 
evolution of the system. 

3. The thermodynamics of the gravitational degrees of freedom 

This direction is, perhaps, one of the most challenging from a foundational point-of- 
view. It stems from the fact that there seems to be an intimate link between gravity 
and thermodynamics. This follows immediately from the equivalence of energy and 
mass, which implies that heat must affect the gravitational field. Nevertheless, from 
a conceptual point of view it is not understood to what extent the gravitational field 
is 'hot'. Basically, we do not yet have a clear definition of the entropy associated 
with an evolving gravitational field. 

The variational model for heat accounts for the coupling to (and evolution of) the 
gravitational field via the Einstein field equations. However, so far our main focus 
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has been on the matter sector of the problem. In order to explore the role of the 
gravitational field one can consider the coupled system of equations in a number of 
model scenarios (plane waves and simplistic matter models). The aim is to established 
how the second law of thermodynamics feeds into the evolution of the gravitational 
field, and (conversely) how variations in the gravitational field affect the entropy and 
the heat flow. An interesting conceptual issue concerns the link between observers 
and the increase of gravitational entropy. 

8.1.2 The 'near' and 'far' from equilibrium entropy conditions 

We have argued in chapter 2 that, although the evolution towards equilibrium of a ther- 
modynamic system is governed by the second law, there is no criterion to 'select the path' 
towards equilibrium, i.e. there is no fundamental principle of extremal entropy production. 
It would be interesting to explore how systems thermalise once the entropy production is 
constrained. The variational model provides us with a useful context for such work. How- 
ever, we would need to develop a clearer understanding of the entropy production associated 
with a given initial configuration. 

To make progress in this direction, one can consider two extreme cases; the near and 
far from equilibrium regimes. In the former case, we can rely on an expansion away from 
thermal equilibrium (as in the moment expansion of kinetic theory), while in the latter 
we need to account for the non-linearities of the variational model. Non-linear systems 
are, in fact, very interesting because the available models for relativistic heat differ in 
the quadratic (and higher order) terms. Thus, we need to explore to what extent these 
differences lead to distinguishable differences in the evolution. To this end, one can build 
on experience of non-linear evolutions of relativistic material systems in general relativity, 
accounting for a suitably constrained thermodynamics. 

1. Near-equilibrium regime. In chapter 2, we gave an example of the principle of min- 
imum entropy production for systems close to thermal equilibrium. However, this 
notion is closely related to the Onsager's symmetry principle of the Linear Irreversible 
Thermodynamics programme. It is, of course, not at all obvious that usual derivation 
can be extended to the relativistic problem. However, if the argument fails then we 
will have learned something fundamental about non-equilibrium physics. 

An understanding of the near-equilibrium behaviour gained from a perturbative ap- 
proach would obviously be incomplete. Moreover, as we have shown at the end of 
chapter 6, even at linear order in a perturbation expansion we need information about 
the full non-linear regime. Nevertheless, such analysis gives us the causality criterion 
and the onset of thermodynamic instability. In the case of physical unstable modes, 
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these criteria would be serve to set limit cases for the evolutions in the non-linear 
regime. 

2. Far- from- equilibrium regime. To fully understand the role of the second law and its 
association with pattern and structure formation we need a non-linear framework 
which must account for (as in numerical relativity) the coupling between matter and 
spacetime. In order to explore the non-linear sector, we may consider the maximum 
entropy principle. This idea is based on the notion that when estimating a prob- 
ability distribution, it is natural to select the distribution that leaves the largest 
remaining uncertainty (maximum entropy) consistent with the system constraints. 
The advantage of this approach is that it minimises the information required (which 
indeed was Jaynes' motivation for the approach) in the calculation. The maximum 
entropy principle is intimately linked to information theory, identifying the entropy 
in that context with the statistical physics. This connection is relevant for the grav- 
ity problem, especially as it may provide insight into the issues like information loss 
associated with black-hole horizons. 

To sum up, one could extend the variational model in such a way that the evolution 
ensures maximum entropy production. This can be done using Lagrange multipliers, 
but in order to ensure computational efficiency we need to learn from other contexts 
where this principle has been implemented {e.g. for biological systems). Successful 
implementation of this idea is likely to lead to a real breakthrough in non-linear 
relativity simulations of heat. 

8.1.3 Classical origins of the arrow of time 

The research described in this thesis is entirely within the realm of classical physics. Yet, it 
is clear that any real understanding of the nature of time will have to account for quantum 
aspects. While we make progress at the classical level, we should prepare the ground for 
future explorations of the quantum arena. In absence of a theory for quantum gravity this 
is obviously a huge step, but some basic principles seem clear. For example, if we want to 
make direct contact (and build upon) Rovelli's notion of thermal time (associated with the 
evolution of pre-symplectic systems), then we need to develop a Hamiltonian description 
for the relativistic thermo/hydrodynamics. This is known to be a challenging problem, but 
Dirac's procedure for developing a Hamiltonian system from a given Lagrangian is (at least 
in principle) clearly laid out. However, the steps involved are far from straightforward in 
practice. This is particularly true in the case of a constrained variational model, as in the 
present case. 
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8.2 Conclusion 

The work presented in this thesis represents a few fundamental steps towards a deeper 
understanding of a relativistic theory of irreversible processes. Thus, we have learned that 
the dynamical contribution of heat plays a crucial role in the evolution of a material system 
in spacetime. Moreover, we have seen that, to correctly account for such dynamics, we 
need to include the contribution of heat directly into the variational principle through the 
entropy current. In this manner, the equations of motion for an internal observer (one 
which is moving with the material particles) include a thermodynamic description where 
the inertial properties of heat are an essential feature for the system's evolution towards 
equilibrium. To show this, we have made a linear perturbation analysis for both the 
adiabatic and non-adiabatic cases. In the former, we generalise the two-stream instability 
analysis for linear perturbations of equilibrium to a relativistic setting. The consequences 
of this analysis are yet to be explored. In the latter case, we obtained a consistent first 
order perturbative analysis of deviations from thermal equilibrium. This is a novel and 
unexpected result which provides new light on the old second-order paradigm. However, 
let us emphasise that such a consistent first-order model can only stem from an underlying 
second order dynamical theory which contains the complete thermodynamic information 
for the propagation of linear waves. Finally, we have presented a number of exciting 
directions in which this programme can be extended. With some luck, one of them may 
give us a hint in a quest that, just like entropy itself, grows as we walk ahead. 
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Figure 8.1: Gravity and entropy production: What can our intuition tell us about the 
future of this initial configuration? [This image is property of Steve Strawn, splutphoto.com] 
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Regular and Anomalous Heat Conduction 



This appendix contains an extract of a set of notes by Brandon Carter \l^ . We 
have included this transcript here in order to make more accessible a reference 
which, otherwise, would be a bit less than impossible to ge\j The material be- 
low is a collection of the original ideas which motivated the later developments 
presented in this work. 

The purpose of this work is to offer a fresh approach to an old problem by providing an 
appropriate relativistic formulation of classical thermodynamics. In particular, of the law 
of heat conduction. The essential feature of this approach is the prominence given to the 
canonical 4-momentum covector associated with the entropy flux. The main result of this 
work is to provide a more satisfactory alternative to the traditional Eckart heat conduction 
equation. 



Eckart's equation [c.f. Eckart's hypothesis, equation (4.31)] differs from the non 



relativistic Fourier's equation only by the second term on the left hand side involving 



the acceleration u". Many authors have found that (4.31) is unsatisfactory because it 
leads to a causality paradox (heat conduction faster than light), and it has been suggested, 
notably by Cattaneo [18] and Kranys |40j that this should be remedied by replacing the 
LHS of Eckart's equation by the ad-hoc expression q°'-\-Th'^^q^ , where the coefficient r of the 
additional term is a new phenomenological parameter, to be interpreted as representing a 
microscopic relaxation time. 

In the present approach, based purely on the most intrinsically plausible thermody- 
namic axioms, we find that an extra term having essentially the same effect as one of 



the Kranys type turns up automatically on the LHS of (4.31). Part of our purpose is 



^The author is thankful to Greg Comer for this communication. 
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to point out that it is not surprising that the original form of (4.31) gave unsatisfactory 
results. The reason for this, is that the traditional textbook derivation of Eckart's conduc- 
tion equation is based on a rather incoherent set of simplifying assumptions. The weak 



point in the traditional derivation of (4.31) is the (usually implicit rather than explicit) 
assumption that inertia associated with the heat current is zero. A priori, this might not 
seem unreasonable in a course-tuned treatment aimed at maximum simplicity, but its utter 
inappropriateness is made evident post-facto by the fact that the Eckart acceleration term 
in is naturally interpretable (c./. Ehlers |26| ) as being due to an effective inertia which 
amounts quantitatively (in units for which Boltzmanns constant is unity) to the rest mass 
associated with the energy per unit entropy. It follows that instead of ignoring the inertia 
associated with the heat current, a much more reasonable ansatz would be to suppose 
consistently from the outset that there is an effective mass per unit entropy given precisely 
by the temperature 9. When this condition is satisfied, we shall describe the heat flux as 
regular. In the regular case, we shall show that an appropriately modified version of the 



standard textbook derivation of (4.31 ) leads unambiguously to the appearance of an extra 
term C^{q"'/s) on the LHS, where s is the rest frame entropy density, and £u denotes Lie 
differentiation with respect to the fundamental flow vector u'^. 

Although it is certainly not fair to suggest that this regularity condition is 'ad-hoc' (as 
was implied by Israel with reference to an abbreviated oral communication of this work) it 
would nevertheless be entirely justifiable to describe it as 'naive', since it is certainly not a 
logical necessity that inertia appearing in different contexts should always be strictly the 
same. In defence of the regularity condition all that one can say is the following. 

Firstly, it will be accurately valid in the limit when thermal energy is entirely in the 
form of a perfect gas of black body equilibrium radiation with negligible coupling to the 
background medium, and hence it will be a good approximation in many problems in- 
volving radiative transfer in radiation pressure dominated stellar or cosmological media. 
Secondly, in many other contexts the difference between the various conceivable correction 
terms will be important only for variations of such high frequency that a straightforward 
thermodynamic treatment would fail anyway, so that one is in practice free to choose the 
most mathematically satisfactory hypothesis, i.e. that of regularity. Nevertheless, in addi- 
tion to these two kinds of physical situation (i.e. those where regularity holds accurately, 
and those where there is no practical loss of accuracy in imposing it for convenience) there 
are also situations where anomalous behaviour (i.e. deviations from regularity) can give 
important measurable effects, the most obvious example being that of superfluidity. One 
may draw a close analogy with the case of regular and anomalous magnetic moments: 
for the electron, the 'naive' Dirac's theory applied very well. However, for more strongly 
interacting particles, such as the proton, a significant anomaly occurs. In a somewhat ana- 
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logous manner, the strong correlations between atoms of superfluid helium make possible 
the occurrence an effective inertia per unit entropy that is enormously greater than the 
regular value (which is of course negligibly small in such a highly non-relativistic context). 
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